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Abstract 



In this thesis, I study the phase diagram of dense, locally neutral three-flavor quark matter as a 
function of the strange quark mass, the quark chemical potential, and the temperature, employing 
a general nine-parameter ansatz for the gap matrix. At zero temperature and small values of 
the strange quark mass, the ground state of quark matter corresponds to the color-flavor-locked 
(CFL) phase. At some critical value of the strange quark mass, this is replaced by the recently 
proposed gapless CFL (gCFL) phase. I also find several other phases, for instance, a metallic 
CFL (mCFL) phase, a so-called uSC phase where all colors of up quarks are paired, as well as the 
standard two-flavor color-superconducting (2SC) phase and the gapless 2SC (g2SC) phase. 

I also study the phase diagram of dense, locally neutral three-flavor quark matter within the 
framework of a Nambu-Jona-Lasinio (NJL) model. In the analysis, dynamically generated quark 
masses are taken into account self-consistently. The phase diagram in the plane of temperature 
and quark chemical potential is presented. The results for two qualitatively different regimes, 
intermediate and strong diquark coupling strength, are presented. It is shown that the role of 
gapless phases diminishes with increasing diquark coupling strength. 

In addition, I study the effect of neutrino trapping on the phase diagram of dense, locally 
neutral three-flavor quark matter within the same NJL model. The phase diagrams in the plane 
of temperature and quark chemical potential, as well as in the plane of temperature and lepton- 
numbcr chemical potential arc presented. I show that neutrino trapping favors two-flavor color 
superconductivity and disfavors the color-flavor-locked phase at intermediate densities of matter. 
At the same time, the location of the critical line separating the two-flavor color-superconducting 
phase and the normal phase of quark matter is little affected by the presence of neutrinos. The 
implications of these results for the evolution of protoneutron stars are briefly discussed. 
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Chapter 1 

Introduction 



The phase diagram of neutral quark matter was poorly understood as I began with the research 
on this topic in 2003. The task of my thesis was therefore to illuminate the phase structure of 
neutral quark matter. A phase diagram is a two-dimensional diagram with axes representing the 
temperature and the chemical potential, the density, or other similar quantities. Therefore, phase 
diagrams tell us in which state is a system for a given temperature and a given chemical poten- 
tial. Besides, they contain the information at which temperatures and which chemical potentials 
transitions to other phases occur. Such phase transitions can be of first or second order, or simply 
crossovers. This depends on the order parameter of the system. If it changes discontinuously, then 
a first, otherwise a second-order phase transition or a crossover appears. 

In Sec. 11.11 I show the status of knowledge of the phase diagram of strongly interacting matter 
before I began with my research for this thesis in 2003. As one can see, the phase diagram of neutral 
quark matter was indeed poorly understood at that time. In Sec. 11.21 I discuss the behavior of 
sufficiently cold and dense quark matter, namely that quark matter is color-superconducting and 
that color-superconducting quark matter appears in several phases. The most important color- 
superconducting phases are presented in Sec. 11.21 In Sec. 11.31 I give a short but comprehensive 
introduction into stellar evolution because the cores of neutron stars are the only places in nature 
where one expects neutral color-superconducting quark matter. Thereby, it will become clear 
what a star is, how it is formed, which processes happen in a star, and what are the final stages 
of stellar evolution. In Sec. 11.41 I focus on neutron stars. It is explained how neutron stars are 
formed, of which matter they consist, and I present the structure of neutron stars. In Sec. 11.51 
it is argued why neutral color-superconducting quark matter is expected to occur in neutron star 
cores. 

In Chapter^ I present the calculations to obtain the pressure for neutral color-superconducting 
quark matter, and I show the phase diagram of neutral quark matter. With this one can predict 
in which state neutral quark matter is in the cores of neutron stars. 

In Chapter [3J I summarize the results and conclude my thesis. 

Important definitions and useful formulae can be found in the Appendix. 

1.1 The phase diagram of strongly interacting matter 

The fundamental theory of the strong interaction is called quantum chromodynamics (QCD). 
The participants of the strong interaction are the elements of hadronic matter, the quarks and 
gluons. Quarks interact via gluons because both particle species carry so-called color charges (red, 
green, and blue) which are responsible for the interaction. In our everyday life, we do not see 
quarks or gluons because they are confined into hadrons. This is because the quark interaction 
caused by gluons is so strong that they cannot exist as free particles. QCD is an asymptotically 
free theory [1]. At high temperatures or densities, the quarks are deconfined because their mutual 
distances decrease and the exchanged momenta increase so that the interaction becomes sufficiently 
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weak [2]. The state of deconfined quarks is called the quark-gluon plasma (QGP). Such a phase 
certainly existed in the early universe which was very hot, but close to net-baryon free. Nowadays, 
the only place in nature where a QGP may exist is in the interior of neutron stars. Here, the density 
is extremely high and the temperature low. The third place where an artificially created QGP 
could appear is in heavy-ion collisions. The temperatures and densities which are reached by the 
collisions depend on the bombarding energies. 
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Figure 1.1: The status of knowledge of the phase diagram of strongly interacting matter as I began with 
the research on this topic for my thesis in 2003 [3]. Note that this phase diagram is drawn 
schematically. 



In Fig. 11.11 I show the status of knowledge of the phase diagram of strongly interacting matter 
in the plane of temperature T and quark chemical potential before I began with the research for 
my thesis in 2003 [3] . There is a phase transition at the point (T, /i) ~ (0, 308) MeV which separates 
the gaseous nuclear phase at lower a from the liquid nuclear phase at higher \i. The nuclear liquid- 
gas transition [4] which is a first-order phase transition starts from this point and disappears in 
a critical endpoint at T ~ 10 MeV and slightly lower quark chemical potential. In this endpoint, 
the transition is of second order. Above the endpoint, there is no distinction between these two 
phases. The transition point at zero temperature is easily found because nucleons of infinite and 
isospin-symmetric nuclear matter in the ground state at normal nuclear density no ~ 0.15 fm 
are bound by 16 MeV (if one neglects the repulsive Coulomb forces). The energy per baryon or 
the energy density per baryon density respectively is given by E/Nb = e/ns = rn n — 16 MeV 
where m n ~ 940 MeV is the rest mass of the baryon. By using the thermodynamic relation 
e = Ts — p + nn, one obtains for the ground state of nuclear matter where the pressure p = 
that the baryon chemical potential is identical to the energy per baryon [1b = £/n>B — 924 MeV. 
Since a baryon contains three quarks, the quark chemical potential is one third of the baryon 
chemical potential, fi = /is/3. This leads to the result fi ~ 308 MeV for the transition point at 
zero temperature. 

Nuclear matter consists of droplets which is the most energetically preferred form for nuclear 
matter in the low-density and low-temperature regime of the phase diagram. In the gaseous phase 
at low density and nonzero temperature, nucleons will be evaporated from the surface of the 
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droplets so that there is a mixture of droplets and nucleons. As soon as the chemical potential 
exceeds the one corresponding to the nuclear liquid-gas phase transition, only droplets of nuclear 
matter but no evaporated nucleons will appear. 

At low quark chemical potentials, strongly interacting matter is in the hadronic phase. At 
nonzero temperature, nuclear matter not only consists of nucleons but also of thermally excited 
hadrons. Therefore, at low quark chemical potentials, a large amount of pions can be found. By 
increasing the temperature, the system passes the quark-hadron transition line and enters the 
regime of the QGP. The critical endpoint of the transition line at (T, fi) = (162 ± 2, 120 ± 14) MeV 
obtained by lattice QCD calculations [5] depends on the value of the quark masses and is of second 
order. For smaller quark chemical potentials, the transition becomes a crossover, and there is no 
real distinction between hadronic matter and the QGP. For larger quark chemical potentials, the 
transition is a line of first-order phase transitions that separates the hadronic phase from the QGP. 
But it is not known whether the line of first-order phase transitions goes all the way down to zero 
temperature. If so, also the precise value of the quark chemical potential for the phase transition 
at zero temperature is unknown. One should mention that these lattice QCD calculations are 
not very reliable at nonzero quark chemical potential. In addition, these calculations are done 
with probably unrealistic large quark masses and on fairly small lattice sizes. For smaller quark 
masses, the endpoint should move towards the temperature axis. By increasing the quark chemical 
potential, the nucleons will be packed denser and denser until the quark-hadron phase transition is 
reached. As soon as the quark chemical potential exceeds this transition line, the system becomes 
a color superconductor at low temperatures or a QGP which is in the normal conducting phase 
(NQ) at high temperatures. 

At large quark chemical potentials and low temperatures, quark matter becomes a color super- 
conductor. In this thesis, I shall focus on this part of the phase diagram and its transitions to the 
hadronic phase and to normal quark matter. The reader will see that there exist various color- 
superconducting phases and that this thesis will update the phase diagram of strongly interacting 
matter in the quark regime. 



1.2 Color superconductivity 

Quarks are spin-i fermions and therefore obey the Pauli principle which requires that one quantum 
state is occupied by only one fcrmion. At zero temperature, non-interacting quarks occupy all 
available quantum states with lowest possible energies. This behavior is expressed with the Fermi- 
Dirac distribution function for zero temperature, 

f F (k) = 6(f,-E k ) , (1.1) 

where E k — \^k 2 + rn 2 is the energy of a free massive quark. All states with momenta k = |k| 
which are less than the Fermi momentum kp = \/ [i 2 — m 2 are occupied. The states with momenta 
larger than the Fermi momentum are empty. The pressure for massive non-interacting quarks at 
zero temperature is given by 

n / >oc ' k 4 n f k " k 4 

P=A/ dk ^rfF( k )- B =7^ dk—-B, (1.2) 



67T 2 J E k 67T 2 J ' Ey, 

where g = 2N c Nf is the degeneracy factor in which N c is the number of colors, and Nf the number 
of flavors. The factor two comes because of spin degeneracy. The bag constant B assigns a nonzero 
contribution to the vacuum pressure and, in this way, provides the simplest modelling of quark 
confinement in QCD [6]. A typical value for the bag pressure is B x l 4 = 0.17 GeV which is used in 
Rcfs. [7,8]. In the limit of zero mass and zero temperature, the pressure of non-interacting quarks 
reads 

p-dpS-a- (") 

This is the case of non-interacting quarks. What happens when the quark interaction is switched 
on? 
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At asymptotically large quark chemical potentials, the strong coupling constant becomes small 
so that the dominant interaction between quarks is given by single-gluon exchange. The quark- 
quark scattering amplitude in the one-gluon exchange approximation is proportional to 



ErpArpA _ 
J-ki 1 !] - 

A=l 



- W1 -(f ifc fc- 6^ + ^1(6^-5^), 



(1.4) 



where i, j are the colors of the incoming, and fc, I those of the outgoing channel. The first term in 
this equation is antisymmetric and corresponds to the antitriplet channel which is responsible for 
the dominant attractive interaction while the second term is symmetric and corresponds to the 
sextet channel which is responsible for the repulsive interaction, see Fig. 11.21 




+ 



a ' ~s 

attractive repulsive 



Figure 1.2: The diagrammatic representation of the one-gluon exchange interaction between two quarks 
in QCD [9]. 



Therefore, the antitriplet channel with its dominant attractive interaction ensures that quarks 
with large momenta (quarks near the Fermi surface) form bosonic quark Cooper pairs [10-12]. 
This state is called a color superconductor in analogy to superconductivity of electrons [13]. The 
arguments how color superconductivity is created hold rigorously at asymptotically large densities. 
The highest densities of nuclear matter that can be achieved in colliders or that occur in nature 
in the cores of neutron stars are of the order of ten times the nuclear matter ground state density 
at which the quark chemical potential approximately amounts to \i ~ 500 MeV. Nevertheless, 
calculations in the framework of an NJL model [14] show that color superconductivity also occurs 
at moderate densities and is not limited to asymptotically large densities [15]. 

As a consequence of (color) superconductivity, there exists at least one gap in the quasiparticle 
spectra. Such a color-superconducting gap (parameter) A is a diquark condensate which is defined 
as an expectation value, 

A cx (iFOip) , (1.5) 

where the operator, 

O = CUor <g> ©flavor <g) Chirac , (1-6) 

acts on the quark spinor field i\) in color, flavor, and Dirac space. The color-superconducting 
gap parameters A, are zero in normal quark matter and nonzero (A^ ~ 100 MeV) in color- 
superconducting quark matter, and they are equal to one half of the binding energy of a quark 
Cooper pair. The values of the gap parameters can be obtained by solving the gap equations 

dpcsc = Q (17) 

where pese is the pressure of color-superconducting quark matter. 

In ordinary superconductors, the [£/(l) cm ] gauge symmetry is broken so that the photons 
become massive. (Throughout this thesis, I indicate local, i.e., gauged, symmetries by square 
brackets.) This leads to the so-called Meissner effect, the expulsion of magnetic fields from the 
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superconducting region. In a color superconductor, the [S'6 r (3) c ] color gauge symmetry is broken 
so that some of the eight types of gluons become massive. In an ordinary superconductor as well 
as in a color superconductor, thermal motion will break up Cooper pairs and therefore destroy the 
(color-)superconducting state. In a color superconductor, this transition is of second order and 
happens at the critical temperature 

T BCS = 2l Ao _ o 567Ao ( L g) 
TT 

where Ao is the color-superconducting gap at zero temperature, and 7 ~ 0.577 the Euler- 
Mascheroni constant. 

But there are also differences when one compares superconductivity with color superconduc- 
tivity: the electrons in superconductors first have to overcome their repulsive Coulomb forces in 
order to form Cooper pairs while in color superconductors, the formation of quark Cooper pairs 
is much simpler because there already exists the attractive interaction in the antitriplet channel. 
Quarks, unlike electrons, come in various flavors, see Table fTTTl and carry color charges. Because 
of this latter quark property, superconductivity of quarks is called color superconductivity. 



Flavor 


Mass [McV] 


Q [e] 


up 


2. ..8 


2/3 


down 


5... 15 


-1/3 


strange 


100... 300 


-1/3 


charm 


1300 . . . 1700 


2/3 


bottom 


4700 . . . 5300 


-1/3 


top 


174000 ± 17000 


2/3 



Table 1.1: The masses and electric charges of the quark flavors [16]. The abbreviation for the respective 
quark is written in the bold Roman font. There also exists the corresponding antiquark for 
each flavor with charge —Q. They are not shown here. 

The quarks in color-superconducting quark matter are called quasiparticles or quasiquarks, 
respectively. One only needs to consider the lightest quarks (up, down, and strange) for color- 
superconducting quark matter, because the heavy quarks (charm, bottom, and top) are so massive 
that their occurrence in quark matter is extremely unlikely for the densities and temperatures 
under consideration. 

The color and flavor structure of the condensate of quark Cooper pairs which is also called 
the color-flavor gap matrix depends on which quark flavors participate in pairing and which total 
spin J the Cooper pairs have. For J — 0, the spin part of the quark Cooper pair wavefunction 
is antisymmetric, and therefore the color-flavor part has to be symmetric in order to fulfill the 
requirement of overall antisymmetry. Since quarks pair in the antisymmetric color-antitriplet 
channel, the flavor part of the wavefunction also has to be antisymmetric. 

1.2.1 The 2SC phase 

In order to fulfill the symmetry requirements, at least two quarks of different flavor are needed for 
the J = case. Therefore, the simplest ansatz for the color-flavor gap matrix has the form: 

= ^ke f9 A k , (1.9) 

where the order parameter 

A fc = <5 fe3 A (1.10) 

conventionally points in anti-blue color direction. Together with the Dirac part, the gap matrix 
reads 

*'? = -ysA#. (1.11) 
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This is the so-called 2SC phase which is an abbreviation for 2-flavor color superconductor. The 
color indices i and j run from one to the number of colors iV c which is equal to three. The flavor 
indices / and g run from one to the number of flavors Nf participating in pairing which is equal 
to two in the 2SC phase. In this phase, red up quarks pair with green down quarks, and red down 
quarks pair with green up quarks, and form anti-blue quark Cooper pairs. 

The blue quarks remain unpaired and therefore cause gapless quasiparticles. These quasi- 
particles give dominant contributions to the specific heat and to the electrical and heat conduc- 
tivities. They are also responsible for a large neutrino emissivity produced by the (3 processes 
db — > Ub + e~ + D c and Ub + c~ — > db + v e . The other four quasiparticles and quasiantiparticles 
fulfill the dispersion relation 

'(E k - etf + |A| 2 , (1.12) 

where A is the gap, and e = ± stands for quasiparticles and quasiantiparticles, respectively. At 
small temperatures (T <C A) , the contributions of these quasiparticles to all transport and many 
thermodynamic quantities are suppressed by the exponentially small factor exp (— A/T) [9]. The 
gluons are bosons and therefore, their number density is small at low temperature. In the 2SC 
phase, the [S'C/(3) C ] gauge symmetry is broken to [SU(2) C ]. Consequently, there are 8 — 3 = 5 
broken generators. They represent five gluons which are gapped because of the color Meissner 
mass. Therefore, gluons have only tiny influence on the properties of quark matter in the 2SC 
phase. The unpaired blue quarks are responsible for the absence of baryon superfluidity. Only 
the anti-blue quasiparticles carry a nonzero baryon number. This can be seen by the generator of 
baryon number conservation, 

S = 5--^=T 8 = diag c (0,0,l) > (113) 

where Tg is the eighth generator of the [SU(3) C ] group. The electromagnetic generator of the 
unbroken [£/(l) cm ] gauge symmetry in the 2SC phase is 

Q = Q--L Ta , Ql = Ql = l, Q b u = l, CT d = Ql = - 1 -, Q b d = 0. (1.14) 

where Q — diagy (§,—5) is the electromagnetic generator of the [t/(l) cm ] gauge symmetry in 
vacuum. The gauge boson of [f/(l) em ] is the medium photon. Therefore, there exists no elec- 
tromagnetic Meissner effect in the 2SC phase and that is the reason why a magnetic field would 
not be expelled from the color-superconducting region. The 2SC phase is a so-called Q-conductor 
because its electrical conductivity is large. The Q charge of the blue up quasiparticle is responsible 
for this behavior, see Eq. (|1.14p . 

At weak coupling, the difference of the color-superconducting pressure to the pressure of 
normal-conducting quark matter for massless non-interacting quarks at zero temperature (| 1 .3(1 
amounts to 

u 2 A 2 

per each gapped quasiparticle [17]. In the 2SC phase without strange quarks, there are six quarks 
from which four of them are gapped so that the pressure for color-superconducting quark matter 
in the 2SC phase at zero temperature approximately reads, 

u 4 a 2 A 2 

Pmc^^ + ^-B. (1.16) 

1.2.2 The CFL phase 

If the strange quark chemical potential exceeds the mass of strange quarks, then there are also 
strange quarks in quark matter at zero temperature, see Eq. flTTT]) . Therefore, it is possible that 
also the strange quarks participate in pairing. Since the spin part of the quark Cooper pair 
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wavefunction is antisymmetric for J = and quarks pair in the antisymmetric color-antitriplet 
channel, also the flavor part has to be antisymmetric in order to fulfill the requirement of overall 
antisymmetry This leads to the following color-flavor gap matrix in the three-flavor case: 

A{? = e ijk e^ h A h k , (1.17) 

where the order parameter is given by 

A£ = #A. (1.18) 
Together with the Dirac part, the gap matrix reads 

$(; = 75 A{/. (1.19) 

There is one difference in comparison to the 2SC phase: the antisymmetric tensor in Eq. (11.17P 
now possesses three instead of two flavor indices because three instead of two quarks participate in 
pairing, Nf = 3. The condensate breaks [SU(3) C ] x SU (3) r+ £ to the vectorial subgroup SU (3) c+r+ £ 
and is still invariant under vector transformations in color and flavor space. This means that a 
transformation in color requires a simultaneous transformation in flavor to preserve the invariance 
of the condensate. Therefore, the discoverers [18] of this three-flavor color-superconducting quark 
state termed it the CFL phase which is an abbreviation for color-flavor-locked phase. The CFL 
phase is the true ground state of quark matter because all quarks are paired which leads to the 
highest pressure of all color-superconducting phases. 

In contrast to the 2SC phase, the CFL phase also has gaps in the repulsive sextet channel. 
The color-flavor gap matrix (|1.17[) can be extended by the (small) symmetric sextet gaps, 

A{/ = A (3 ,3) - SfS?) + A (6 , 6) (5{6? + 5° 8$) , (1.20) 

where A (3 3) is the antitriplet and A(g,6) is the sextet gap. This can be rewritten as 

A{/ = A[6{ 6 9 j + A' 2 5f5 f 3 , (1.21) 

where A^ = A^g gj + A( 6 6 ) and A' 2 = — A^g) + A( 66 ). By introducing the color-flavor projec- 
tors [19], 

M! = \ 6 ( 69 3 > M! = \^ fa - \^ , M! = l s « sh + - , (1-22) 

which fulfill the properties of completeness, ~Yi, n Vn — 1, and orthogonality, PiPj = SijPj, the 
color-flavor gap matrix assumes the form, 

A{? = An [Pn]{- = \ (Ax + A 2 ) 8{Vl - A 2 5?Sf , (1.23) 

n=l 6 

where A3 = — A 2 - The singlet gap Ai and the octet gap A2 appear in the spectra of quasiparticles 
and quasiantiparticles, 

el (k) = ^(k-etf + \A n \ 2 , (1.24) 

where n — 1,2, and all quarks are treated as massless for simplicity. By neglecting the small 
repulsive sextet gap, one finds that Ai = 2A 2 = 2A. 

In the CFL phase, there are no gapless quasiparticles. At small temperatures (T <C A), the 
contributions of the quark quasiparticles to all transport and many thermodynamic quantities are 
suppressed by the small exponential factor exp(— A/T) [9]. The influence of the gluons on the CFL 
phase is negligible because all of them are massive because of the color Meissner effect. In contrast 
to the 2SC phase, the CFL phase is superfluid because the U(1)b baryon number symmetry is 
broken, but it has an unbroken [J7(l) om ] gauge symmetry and therefore it is, like the 2SC phase, 
not an electromagnetic superconductor. This is the reason why the CFL phase does not expell a 
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magnetic field from its color-superconducting interior. The electromagnetic generator of the CFL 
phase reads, 

Q = Q-T 3 -^=T S , Q r u = Q 9 d = Q 9 s = Q h d = Q h s = Q, Ql = Q b u = l, Q r d = Q r s = -l. 

(1-25) 

The CFL phase is a Q-insulator because all quarks are gapped and there is no remaining electric 
charge as in the 2SC phase. Therefore, the CFL phase is electrically charge neutral. At zero 
temperature, there are no electrons present [20]. At small temperatures, the electrical conductivity 
of the CFL phase is dominated by thermally excited electrons and positrons [21] and becomes 
transparent to light [21,22]. 

As in the case of chiral perturbation theory in vacuum QCD [23], one could write down an 
effective low-energy theory in the CFL phase. From the symmetry breaking pattern, it is known 
that there are nine Nambu-Goldstone bosons and one pseudo-Nambu-Goldstone boson in the low- 
energy spectrum of the theory [24]. Eight of the Nambu-Goldstone bosons are similar to those 
in vacuum QCD: three pions (n and ir^), four kaons (K°, K°, and K^) and the eta- meson (77). 
The additional Nambu-Goldstone boson (</>) comes from breaking the U{\)b baryon symmetry. 
In absence of gapless quark quasiparticles, this Nambu-Goldstone boson turns out to play an 
important role in many transport properties of cold CFL matter [21,25]. Finally, the pseudo- 
Nambu-Goldstone boson (77') results from breaking of the approximate axial U(1)a symmetry. A 
possible phase transition to the CFL phase with a meson (e.g., kaon or eta) condensate could 
happen if M a > M^ 3 A 2 / 3 [26-30], where M u is the up and M s the strange quark mass. 

In the CFL phase, there are nine quarks, and all of them are gapped: one quasiparticle with gap 
Ai and eight quasiparticles with gap A2. Therefore, the pressure of the CFL phase for massless 
quarks at zero temperature approximately is given by 

3^ 4 fi 2 Al 2/i 2 A| „ 3^ 4 3/i 2 A 2 „ , „„, 

where Eqs. (fOj) . (fl~T3| . and A x = 2A 2 = 2 A are used. 



1.2.3 Spin-one color superconductivity 

Since quarks pair in the antisymmetric color-antitriplet channel and the spin part of the quark 
Cooper pair wavefunction is antisymmetric for J = 0, condensation with only one flavor is for- 
bidden by the Pauli principle, but it is possible for the J = 1 channel, where the spin part of 
the wavefunction is symmetric. Thus, the Cooper pair wavefunction is, as required, overall an- 
tisymmetric. Spin-one color superconductivity [12, 31-37] is much weaker than spin-zero color 
superconductivity. The gap in spin-one color-superconducting systems is of the order of 100 keV. 
Such a small gap will not have big influences on the transport and many thermodynamic proper- 
ties of the quark matter [9] . Spin-one color superconductivity is less favored than spin-zero color 
superconductivity since the latter has a higher pressure because of the larger gap. This is why one 
does not expect that spin-one color-superconducting quark phases dominate in the phase diagram 
of neutral quark matter. But they could be favored if it is not possible to form a spin-zero color- 
superconducting state because of a too large mismatch between the Fermi surfaces of different 
quark flavors [38]. 

The general structure of the gap matrix for spin-one color-superconducting systems reads [36, 

37], 

3 



A ab = iA eabcc c 

c,i—l 



k cos 8 + 7j_ sin < 



(1.27) 



where 7^ = 7* — k l (7 • k), k = k/fc, k is the momentum vector, and k its absolute value. Spin-one 
color-superconducting phases are called longitudinal if 8 = and transverse if 8 = ir/2. Many 
different spin-one color-superconducting phases can be constructed by choosing various specific 
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3 x 3-matrices C. The most important of them are the A-phase, the color-spin-lockcd (CSL) 
phase, the polar phase, and the planar phase [36,37], 



/0 o\ 1 (1 0' 

c (A-phasc) = _L o o C (CSL) = _L 1 



\i * / ^ 



,0 1 



/o o o\ 1 /l 0' 

C (polar) _ q q Q J C (planar) _ J_ Q j q 

VO 1 / ^2 1 o 



(1.28) 



which are characterized by different symmetries of their ground state. 

The original group [SU(3) C ] x 50(3) j x [U(l) em ] breaks down as follows [36,37]: 

A-phase: SU(2) C x SO(2)j x U(T) em , 

CSL: 50(3)j , 

Polar: SU(2) C x SO(2)j x C/(l) cm , 
Planar: SO(2)j x f7(l) cm . 

In spin-one color superconductors, there can exist an electromagnetic Meissner effect in con- 
trast to spin-zero color superconductors. If so, magnetic fields will be expelled from the color- 
superconducting region. This is, for example, the case in the CSL phase. The most energetically 
preferred spin-one color-superconducting quark phase is the transverse CSL phase because it has 
the highest pressure [39]. 



1.3 Stellar evolution 

Stars begin their life as objects which are formed out of contracted interstellar gas clouds in 
galaxies. They are nuclear burning factories: light nuclei, such as hydrogen, will be burned to 
heavier nuclei by fusion reactions. After all fusion reactions are completed, stars end their life as 
compact stars: white dwarfs, neutron stars, or black holes. 

1.3.1 The formation of stars 

Interstellar clouds which mainly consist of hydrogen contract if their gravitation exceeds the pres- 
sure from inside caused by turbulence and temperature. In order to obtain such a large gravitation, 
interstellar clouds have to possess big masses. The so-called Jeans criterion has to be fulfilled so 
that an interstellar cloud is able to contract [40], 

M > M CI (T, Pc ) ~ 28TI/J-3 , (1.29) 

where T is the temperature in Kelvin and p c the central density of the interstellar cloud in cm -3 . 
The critical mass M cr is obtained in units of the solar mass. In Table 11.21 some values for the 
critical mass are shown. 





1 cm 3 


100 cur 3 


10 4 cm" 3 


10 K 
100 K 


880 Mq 
28000 M Q 


88 M 
2800 M 


8.8 Mq 

280 M 



Table 1.2: Values for the critical mass of interstellar clouds [40]. 



By the contraction of the interstellar cloud, first stars of spectral type O are created in its 
center. Their ultraviolet radiation ionizes the hydrogen gas around them so that it becomes hot 
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and is shining. This so-called H II region has a temperature T ~ 10000 K and expands into the 
cool outer regions of the interstellar cloud which consist of cold hydrogen gas, so-called H I regions 
which have a temperature T ~ 100 K. Wavy bays and globules are created by this expansion. The 
globules have diameters of up to one parsec (pc) and masses of up to 70 M . Within 500000 years, 
they fragment and collapse to protostars which emit only infrared radiation because there are dense 
dust clouds around them which fall on them within several million years. In the meantime, the 
contraction of the star continues until the pressure from inside becomes as large as the gravitation. 
Then, the star is on the main sequence in the Hertzsprung-Russell diagram (HRD), see Fig. 11.31 
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Figure 1.3: The Hertzsprung-Russell diagram (HRD) [41]. 



The massive and bright stars are in the upper and the darker stars with lower mass are on the 
lower part of the HRD. Hot blue stars are in the left and cooler red stars are in the right part 
of the HRD. During the contraction of gas balls to stars, their rotation becomes faster and faster 
because of angular momentum conservation. This could lead to a splitting of the gas balls because 
of large centrifugal forces so that gravitationally bound narrow double, multi-star, or solar systems 
will be formed. Wider gravitationally bound double or multi-star systems are formed if the stars 
of the interstellar cloud come close to each other by their movement and by gravitational forces. 
In this way, many stars and star systems are created inside interstellar clouds so that finally there 
exists an open cluster in their center from which an O association, a star cluster with hundreds of 
stars of spectral type O until B2, expands into the outer regions. 



1.3.2 Main sequence stars 

During the contraction, the central temperature T c of the stars becomes higher and higher. At 
T c ~ 5 ■ 10 6 K, hydrogen burning is initiated in the cores of the stars. For producing such 
large temperatures in the cores, the stars need at least a mass M ~ 0.1 M Q . Otherwise brown 
dwarfs will be created in which hydrogen burning never occurs. Hydrogen-burning stars are on 
the main sequence of the HRD. The contraction is finished because the pressure is as large as 
the gravitational attraction. There are two possibilities to burn hydrogen into helium by fusion 
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reactions [41,42]: these are the proton-proton cycles (pp cycles), see Fig. I1.4[ and the carbon- 
nitrogen-oxygen cycles (CNO cycles), see Fig. 11.51 
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Figure 1.4: The three proton-proton cycles (pp cycles). Percent numbers specify how large the probabil- 
ity of a reaction is. The star at 8 Be means that it is not stable and decays into two a-particles. 
Cycle I was discovered in Ref. [43]. Positrons which are formed by fusion reactions will be 
annihilated. 
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Figure 1.5: The carbon-nitrogen-oxygen cycles (CNO cycles) [44]. Cycle I is named after their dis- 
coverers: Bethe-Weizsacker cycle [45]. Cycle II is called CNO bi-cycle, and cycle III is 
termed CNO tri-cycle. The widths of the arrows illustrate the significance of the reactions 
in determining the nuclear fusion rates in the stellar CNO cycles. Certain hot CNO pro- 
cesses are indicated by dashed arrows. Positrons which are formed by fusion reactions will 
be annihilated. 



The main difference between these two possibilities is that in the pp cycles protons fuse directly 
with each other while in the CNO cycles, they are always burnt into helium by using carbon, 
nitrogen, oxygen, etc. as catalysts. In cooler yellow and red stars, where T c ~ 5 • 10 6 K, only 
the pp cycles are present. The CNO cycles are initiated not until T c ~ 9.5 • 10 6 K, but they are 
dominant in hot blue and white stars where T c > 1.5 • 10 7 K. The energy production of the CNO 
cycles is then much larger as the energy production of the pp cycles. The larger the mass of a 
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(main sequence) star the higher is its central temperature. Stars with big masses use the CNO 
cycles to burn hydrogen into helium. But only twelve percent of the total hydrogen is in the core 
of a main sequence star which can be burnt into helium. Because of the much higher central 
temperatures, blue and white stars with large masses burn the hydrogen in their cores much faster 
into helium as cold yellow and red stars with low masses, see Table [TBI 



Type 


M [Mq] 


Million years 


05 


39 


0.5 


BO 


20 


4.6 


B5 


6.7 


46 


AO 


3.5 


319 


A5 


2.2 


1160 


FO 


1.7 


2700 


F3 


1.26 


3800 


F6 


1.13 


6000 



Table 1.3: Burning period of hydrogen in the cores of main sequence stars [40]. 

1.3.3 Red giants and red super giants 

When the hydrogen in the core is used up, gravitation becomes dominant so that the star contracts 
because there are no fusion reactions any more which can maintain the pressure from inside. This 
leads to higher central temperatures, T c > 100 million Kelvin. Such high temperatures are needed 
to initiate helium burning in the core of a star, 

a (a, 7) 8 Be* (a, 7) 12 C (a, 7) 16 . (1.30) 

To produce such high temperatures in the core, a star needs at least a mass M ~ 0.4 M Q . 
Otherwise the star is only able to burn hydrogen and then ends as a white dwarf. Helium burning 
happens in a process of two stages, a so-called Salpeter process [46] : two helium nuclei form an 
intermediate stage, a 8 Be nucleus which is unstable and is therefore marked by a star in the fusion 
reaction (|1.30p . It reacts with a third helium nucleus to produce 12 C. The Salpeter process of 
helium burning is called triple-a process, and it is much more likely as a direct reaction of three 
helium nuclei to form carbon. A further reaction of 12 C with a helium nucleus produces 16 0. 

Above the core of helium burning, there remains a shell of hydrogen burning which expands 
until it reaches the star surface. The star blows up because of the high pressure caused by helium 
burning. Its diameter and luminosity are much larger as in the period of a main sequence star but 
the surface temperature is colder so that the star appears red. It has become a red giant. When 
hydrogen burning ends and helium burning is initiated, the stars go from the main sequence into 
the giant region in the HRD. Red giants are not as stable as main sequence stars because the 
pressure from inside is in imbalance with the gravity of the star. This leads to oscillations of red 
giants: the star contracts, the radius becomes smaller and the star whitens a little. The interior 
pressure increases, the star expands, the radius becomes larger and the star becomes a little more 
red. This procedure repeats again and again. When the helium in the core is used up, the red 
giant contracts. If its mass exceeds eight solar masses, the central temperature is hot enough to 
burn carbon in the core of the star, 

12 C + 12 C — » 20 Nc + a, 12 C + 12 C — > 23 Na + p, 12 C + 12 C — > 23 Mg + n , (1.31) 

otherwise the red giant ends up as a white dwarf. By carbon burning, the red giant blows up to 
a red super giant because the interior pressure becomes very high. If all carbon is used up, the 
red super giant contracts. If the mass of the red super giant exceeds ten solar masses, its central 
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temperature becomes high enough so that neon, oxygen, and silicon will be burnt in the core, 

20 Ne + 7 — > 16 + a, 20 Nc + a — > 24 Mg + 7 , 24 Mg + a — > 28 Si + 7 , (1.32) 
16 + 16 — > 31 P + p, 16 + 16 — > 28 Si + a, (1.33) 
7 + 28 Si — » (p, a, n) , 28 Si (a, 7) 32 S (a, 7) 36 Ar . . . 52 Fe (a, 7) 56 Ni , (1.34) 

otherwise the fusion reactions end with the fusion of carbon. After all carbon is used up in the 
core of the red super giant and its mass is lower than ten solar masses so that the red super giant 
is not able to proceed with further fusion reactions, it finally explodes by a supernova of type II 
and forms a neutron star. But if the mass of the red super giant exceeds ten solar masses, fusion 
reactions can happen until nuclei with mass number A — 56 are produced. These are iron and 
nickel. For larger A, fusion reactions would need energy to produce heavier nuclei. So, there is 
no other way for a massive red super giant with an iron core as to collapse because of the missing 
interior pressure. It will cause a supernova of type II and form a neutron star. If the mass is even 
larger than approximately twenty solar masses, the star will collapse into a black hole. 



1.3.4 Compact stars 

If the mass of a star is smaller than approximately 0.1 M©, hydrogen burning is never initiated 
and the star will end up as a brown dwarf. Otherwise, if the mass is larger, the star will end 
up as a compact star: a white dwarf, a neutron star, or a black hole. If the mass is larger than 
approximately 0.1 Mq but does not exceed eight solar masses, the star will end up as a white 
dwarf. The outer shells of red giants are blown away and form a planetary nebula. There are no 
fusion reactions any more so that the core of the red giant collapses to a white dwarf in the center 
of the planetary nebula. Only the pressure of degenerate electrons is able to stop the collapse 
which is caused by gravity. White dwarfs have densities of about one million grams per cubic 
centimeter, radii of a few thousand kilometers, and cool approximately within ten billion years to 
black dwarfs and become invisible. 

If a star A is a star with a small mass in a narrow binary system, matter from its massive 
companion star B is accreted by star A if star B exceeds its Roche volume when it becomes a red 
giant. This phenomenon has a big influence on the evolution of stars in narrow binary or multi-star 
systems: star A becomes much more massive than star B. If star A is a white dwarf and star B is a 
red giant, the white dwarf accretes matter which can be seen as nova bursts because the accreted 
matter causes nuclear reactions. If the mass of the white dwarf exceeds the Chandrasekhar mass 
which is approximately 1.44 Mq, the white dwarf explodes in a supernova of type I because its 
gravity becomes so large that the degenerate electron pressure is not able to stabilize the star any 
more. Expanding gas shreds are the only remnants of the white dwarf. 

If the mass of the red super giant is larger than eight solar masses and there are no further 
fusion processes in the star, the pressure of degenerate electrons is not able to stop the collapse 
any more. Such red super giants collapse in supernova explosions of type II to neutron stars. More 
details about neutron stars will be presented in Sec. 11.41 

A black hole will be formed if the progenitor star, a red super giant, possesses a mass of about 
20-25 solar masses. No internal force is able to stop the collapse of the red super giant because of 
the huge gravitation. A black hole is defined as a region of space-time which cannot communicate 
with the external universe which means that there is no chance for a particle or even for light to 
escape from a black hole if it has reached the region beyond its event horizon. The event horizon is 
the boundary of a black hole at which gravity is so strong that nothing has a chance to escape. The 
radius of the event horizon for a spherical mass is called the Schwarzschild radius, r$ = 2GM/c 2 . 



1.4 Neutron stars 

At the end of its life, a red super giant with a mass larger than eight solar masses consists of many 
shells of different nuclei which are created by fusion processes. They are arranged like onion shells: 
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the surface shell of the red super giant consists of the lightest nucleus, hydrogen. Towards the 
center, the nuclei which were created by fusion processes get heavier. In the center of red super 
giants with masses larger than ten solar masses, there are iron nuclei. If all fusion processes are 
finished, the red super giant will collapse. Since gravity of such massive stars is so dominant, not 
even the pressure of degenerate electrons is able to stop the collapse. The iron nuclei in the center 
of the star break up because of the high pressure and temperature. This process is called photo 
dissociation, 

7 + 56 Fe — >13a + 4n, (1.35) 

costs energy, and the thermal pressure from inside reduces so that the star collapses in approxi- 
mately 0.1 seconds. Electrons are captured by protons so that matter in the core mostly consists 
of neutrons, 

p + c~ — ► n + v e . (1.36) 

The pressure of these degenerate neutrons stops the collapse. The process (|1.36[) is called inverse 
(3 decay in which a huge amount of neutrinos is produced. These neutrinos arc trapped for a while 
in the collapsed hot star core which is called a protoneutron star. The outer shells fall down on 
this protoneutron star, bounce back and thereby produce an outgoing shock wave which is seen 
as a supernova explosion of type II in which nuclei with mass numbers A > 56 are formed. 
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Figure 1.6: The main stages of the evolution of neutron stars [49] . Shading indicates approximate relative 
temperatures. 



The evolution of protoneutron stars goes through several stages [47-49] as is shown in Fig. 11.61 
During the supernova explosion, there goes a shock through the outer mantle of the protoneutron 
star. The outer mantle is of low density and high entropy, accretes matter and loses energy by 
f3 decays and neutrino emission. The core has a mass of about 0.7 Mq in which neutrinos are 
trapped. The lepton to baryon fraction is Yl — 0.4. The whole protoneutron star in this first 
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stage of evolution has an approximate radius of about 200 kilometers. After approximately five 
seconds, accretion becomes less important. The mantle collapses because of the loss of lepton 
pressure caused by deleptonization. In this stage of evolution, the lepton to baryon fraction is 
in the range < Yj, < 0.4, and the radius amounts to approximately 20 kilometers. If a lot 
of matter accretes onto the protoneutron star so that it exceeds its maximum mass then the 
protoneutron star collapses to a black hole. At approximately 15 seconds after the supernova 
explosion, the protoneutron star is dominated by neutrino diffusion causing deleptonization and 
heating of the protoneutron star. In this stage of evolution, the lepton to baryon fraction is 
Yl = 0, the radius of the protoneutron star is approximately 15 kilometers, and its temperature 
is heating up to 30 MeV < T < 60 MeV. There is also the possibility of forming a black hole 
by deleptonization. Fifty seconds after the supernova explosion, the protoneutron star becomes 
transparent to neutrinos so that the inner part of the star cools down. But the crust remains 
warm because of its lower neutrino emissivity, T ~ 3 • 10 6 K. After 50-100 years, also the crust 
cools down by neutrino emission and the star becomes isothermal. In later stages, the star cools 
down by direct URCA processes, 

n — ► p + e~ + D c , p + c~ — ► n + u e , (1-37) 

or modified URCA processes, 

n + n — ► n + p + c~ + v c , n + p + c~ — > n + n + v c , (1.38) 

neutrino and photo emission. A cold neutron star has been formed. 

1.4.1 Pulsars 

Before neutron stars were discovered, theorists speculated about the existence of neutron stars. 
In 1932, Landau called them weird stars. In 1934, Baade and Zwicky realized that there is a 
connection between supernovae of type II and neutron stars. The first neutron star calculations 
were done by Tolman, Oppenheimer, and Volkoff in 1939 who created mass-radius diagrams of 
neutron stars [50]. Further work has been done by Wheeler et al. (1960-1966) and Pacini (1967). 

In the summer of 1967 in Cambridge/England, Jocelyn Bell, a student of Anthony Hewish who 
got the Nobel prize in 1974, detected a neutron star as a pulsar for the first time. Neutron stars 
can be observed as pulsars which are pulsating sources of radiation because they can be identified 
by their very precise radio pulses. Pulsars possess a strong magnetic field of about 10 12 G in which 
highly energetic electrons gyrate and thereby produce synchroton radiation which is emitted at 
the magnetic poles of the pulsar, see Fig. 11.71 Usually, the rotation axis of a pulsar is inclined to 
the axis of the magnetic field. Thereby, the cone of the synchroton radiation can be detected only 
once in a rotation period of the pulsar. This is the pulsation phenomenon of pulsars which is also 
called the lighthouse effect. 

From observations with radio telescopes one knows that pulsars have periods in the range of 1.6 
milliseconds to several seconds. They rotate so fast because of angular-momentum conservation. 
In Sec. 11.3.11 it was mentioned that stars (gas balls) rotate. This rotation is kept during the life 
of the stars. When red super giants collapse to pulsars, the rotation velocity increases very much. 
The pulsar periods are very stable and therefore increase not much in time. The Crab pulsar for 
example has a rotation period of 33 milliseconds and this changes only about 0.036% per year. 
From the present pulsar period P and its time derivative one is able to determine the characteristic 
age t = P/P of a pulsar. Because of the rapid rotation, a pulsar has an oblate shape and therefore 
is not exactly spherically symmetric. Observations like gravitational redshift measurements and 
mass determinations in binary systems as well as theoretical calculations show that the masses 
of neutron stars or pulsars, respectively, approximately amount to 1.5 Mq, and that they have 
radii of about ten kilometers. Therefore such compressed, massive objects have an unbelievable 
mass density which is of about 10 14 g/cm 3 . A further observation in pulsars are so-called glitches, 
sudden small jumps in the rotation period of pulsars. They are most probably caused by vortices 
and rearrangements in the crust of the pulsar. Because of that, it comes to a decrease of the 
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angular momentum in the superfluid and an increase of the angular momentum in the crust [51]. 
Pulsars show proper motions which originate in so-called pulsar kicks. A possible explanation 
is that they are created if neutrinos emit asymmetrically during the supernova explosion which 
leads to a propulsion of the pulsar. Another observation in pulsars of binary systems are x-ray 
emissions and bursts: mass is transferred from companion stars onto the accretion disc of the 
pulsars. Because of the strong magnetic fields of the pulsars, matter from the accretion disc is 
diverted to the poles of the magnetic field. At this places, nuclear fusion processes are initiated, 
which emit x-rays. 




Figure 1.7: A schematic representation of a pulsar [7]. 



1.4.2 Structure of neutron stars 

In Sec. 11.41 it was shown that neutron stars consist of neutrons because of inverse /3 decays (|1.36[) . 
This does not mean that neutron stars only consist of neutrons. For different mass densities, there 
are different phases of matter inside neutron stars, see Fig. 11.81 

Neutron stars consist of an atmosphere of electrons, nuclei, and atoms. Only a fraction of the 
electrons are bound to nuclei. The equation of state was calculated by Feynman, Metropolis, and 
Teller [52]. They found out that the nuclei in this regime, where the mass density p < 10 4 g/cm 3 , 
are mainly 56 Fe. 

In neutron stars with temperatures above typically 100 eV, between the atmosphere and the 
solid crust a layer is present, where nuclei and electrons are in a liquid phase called the ocean [53]. 

One assumes complete ionization of the atoms, when the spacing between the nuclei becomes 
small compared to the Thomas-Fermi radius ttf — clqZ _1 / 3 of an isolated neutral atom. In 
this equation, ao is the Bohr radius and Z the charge number. The mass density approximately 
amounts to p ~ Am u nN, where A is the mass number, m u the atomic mass unit, and njv the 
number density of nuclei, which depends on the radius of a spherical nucleus whose volume is 
the average volume per nucleus [54], §7rr 3 = l/n N . By combining the last three equations, one 
finds that the outer crust of cold neutron stars begins when p ~ 10 4 g/cm 3 3> 3AZ g/cm 3 . 
This shell consists of nuclei and free electrons. The equation of state was originally calculated 
by Baym, Pethick, and Sutherland (BPS) [55] and improved in Refs. [56, 57] using up to date 
nuclear data. The BPS model is valid for zero temperature which is a good approximation for 
the crust of nonaccreting cold neutron stars. The outer crust of nonaccreting cold neutron stars 
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contains nuclei and free electrons. The latter become relativistic above p ~ 10 7 g/cm 3 . The nuclei 
are arranged in a body-centered cubic (bcc) lattice. The contribution of the lattice has a small 
effect on the equation of state but it changes the equilibrium nucleus to a larger mass number and 
lowers the total energy of the system because it will minimize the Coulomb interaction energy of 
the nuclei. The latter are stabilized against (3 decay by the filled electron sea. At p ~ 10 4 g/cm 3 , 
56 Fe is the true ground state. With increasing mass density, it is not the true ground state any 
more because the nuclei capture electrons, emit neutrinos and become neutron richer. When the 
mass density p ~ 4.3 • 10 11 g/cm 3 , the so-called neutron drip line is reached. Neutrons begin to 
drip out of the nuclei and become free. This happens because the equilibrium nuclei become more 
and more neutron-rich, and finally no more neutrons can be bound to nuclei. At the neutron drip 
point, the inner crust of neutron stars begins. At p ~ 2 • 10 14 g/cm 3 nuclei do not exist anymore, 
signalling the end of the neutron star crust. 




R ~ 10 km 



Figure 1.8: Cross-sections of neutron stars (schematic) [58]. Each sector shows another possible type of 
neutron stars. 



The equation of state of the inner crust was calculated by Baym, Bethe, and Pethick [59] and 
another equation of state for this regime was derived by Negele and Vautherin [60] . Also a relativis- 
tic mean field model has been used to describe the density regime of the neutron star crust within 
the Thomas- Fermi approximation (see Ref. [61] and references therein). For higher densities, the 
nuclei disintegrate and their constituents, the protons and neutrons, become superfluid. Muons 
also begin to appear in these shells. The equation of state of this regime can be calculated by 
using non-relativistic many-body theories [62] or relativistic nuclear field theories [63-66]. What 
kind of matter exists in the cores of neutron stars depends on their central densities. Neutron 
stars with lower central densities consist of protons, neutrons, electrons, and muons while others 
with larger central densities can consist of hyperons, a pion or kaon condensate. In neutron stars 
with huge central densities, n ss 10no, even the protons, neutrons, and hyperons can disintegrate 
into their constituents: quarks which are deconfined. Such neutron stars contain a quark core [58] 
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whose true ground state is strange quark matter [67] so that not only the light up and down quarks 
but also the strange quarks occur if their mass is low compared to the strange quark chemical 
potential. The remaining three quark flavors are too heavy to participate in the quark matter 
of neutron stars. Because of the dominant attractive interaction in the antitriplet channel, the 
true ground state of strange quark matter in the cores of neutron stars is color-superconducting 
strange quark matter [3,9,68-73]. 




2 rTkml 8 11 



Figure 1.9: A schematic representation of the quark-hadron phase transition in neutron stars with a 
mass of about 1.4 M Q [74]. 



In nature, phase transitions can take place either through sharp boundaries between pure 
phases which are located next to each other or through mixed-phase regions. Thus, it does not 
necessarily mean that the phase transitions in neutron stars have to take place as sudden as 
represented in Fig. 11.81 In Fig. 11.91 one can see that the quark-hadron phase transition is not 
separated into two distinct shells but that there is a smooth transition from the hadronic phase 
into the quark phase. This transition goes in several stages, and by the effects of Coulomb forces 
and surface tensions, some interesting structures are formed. Quarks at first appear in a structure 
of drops in a shell under the pure hadronic phase in which the amount of hadrons exceeds the 
amount of quarks. With increasing density, quarks form a structure of rods, and finally a structure 
of slabs and are surrounded by hadronic matter. At higher density, the amount of quarks exceeds 
the amount of hadrons which form a structure of slabs. With increasing density, hadrons form a 
structure of rods, and finally a structure of drops that are surrounded by quark matter. Ultimately, 
the cores of neutron stars with huge central densities consist of pure quark matter. 

1.4.3 Properties of neutron star matter 

The matter in neutron stars is in its ground state. It is in nuclear equilibrium which means that 
the energy cannot be lowered by strong, weak, or electromagnetic interactions. When matter is 
in equilibrium concerning weak interactions, one calls it /^-equilibrated matter or one says that 
matter is in f3 equilibrium. It means that the reactions for (inverse) (3 decay (|1.36[) or the reactions 
for weak interactions, respectively, are in equilibrium for all lepton families, 

p + /r< — >n + v . (1.39) 

In this equilibrated reaction, (3 = (e,^t, r). The leptons are listed in Table fOl They are spin-i 
particles, and therefore fermions. 
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Lepton 


TV jT IT t rl 

Mass [MeVJ 


Q [e] 


9 


electron (e) 


0.51099906 ±0.00000015 


-1 


2 


electron neutrino (u e ) 


< 7.3 ■ 10" 6 





1 


muon (/i) 


105. 658389 ± 0.000034 


-1 


2 


muon neutrino (v^) 


< 0.27 





1 


t lepton (r) 


1776.3 ±2.4 


-1 


2 


r neutrino (y T ) 


< 31 





1 



Table 1.4: The masses, electric charges, and degeneracy factors g of the leptons [16]. The abbreviation 
for the respective lepton is shown in brackets. There also exists the corresponding antiparticle 
for each lepton with charge — Q. They are not shown here. 

Since the quark flavor content of protons and neutrons is 

p=(u,u,d), n=(u,d,d), (1.40) 

one can express the reaction (|1.39|) in terms of quark flavors, 

u + /r< — >d + v p . (1.41) 

By weak interactions, also the transformation from up into strange quarks is possible, 

u + /T< — >s + vp. (1.42) 

Neutrinos carry only lepton number. That is why the chemical potential of each neutrino family 
is equal to the chemical potential of each lepton family, 

Vup = fJ-L/3 ■ (1-43) 

From the equilibrated reactions (|1.41[) and (|1.42p . one can directly write down the equation for 
the corresponding chemical potentials, 

Mu + M/3 = ^d,s + fJ-Lff , (I- 44 ) 

where [id,s means that fi s — fi^, a fact that directly comes out of the equilibrated reactions (| 1 .4 1 1) 
and (|1.42p . The last equation can be rewritten as 

Mu = Md,s + . (!-45) 

where 

MQ = MLp - M/3 (1-46) 

is the chemical potential of electric charge. This means that electrons, muons, and r leptons carry 
both, lepton number and electric charge, 

M/3 = VL (i - I^q ■ (1-47) 

Because of the quark content of the neutron (|1.40p , one can define the baryon chemical potential 

fi B = 3/i = fi n = /i„ + 2fi d , (1-48) 

where fi is the quark chemical potential and fi n the chemical potential of neutrons. This relation 
together with the fact that /i s = Hd is used to solve the equation (|1.45[) for each quark flavor. One 
obtains 

t*ff' = fAff + VQQff' , (!-49) 

where 

e = diag(|4-i) (1.50) 
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is the matrix of electric charge in flavor space. 

With Eq. (|1.49p . one automatically satisfies (3 equilibrium in normal quark matter where the 
color symmetry [SU (3) c ] is not broken. This is not the case in color-superconducting quark matter. 
There, one has to know the quark chemical potential of each quark color and flavor. In order to 
satisfy equilibrium for each quark color and flavor, one starts with the equation for (3 equilibrium 
for each quark flavor (|1 .49|) and adds to it the terms for each color which consist of color chemical 
potentials and the generators of the [SU(3) C ] group [75], 

8 

»%, = {jtSff, + MqQ//') S"' + 5 ff > ]T /x c Tf . (1.51) 

c=l 

In this equation, the color indices i and i' are superscripted while the flavor indices / and /' are 
subscripted. But not all of the eight color chemical potentials are nonzero. This can be proven by 
calculating the tadpoles and was done for the 2SC phase in Ref. [76] where only a nonzero color 
chemical potential fig is present. I extended this calculation for the color-superconducting phases 
which I had to investigate for this thesis and got the result that there are nonzero color chemical 
potentials ^ and fig. Therefore, the chemical potential for each color and flavor can be simplified 
to 

4 f , = {fiSff, + n Q Q fr ) 5 U ' + (^ 3 Tf + M8 7f ) S Jf , (1.52) 

for the purposes of this thesis. Later on, I omit the double color and flavor indices of the quark 
chemical potential matrix (| 1.520 and denote fj^f, by fjij because the quark chemical potential 
matrix (|1.52[) is diagonal in color-flavor space. Later on, also the double flavor indices of the 
matrix of electric charge (|1.50p will be omitted because it is diagonal in flavor space, e.g. Q/f is 
denoted by Qf. Also, the double color indices of the generators of the [SU(3) C ] group, Tf and 
Tg* , can be omitted because they are diagonal in color space, e.g. Tf can be denoted by T£, and 
If' can be denoted by T 8 \ 

Stars are bound by gravity and have to be electrically charge neutral, otherwise they would 
be unstable and explode because of repulsive Coulomb forces. The number density of electrically 
charged particles is given by 

riQ = (ip f Qijj) - n e - - n T , (1.53) 

where ip is the quark spinor, and n e , n^, and n T are the number densities of electrons, muons, and 
t leptons, respectively. The electrical charge neutrality condition, 

n Q = p- = , (1.54) 

demands that the number density of electrically charged particles uq , which can be calculated by 
taking the derivative of the pressure p with respect to the chemical potential of electric charge, is 
equal to zero. 

Stars without color-superconducting quark matter are automatically color neutral because the 
[SU (3) c ] color symmetry is not broken. This is not the case in color-superconducting quark matter. 
If stars consist of color-superconducting quark matter, then they have to be color neutral because 
on the one hand, one is not able to observe color charges in nature, on the other hand, stars will 
not be stable if they are not color neutral. In the following, I show that the color number densities 
ri3 and n% have to be equal to zero in order to fulfill color neutrality. 

The spinor of quark colors is defined by 

ijj = (■tp r7 tfj g ,ip b ) T . (1.55) 
Herewith, the number densities of quarks read, 

n r ee (V4V>r) , n g ee (ipliJg) , n b = , (1.56) 
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so that 



n = (ip'ip) — n r + n g + n b 
n 3 = (^T 3 ip) = § (n r - n g ) 
1 

2V3 



n 8 = (^T 8 ip) 



2n b ) 



(1.57a) 
(1.57b) 
(1.57c) 



I only need to consider the color number densities n 3 and n$ because there are only nonzero color 
chemical potentials [i 3 and m the color-superconducting phases which I had to investigate for 
this thesis. In order to fulfill color neutrality equal number densities of red, green, and blue quarks 
are necessary, 

n r = rig = rib ■ (1.58) 
By inserting this into Eqs. (|1.57p . one obtains the conditions for color neutrality 

ri3 = n$ = . 



The color neutrality conditions 



n 3 



dp 







n 8 



dp 







(1.59) 



(1.60) 



demand that the number densities of color charges 713 and ng, which can be calculated by taking 
the derivative of the pressure p with respect to the corresponding color chemical potential, are 
equal to zero. 

In QCD, color neutrality is realized dynamically due to the generation of gluon condensates 
(Aq) and (Aq) ^ [76-78]. The appearance of such condensates is equivalent to nonzero 
values of the corresponding color chemical potentials fi 3 and fi$ . 

1.4.4 Toy models of neutral normal quark matter 

In this thesis, I present the phase diagram of neutral quark matter. This will be done in Chap- 
ter O For a better understanding of the properties of neutral quark matter, it is advantageous to 
introduce some simple toy models. Therefore, some formulae of thermodynamics and statistical 
mechanics are needed [79,80]. The pressure for non- interacting massive fermions and antifermions 
at nonzero temperature T reads [81,82], 



2tt 2 



dk k 2 { In 



1 + exp - 



E-fj, 
T 



In 



1 + exp - 



E 



T 



(1.61) 



where g is the degeneracy factor, k = |k| is the momentum, E = \/k 2 + m 2 is the relativistic total 
energy m is the mass, and /1 is the chemical potential of the fermions. The first term corresponds 
to the pressure of fermions while the second one corresponds to the pressure of antifermions. By 
integration by parts, one obtains, 



6tt 2 



dfc- 



k 
~E 



4 r 



riF 



E 



T 



riF 



E 



T 



where 



riF (x) = 



1 



(1.62) 



(1.63) 



c x + 1 

is the Fermi-Dirac distribution function. At zero temperature, the Fermi-Dirac distribution func- 
tion becomes a Heaviside function, cf. Eq. (|1.1[) . 



(1.64) 



so that one has to integrate from k = up to the Fermi momentum, 



k F = \J\F 



(1.65) 
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Herewith, one obtains from Eq. (|1.62p the pressure of non-interacting massive fermions for zero 
temperature, 

fc 4 



P 



9 
6tt 2 



dfc — . 

E 



(1.66) 



Only the contribution from fermions survives. The pressure of non-interacting massive fermions 
for zero temperature can also be obtained by using the equation, 



This result can be calculated directly from the pressure (jl.6ip by using the relation, 



lim In 



-xO{-x) , 



(1.67) 



(1.68) 



or by integration by parts of Eq. Q1.66p . The integrals in Eqs. (|1.66[) and (|1.67p have an analytical 
solution, 

m 



P 



24tt 2 



j 3 ^27 ^ 4 i 

kp/j, — —m kp[L — —m in 



k F + /i 



(1.69) 



while the integrals in Eqs. (|1.61|) and (|1.62p can only be solved numerically. In the limit of 
vanishing mass, one obtains from Eq. (|1.62[) the pressure of non-interacting massless fermions and 
antifermions at nonzero temperature, 



P : 



24tt 2 



.,4 i o„ 22^2 i ' „4rp4 

15 



(1.70) 



Details how to get this result are shown in Sec. IB. II in the Appendix. In the limit of zero mass 
and zero temperature, the pressure of fermions reads, 



P 



9 4 
A 4 



24tt 2 



The number density of non-interacting fermions can be obtained by 

dp 
dfx 

which, in the case of massive fermions at nonzero temperature, leads to the result, 



9 
2vr 2 



dfc k 2 



E- n 
T 



E + n 
T 



(1.71) 



(1.72) 



(1.73) 



The first term is the contribution of fermions while the second one is the contribution of an- 
tifermions. The integral can only be solved numerically, but in the case of zero temperature, one 
obtaincs an analytical result by using Eq. (|1.64[) . 



The number density of massless fermions and antifermions at nonzero temperature reads 



so that the number density of massless fermions at zero temperature is 

9 3 
n= 6^ ' 



(1.74) 



(1.75) 



(1.76) 



From these formulae of thermodynamics and statistical mechanics one is able to construct 
simple toy models for quark matter. In the following, I present toy models of non-interacting 
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normal quark matter in neutron stars. Toy models for color-superconducting quark matter will be 
presented in Sec. 11.5.11 If a protoneutron star consists of normal quark matter, then one has to 
consider neutral /3-equilibrated quark matter at nonzero temperature. Electrons and muons are 
present in quark matter in protoneutron stars in order to make them electrically neutral. Charm, 
bottom, and top quarks as well as r leptons are too heavy so that they do not exist in the cores 
of compact stars where the quark chemical potential /i ~ 500 MeV. But in protoneutron stars, 
neutrinos are trapped. They can be treated as massless in good approximation. The pressure of 
a simple toy model of normal quark matter in protoneutron stars reads 



Pnq 



O ft 1 J P CO f 

= ±!-Y dkk 2 \\n 

rp M n CO f 




(1.77) 



where the first line in this equation is the contribution of quarks, the second line is the contribution 
of electrons and muons, and the third line is the contribution of massless electron and muon 
neutrinos to the pressure. In the third line, also the bag pressure is subtracted. The relativistic 
energies of quarks, electrons, and muons are given by Ef = (k 2 + M 2 ) 1 / 2 and Ep = (k 2 + nip) 1 / 2 , 
where Mf and nip are the masses of the quark flavors (/ = u,d,s) and leptons (/3 = e,/i). 
The chemical potentials of the quark flavors and the leptons are denoted by [if (|1.49|) or [j,p, 
respectively, where 

M/5 = Hl i3 - MQ . (1-78) 



cf. Eq. p.46[) . Because of color symmetry and spin degeneracy, the degeneracy factor of the quark 
contribution to the pressure is g = 2N C — 6. The degeneracy factor of electrons and muons is 
g = 2 because of spin degeneracy. For neutrinos the spin is always opposite the momentum and 
this is referred to as left-handed, whereas the antineutrinos are always right-handed. That is why 
the degeneracy factor of neutrinos is g = 1. The chemical potential of muon neutrinos can be set 
equal to zero which is a good approximation for matter in protoneutron stars. Electrical neutrality 
can be achieved by using Eq. (jl.54p . 



E 



E 



E Qf n f "E^ = 



(1.79) 



where 



n f 



-7T 2 , 

™ -10 



dk k 2 



Up 



E f - M/ 
T 



is the net number density of each quark flavor, and 



np 



T 



(1.80) 



np 



— dkk 2 

7T 2 Jo 



np 



E 



T 



np 



Ep + np 
T 



(1.81) 



is the number density of electrons or muons, respectively. As one expects, the condition for 
electrical neutrality of non-interacting normal quark matter in [3 equilibrium (|1.79[) is of the form 
y^ - Qini — 0, where Qi is the electric charge of the particle species i and rij its number density. In 
order to have neutral quark matter, for a given quark chemical potential u, the chemical potential 
of electric charge uq has to be found by solving for it in Eq. (|1.79p . 

A simple toy model for cold normal quark matter in neutron stars can be obtained by taking 
the limit T — ► in the above simple toy model for protoneutron stars. Also, neutrinos can be 



36 



Introduction 



neglected because in cold neutron stars, they are not trapped any more. Therefore, the pressure 
of a simple toy model of cold normal quark matter in neutron stars reads 

«.-?Ey D "^+5?gy 1 (L82) 

where fcpy = (/x 2 — M 2 ) 1 / 2 and fcp^ = (M/3 — rTi^) 1 / 2 are the Fermi momenta of the quark flavors / 
and leptons f3, respectively. Again, electrical neutrality can be achieved by solving Eq. (11.79|) for 
the chemical potential //q for a given quark chemical potential /i. But for T — > 0, the expressions 
for the number density of each quark flavor and for each lepton respectively simplify to 



i- 3 h 3 
n f = -f, nf, = -±. (1.83) 



In the limit of zero temperature, zero up and down quark masses, and massless electrons, and by 
neglecting the contribution of muons, one can solve Eq. (|1.79[> for the chemical potential of electric 
charge [83], 

M 2 

^ Q ~ 4/x" ' (L84) 
where the Taylor expansion of the strange quark Fermi momentum, 

2 j/f-2 



Mi M. 



k Fs = VJj^Ml ~ n, - -± ~ us - -i , (1.85) 

is used, and terms which are of higher order are neglected because their contributions are small. 
If the contribution of the strange quarks is neglected, one obtains, 

»Q-~\- (1-86) 

Inserting this into Eq. (|1.49[) and calculating the respective quark flavor number densities (|1.83[) 
leads to the result that there are nearly twice as many down quarks as up quarks. The simple 
results for the chemical potential of charge (11.84|) and (|1.86p are a good approximation for \i 3> M s 
as one can see by comparing it with Fig. 5 in Ref. [8]. One also realizes that strange quarks help 
neutralizing quark matter and therefore less electrons are needed. For the ideal but unrealistic 
case of zero strange quark mass, no electrons are present in neutral normal quark matter. 



1.5 Color superconductivity in neutron stars 

From the statements made in Sec. 11.4.21 one can expect strange quark matter in the cores of neu- 
tron stars where the density is so large that deconfined quark matter is able to exist. Because of the 
dominant attractive interaction in the antitriplet channel, quarks form Cooper pairs. The typical 
temperatures inside (proto)neutron stars are so low that the diquark condensate is not melted. 
That is why one expects not (only) normal strange quark matter but even color-superconducting 
strange quark matter in the cores of neutron stars. 

The color-superconducting gap affects the transport properties, e.g. conductivities and viscosi- 
ties which have an influence on the cooling rates and on the rotation period of neutron stars. 
It also modifies the thermodynamic properties, e.g. the specific heat and the equation of state 
which have an influence on the mass-radius relation of color-superconducting neutron stars [9] . In 
Ref. [8], the effect of color superconductivity on the mass and the radius of compact stars made 
of pure quark matter is investigated. The authors confirmed the result of Ref. [86] that color 
superconductivity does not alter the mass and the radius of quark stars, if the diquark-coupling 
constant is chosen to reproduce vacuum properties such as the pion-decay constant. The reason is 
that color superconductivity in neutral quark matter has a tiny effect on the equation of state, cf. 
Figs. 5.9 and 5.10 in Ref. [7]. The color-superconducting gap has significant effects on the mass 
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and radius of quark stars only if the diquark-coupling constant is artifically increased whereby the 
value of the gap itself artificially increases. For gaps on the order of 300 MeV, the mass and radius 
of quark stars are approximately twice as large as for normal-conducting quark stars so that such 
quark stars are of the same size and mass as ordinary neutron stars. Therefore, it is impossible 
to decide whether a compact star consists of normal conducting or color-superconducting quark 
matter, or simply of hadronic matter. 




Figure 1.10: The low-energy part of the dispersion relations, i.e. the quasiparticle energies as a function 
of the quark momentum. Left panel: dispersion relations of quark quasiparticles in the 
2SC phase (solid lines) and the normal quark phase (dashed lines) . Right panel: dispersion 
relations of quark quasiparticles in the g2SC phase (solid lines) and the normal quark phase 
(dashed lines). The dispersion relation with the lowest energy in the g2SC phase shows a 
gapless mode [9,84,85]. The value of the averaged chemical potential of the up and down 
quarks which pair is denoted by p, and one half of the difference of their chemical potentials 
is denoted by 8[i. In the massless case, the chemical potentials of the free up and down 
quarks are replaced by fj^ = p, ± [{S/j,) 2 — A 2 ] 1//2 . 



In some cases, color superconductivity is accompanied by baryon superfluidity or the electro- 
magnetic Meissner effect. Baryon superfluidity causes rotational vortices while the electromagnetic 
Meissner effect entails magnetic flux tubes in the cores of neutron stars. 

At large strange quark masses, neutral two-flavor quark matter in f3 equilibrium can have an- 
other rather unusual ground state called the gapless two-flavor color superconductor (g2SC) [84]. 
While the symmetry in the g2SC ground state is the same as that in the conventional 2SC phase, 
the spectrum of the fermionic quasiparticles is different, see Fig. ll.lOl The g2SC phase appears at 
intermediate values of the diquark coupling constant while the 2SC phase appears in the regime of 
strong diquark coupling. Gapless modes are created if the mismatch between the Fermi momenta 
of the quarks which pair becomes large. In the case of the g2SC phase, 5fi > A where <5/i = /x e /2. 
The existence of gapless color-superconducting phases was confirmed in Refs. [8,87,88], and gen- 
eralized to nonzero temperatures in Refs. [85,89]. It is also shown that a gapless CFL (gCFL) 
phase appears in neutral strange quark matter [90,91]. But gapless color-superconducting quark 
phases are unstable in some regions of the phase diagram of neutral quark matter because of 
chromomagnetic instabilities [92] so that another phase will be the preferred state. Chromo- 
magnetic instabilities occur even in regular color-superconducting quark phases. The author of 
Ref. [93] shows that chromomagnetic instabilities occur only at low temperatures in neutral color- 
superconducting quark matter. The author of Ref. [94] points out that the instabilities might 
be caused by using BCS theory in mean-field approximation, where phase fluctuations have been 
neglected. With the increase of the mismatch between the Fermi surfaces of paired fermions, phase 
fluctuations play more and more an important role, and soften the superconductor. Strong phase 
fluctuations will eventually quantum-disorder the superconducting state, and turn the system into 
a phase-decoherent pseudogap state. By using an effective theory of the CFL state, the author of 
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Rcf. [95] demonstrates that the chromomagnetic instability is resolved by the formation of an inho- 
mogeneous meson condensate. The authors of Ref. [96] describe a new phase in neutral two-flavor 
quark matter within the framework of the Ginzburg-Landau approach, in which gluonic degrees 
of freedom play a crucial role. They call it a gluonic phase. In this phase, gluon condensates 
cure a chromomagnetic instability in the 2SC solution and lead to spontaneous breakdown of the 
color gauge symmetry, the [£/(l) em ] and the rotational 5*0(3) group. In other words, the gluonic 
phase describes an anisotropic medium in which color and electric superconductivity coexist. In 
Ref. [97], it was suggested that the chromomagnetic instability in gapless color-superconducting 
phases indicates the formation of the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) phase [98] which 
is discussed in Ref. [99] in the context of quark matter. Other possibilities could be the formation 
of a spin-one color-superconducting quark phase, a mixed phase, or a completely new state. The 
authors of Ref. [100] suggest that a mixed phase composed of the 2SC phase and the normal 
quark phase may be more favored if the surface tension is sufficiently small [101]. The authors 
of Ref. [102] suggest a single first-order phase transition between CFL and nuclear matter. Such 
a transition, in space, could take place either through a mixed phase region or at a single sharp 
interface with electron- free CFL and electron-rich nuclear matter in stable contact. The authors 
of Ref. [102] constructed a model for such an interface. 



1.5.1 Toy models of neutral color-superconducting quark matter 

In this subsection, I show some toy models for color-superconducting quark matter. In these toy 
models, which are valid for zero temperature, up and down quarks are treated as massless, and 
the strange quark mass is incorporated via a shift in the strange quark chemical potential which 
is a good approximation. Neutrinos are not present in cold quark matter in neutron stars so that 
fi e = —(J.Q- Also, muons are not taken into account in the following toy models but electrons which 
are treated as massless for simplicity. In the equations which satisfy the neutrality conditions of 
these toy models, terms of higher order will be neglected. The pressure of the toy model for 
color-superconducting quark matter in the 2SC phase reads [83] , 
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where the first line is the contribution of gapped quasiparticles, the second line is the contribution 
of free blue up and free blue down quarks as well as free strange quarks, and the third line is the 
contribution of massless electrons, the pressure correction due to the four gapped quasiparticles, 
and the bag pressure. The chemical potential for each quark color and flavor is given by Eq. (|1.52j) . 
and their Fermi momenta by 

(k F )) = yV}) 2 ~ Mf ■ (1-88) 

Up and down quarks are treated as massless in this toy model, and the strange quark Fermi 
momenta for each color will be calculated by the simplified expression, 



(k F )l = V(/4) 2 -M2 * /4 - * Ml - . (1-89) 
as mentioned above. All quarks which pair have the same common (averaged) Fermi momenta, 

k< F — EE k Fl EE (k F ) T u = (k F f d = k F2 EE (k F ) r d = (k F )l = fl+ ~/ZQ + ^= M8 . (1.90) 
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In the 2SC phase, the Fermi momenta k Fl and kp 2 are equal. Therefore, only one common Fermi 
momentum fc™ mmon i s used. The neutrality conditions (|1.54|) and (|1 .60|) require [83], 

M 2 

HQ — ~~2/j ' ^ 3 = ' M8 -°- (L91) 
In the case without strange quarks, one obtains 

MQ-~f, M3 = 0, M8-0. (1.92) 

These approximations are in good agreement with the exact results for /j, ^> M s , cf. Figs. 5 and 6 
in Ref. [8]. In this toy model, the pressure difference of the 2SC phase to the normal quark phase 
reads [83] 

16/i 2 A 2 - Mf 
16tt 2 

so that the 2SC phase is preferred to the normal quark phase when 



S P2SC ^ P 1g 2 s (1.93) 



M 2 

A > -f- ■ (1.94) 



Mi 

The pressure of the toy model for the CFL phase reads 

*^ = ^££ / ^(tf-litf + ^ + ^-fl, (1.95) 

i=r f=u J ° 

where the first term is the contribution of the nine gapped quarks, the second term is the con- 
tribution of electrons, the third term is the correction to the pressure due to the gap, and the 
fourth term is the bag pressure. As in the 2SC phase, the quarks which pair have the same Fermi 
momenta, 
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(1.96a) 



t+ 6^ Q + 2^1 M8 ' (1.96b) 

111 M 2 

'' + 6 (lQ + l fl3 -^ S - 4^' (L96C) 

1 1 1 M s 2 . 

! - 3^-4^-^^- 4^- (L96d) 

In the CFL phase, the neutrality condition for fi 3 , Eq. (|1.60p . requires that 

= -MQ • (1-97) 
It is very useful to know that in the CFL phase the following relation is valid: 

nQ - n 3 \= n % + n e = , (1.98) 

V3 

where n e is the number density of electrons. Because of the neutrality conditions, tiq, n 3 , and 
rig have to be equal to zero so that one concludes from Eq. (|1.98[) that also n e and with it fie 
have to be zero which is in agreement with the arguments made in Ref. [20] . It means that in the 
CFL phase there are neither electrons nor muons nor t leptons allowed, otherwise the neutrality 
conditions will be violated. That is why one obtains the simple result 



HQ = M3 = 



(1.99) 
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In order to satisfy the neutrality condition for /xg, Eq. (|1.60p . one gets, 



(1.100) 



Ms - - 



In this toy model, the pressure difference of the CFL phase to the normal quark phase reads [83, 



This is the same requirement as for the 2SC phase, and because the pressure difference of the CFL 
phase to normal quark matter is three times larger than the pressure difference of the 2SC phase 
to normal quark matter, the authors of Ref. [83] claim that the 2SC phase is absent and the CFL 
phase is the preferred state in compact stars. In Chapter ® I shall show that this statement is 
not correct in general so that these simple toy models cannot replace a thorough analysis of the 
preferred quark phases by an NJL model [14, 103-106] for example. That is why it is so important 
to construct the phase diagram of neutral quark matter with a more precise model which is done 
in Chapter [2] 




(1.101) 




(1.102) 



Chapter 2 



The phase diagram of neutral 
quark matter 

At sufficiently high densities and sufficiently low temperatures quark matter is a color supercon- 
ductor [10-12]. This conclusion follows naturally from arguments similar to those employed in 
the case of ordinary low-temperature superconductivity in metals and alloys [13]. Of course, the 
case of quark matter is more complicated because quarks, unlike electrons, come in various flavors 
(e.g., up, down, and strange) and carry non-Abelian color charges. This phenomenon was stud- 
ied in detail by various authors [15, 18, 19,32,33,35,36,99, 107-109]. Many different phases were 
discovered, and recent studies [84,85,87-91,110-119] show that even more new phases exist. (For 
reviews and lectures on color superconductivity see Refs. [3,9,68-73].) 

In nature, the most likely place where color superconductivity may occur is the interior of 
neutron stars. Therefore, it is of great importance to study the phases of dense matter under the 
conditions that are typical for the interior of stars. For example, one should appreciate that matter 
in the bulk of a star is neutral and /3-equilibrated. By making use of rather general arguments, it 
was suggested in Ref. [83] that such conditions favor the CFL phase and disfavor the 2SC phase. 
In trying to refine the validity of this conclusion, it was recently realized that, depending on the 
value of the constituent (medium modified) strange quark mass, the ground state of neutral and 
/3-equilibrated dense quark matter may be different from the CFL phase. In particular, the g2SC 
phase [84, 85] is likely to be the ground state in the case of a large strange quark mass. On the 
other hand, in the case of a moderately large strange quark mass, the CFL and gCFL phases [90] 
are favored. At nonzero temperature, T^0, some other phases were proposed as well [111]. 

I note that the analysis in this thesis is restricted to locally neutral phases only. This automat- 
ically excludes, for example, crystalline [99] and mixed [100, 101, 120] phases. Also, in the mean 
field approximation utilized here, I cannot get any phases with meson condensates [28-30]. 

In this chapter, I present the phase diagram of neutral quark matter. In Sec. 12.11 I show 
the phase diagram of massless neutral three-flavor quark matter. In Sec. 12. 2[ quark masses are 
treated self-consistently within the framework of a three-flavor NJL model [106], and the phase 
diagram with a self-consistent treatment of quark masses is presented. In Sec. l2.3[ this NJL model 
is extended to nonzero neutrino chemical potentials, and the influence of neutrinos on the phase 
diagram is discussed. 

I use the following conventions: I calculate in natural units, h = c = feg = lj and utilize the 
Dirac definition of the 7-matrices which are shown in Sec. IA.2l in the Appendix. Latin indices run 
from one to three while Greek indices run from zero to three. Four-vectors are denoted by capital 
Latin letters while three- vectors are written in the bold upright font. The space-time vector is 
defined as A M = (x°, x 1 , x 2 , x 3 ) — (t,x), the four-momentum vector is denoted as = (fc°,k), 
and the metric tensor is given by g^ u = diag (1, — 1, — 1, — 1). Absolute values of vectors are 
denoted by italic Latin letters, e.g. k = |k|. The direction of a vector is indicated by the hat 
symbol, e.g. k = k/fc. I use the imaginary-time formalism, i.e., the space-time integration is 
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defined as f x = j Q dr J v d 3 x, where r is the Euclidean time coordinate and V the three- volume 
of the system. The delta function is defined as 5^ (X — Y) = 5 (rx — t~y) 5^ (x — y). Energy- 
momentum sums are written as: T/V ~^2 K = J d 3 k/ (2n 3 ) — T/(2n 2 ) J2 n J dkk 2 , where the 
sum runs over the Matsubara frequencies oj n = 2nirT = iko for bosons and tu n = (2n + 1) ttT = iko 
for fermions, respectively. 

2.1 The phase diagram of massless quarks 

In this section, I study the phase diagram of massless neutral three-flavor quark matter at zero 
and nonzero temperature as a function of the strange quark mass as well as a function of the quark 
chemical potential. I start with the QCD Lagrangian density in order to derive the QCD grand 
partition function. Then, the Cornwall- Jackiw-Tomboulis (CJT) formalism [121] is used to calcu- 
late the QCD effective action. The gluon-exchange interaction between quarks is approximated 
by a point-like four-fermion coupling. From the effective action of quarks, I derive the pressure 
of color-superconducting quark matter. In order to allow for the most general ground state, I 
employ a nine-parameter ansatz for the gap matrix [110]. The effects of the strange quark mass 
are incorporated by a shift of the chemical potential of strange quarks, /i* — > /x* — mj:/(2/i), where 
i = r, g, b is the color index, m s is the strange quark mass, and /i is the quark chemical potential. 
This shift reflects the reduction of the Fermi momenta of strange quarks due to their mass. Such 
an approach is certainly reliable at small values of the strange quark mass. I assume that it is 
also qualitatively correct at large values of the strange quark mass. In order to draw the phase 
diagram of massless neutral three-flavor quark matter, I solve nine gap equations together with 
the conditions of electric and color charge neutrality. 

2.1.1 Quantum chromo dynamics 

The QCD Lagrangian density is given by 

L = $ (i0 - m) $ - ~i^,*r + £gau g c . (2.1) 

For N c colors and Nf flavors, ip is the 4iV c jV/-dimensional spinor of quark fields, tp = V^7o is the 
Dirac conjugate spinor, J3 = j^D^, where 7 M are the Dirac matrices, and m is the quark-mass 
matrix. The covariant derivative is defined as 

= a M - igA«T a , (2.2) 

where g is the strong-coupling constant, A a ^ are the gluon fields, and T a are the generators of the 
[SU (N c ) c ] group. The gluon field-strength tensor is defined as 

Fr - - + gUcA^Al , (2.3) 

where / a & c are the structure constants of the [SU(N C ) C ] group. The first term in the QCD La- 
grangian density (|2.ip corresponds to interacting quarks while the second one corresponds to the 
gluons. The term £ gaU gc comprises gauge fixing terms and the contribution from Faddeev-Popov 
ghosts. Up to irrelevant constants, the grand partition function of QCD is given by [122, 123] 

Z = J V^VA" exp {/ 1>, A] } , (2.4) 

where 

/ V>, A] = [ (£ + fjM + H Q M Q + llaNa) , (2.5) 
J X 

is the QCD action. The conserved quantities, 



(2.6) 
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are the quark number density operator, the operator of electric charge density of the quarks, and 
the operators of color charge densities of the quarks, respectively. The QCD action can be split 
into the part of interacting quarks, into the part of gluons, and into the part of gauge fixing terms 
and Faddeev-Popov ghosts, 

I V», A] = 1^ $, A] + I A [A] + 7 gaugc [A] , (2.7) 

where 

1^ if),A] = ( $ (X) (X, Y) yj (Y) , (2.8a) 

JX,Y 

Ia[A] = -\J F; v (X) Fr (X) , (2.8b) 

-^gauge [A] — / -^gauge ■ (2.8c) 
JX 

The inverse tree-level propagator for quarks and charge-conjugate quarks, respectively, is given by 

[G+]- 1 (X, Y) = (%J3 X + A 7 o - rh) 6& {X — Y) , (2.9a) 
(X, Y) = {iP c x - Ato - m) (X — Y) , (2.9b) 

see Sec. IB.2I in the Appendix. The matrix of quark chemical potentials is denoted by /2 and is de- 
fined by Eq. (|1.51[) . The covariant derivative in the inverse tree- level propagator of quarks is given 
by Eq. (12. 2|l while the charge- conjugate covariant derivative in the inverse tree-level propagator of 
charge-conjugate quarks reads, 

JnZ = dr+igApZ. (2.10) 

By using the charge-conjugate spinors (|A.14|) and by introducing the Nambu-Gorkov basis with 
the 2 • 4iV c jV/-dimensional quark spinors, 

one can rewrite the QCD action of interacting quarks, 

1^ =\ [ *(X)S^(X,Y)>S>(Y) , (2.12) 

L JX,Y 

where 

is the inverse tree-level propagator for Nambu-Gorkov quarks. 

By adding a bilocal source term to the QCD action (|2.7p , one obtains [3,37], 

I [V,V,A,tC] =/[*,*, A] + - I ®(X)fC{X,Y)V(Y) , (2.14) 

2 Jx.Y 



(!+!-)■ (2-15) 



where 

The four entries of K, are not independent. Due to charge-conjugation invariance of the action, 
(7~ = C (c + ) T C _1 , and since the action has to be real-valued, ip~ = 70 (<^ + )^ 7o- In the presence 
of an external source IC, the grand partition function of QCD, up to irrelevant constants, reads, 

Z[K\ = JvVDyVAlexp{l[%*,A,lC]} . (2.16) 
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The QCD effective action can be derived from the grand partition function of QCD (12.16|) by 
using the CJT formalism [3,8,37,121,124-126], 

T A,S,V] = / - ^TrlnV- 1 - ^Tr (V^V - l) 

+ iTrln5- 1 + iTr(5 - 1 5-l)+r 2 [^,*,A,5,I?] . (2.17) 

The quantities T> and S are the gluon and quark propagator, respectively. The inverse treedevel 
quark propagator Sq 1 was introduced in Eq. (|2.13p . Correspondingly, "Dq 1 is the inverse tree- 
level gluon propagator. The traces run over space-time, Nambu-Gorkov, color-flavor, and Dirac 
indices. The factor ^ in front of the fermionic one-loop terms compensates the doubling of degrees 
of freedom in the Nambu-Gorkov basis. The functional r 2 is the sum of all two-particle irreducible 
(2PI) diagrams. It is impossible to evaluate all 2PI diagrams exactly. However, the advantage 
of the QCD effective action (|2.17j) is that truncating the sum r 2 after a finite number of terms 
still provides a well-defined many-body approximation. Later on, the gluon-exchange interaction 
between quarks is approximated by a point-like four-fermion coupling. This effectively removes 
dynamical gluon degrees of freedom, such that one does not need to worry about gauge fixing or 
possible ghost degrees of freedom. Therefore, the latter are already omitted in Eq. (|2.17p . The 
stationary points of the effective action (|2.17p determine the expectation values of the one- and 
two-point functions, 

5T 5T 5T 5T 5T , 

M = ' M = ' ^ = ' ^ = ' ss=°- (2 ' 18) 

The first two equations yield the Dirac equation for the quark fields W and in the presence 
of the gluon field A® . The third equation is the Yang-Mills equation for the gluon field, 



(2.19) 



where D® b = d v 8 ah — gj ahc A C V {X) is the covariant derivative in the adjoint representation. The 
first two terms on the right-hand side are the contributions from gluon and fermion tadpoles [76] . 

The functional derivative with respect to A® acting on the trace is nontrivial because of the 
dependence of the inverse tree-level propagators Vq 1 and Sq 1 on the gluon field, cf. Eq. (|2.13|) . 
The last term is nonzero if r 2 contains 2PI diagrams with an explicit dependence on A® . As shown 
in Ref. [76], the solution of the Yang-Mills equation in the 2SC phase is a constant background 
field A a ^ ~ g^oS aS - This background field acts like a color chemical potential us and provides the 
color-charge neutrality of the 2SC phase [76]. In this manner, also the color chemical potential 
/j,3 is present in three-flavor color superconductors where the color chemical potentials /i3 and fis 
ensure color neutrality. Later on, I shall remove the gluon degrees of freedom by approximating 
the non-local gluon exchange with a point-like four-fermion coupling. The constant background 
field A a ^ then disappears from the treatment, and the color chemical potentials 113 and fj,% assume 
the role of the background field to ensure color neutrality. 

The fourth equation (|2. 18|) is the Dyson-Schwinger equation for the gluon propagator, 

V-^(X,Y)=V^Z(X,Y)+U^(X,Y) , (2.20) 

where 

KliX,Y) = -2 w ^- ) (2.21) 

is the gluon self-energy. Since I shall approximate the gluon exchange by a four-fermion coupling, 
I do not need to solve the Dyson-Schwinger equation for the gluon propagator. 

The fifth equation (|2.18p is the Dyson-Schwinger equation for the quark propagator, 

S- 1 (X, Y) = Sv 1 (X, Y) + S (X, Y) , (2.22) 



2.1 The phase diagram of massless quarks 



45 



where 

^ x ^- 2 wh^) (2 - 23) 

is the quark self-energy. According to their definition, the self-energies (|2.21[) and (|2.23p are 
obtained from the set of 2PI diagrams by opening one internal line. The quark self-energy is a 
2 x 2-matrix in Nambu-Gorkov space, 



S=^ +E -J. ( 2 - 24 ) 

where the quantities "E^ are the regular self-energies, and $> ± are the anomalous self-energies. 
For $ , also the term gap matrices is used. The gap matrices in connection with the quark 
self-energy (12.23)) yield the so-called gap equations. By solving these gap equations, one obtains 
the gap parameters. In space-time, S + (X, Y) [E _ (X, Y)] has a quark [charge- conjugate quark] 
entering at X and another quark [charge-conjugate quark] emerging at Y. The anomalous self- 
energies have to be interpreted as follows: a quark [charge-conjugate quark] enters <E> + (X, Y) 
[<£> - (X, Y)] at X and, in contrast to the regular self-energies, here, a charge-conjugate quark [an 
ordinary quark] emerges at Y. This is typical for systems with a fermion-fermion condensate in 
the ground state [127]. The self-energies 4> ± symbolize this condensate. This is why the crucial 
quantities regarding color superconductivity are the gap matrices $ . A nonzero value of $ is 
equivalent to Cooper pairing, or, in other words, to a nonvanishing diquark expectation value. 
The self-energies in Eq. (|2.24| are related in the same way as the bilocal sources in Eq. (|2.15[) . 

S- =C(S+) T C- 1 , *~ = 7o ($+) t 7 o • (2.25) 

The quark propagator in Nambu-Gorkov space can be determined from the Dyson-Schwinger 
equation 1)2.22|) . see Sec. IB.3I in the Appendix, 

g + 5- 



where 

g ± 



S = ( 5+ 0- ) > ( 2 - 26 ) 



{[Gt]- 1 + ^ ± -^{[Q^]- 1 + ^ ¥ ) \ (2-27) 

is the propagator for quasiquarks or charge-conjugate quasiquarks, respectively, and 

s ± = - {[g^r 1 + $±g± = -g*$± ([gtr 1 + z ± y 1 , (2.28) 

are the anomalous propagators. These anomalous propagators are typical for superconducting 
systems [127] and account for the possibility that in the presence of a Cooper-pair condensate, 
symbolized by < i> ± , a fermion can always be absorbed in the condensate, while its charge-conjugate 
counterpart is emitted from the condensate and continues to propagate. 

The QCD pressure is, up to a prefactor T/V, equal to the QCD effective action (|2.17p deter- 
mined at the stationary points (|2.18|) which is denoted by T* , 

P=^T\ (2.29) 
where T is the temperature and V the volume of the system. 

2.1.2 The effective action of quarks 

In this section, I investigate the phase diagram of massless neutral three-flavor quark matter. The 
quark spinor has the following color-flavor structure: 

^=(il>*,tt,M,M,tf>M>4l4ltf) T ■ (2-30) 
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Since I approximate the gluon-exchange interaction between quarks by a point-like four-fermion 
coupling, I only need to consider the contributions of quarks to the QCD effective action (|2.17[l . 

T [S] = ilr In S- 1 + ilr (S^S - l) + T 2 [S] . (2.31) 

Here, the tree- level quark propagator Sq 1 (|2.13|) . which occurs in the QCD effective action (|2.17j) . 
is replaced by the quark propagator, 

So'= ( [G °J [G -°]-i) > ( 2 - 32 ) 

where 

[G*]" 1 (X,Y) ee (i? x ± A70) <5 (4) (X - Y) , (2.33) 

is the massless inverse Dirac propagator for quarks and charge-conjugate quarks, respectively, in 
which the constant background field A a ^ disappears from the treatment, and the color chemical 
potentials /13 and ^ assume the role of the background field to ensure color neutrality. The quark 
chemical potential matrix in color-flavor space is defined as 

ft, = diag(^,/iS,^,^,^,^,^,^,^) , (2.34) 

where the chemical potential of each quark is given by Eq. (|1.52p because quark matter inside 
neutron stars is in (3 equilibrium. 

Since I shall present the phase diagram of massless neutral three-flavor quark matter in this 
section, the mass term of the inverse Dirac propagators (|2.33[) is omitted. At sufficiently large 
quark chemical potential, there is no need to take into account the small up and down quark 
masses because the dynamical effect of such masses around the quark Fermi surfaces is negligible. 
Of course, the situation with the strange quark is different because its mass is not very small 
as compared to the quark chemical potential /1. The most important effect of a nonzero strange 
quark mass is, however, a shift of the strange quark chemical potential due to the reduction of the 
Fermi momentum, 

(k F y s = vuw^l - ^ - ^- - ^ - if - ( 2 - 35 ) 

cf. Eq. p.89p . Here, I have approximated \i % s by /i in the denominator. Quantitatively, this does 
not make a big difference. 



Figure 2.1: Left panel: the sunset-type diagram. Right panel: the double-bubble diagram. 



For the sum of all 2PI diagrams I only include the sunset-type diagram, which is shown in the 
left panel of Fig. O 

r 2 [S] = - £ / Tr [T^S (X, Y) VIS (Y. X)] D% (X, Y) . (2.36) 

4 JX,Y 
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The trace runs over Nambu-Gorkov, color-flavor, and Dirac indices. The Nambu-Gorkov vertex is 
defined as 

The stationary point of the effective action of quarks (|2.3ip , 

§ = , (2.38) 
is the Dyson-Schwinger equation for the quark propagator, 

S- 1 {X, Y) = S^ 1 (X, Y) + E (X, Y) , (2.39) 

where 

£ Y) = 2 6s S { y* x) = -9 2 ^ a S (X, Y) T» b Df v (y, X) (2.40) 

is the quark self-energy. The effective action at the stationary point which is determined by the 
Dyson-Schwinger equation for the quark propagator (12.39|) reads, 

T* = ^Tr Ills'- 1 - -Tr(ES') . (2.41) 

For translationally invariant systems, it is advantageous to work in energy-momentum space 
instead of in space-time, 

S- X (X,Y) = ^^e-^x-^S- 1 (K) , (2.42a) 

K 

S(X,Y) = ^J2 e ~ iK(X ~ Y)s ( K ) > ( 2 - 42b ) 
K 

Df u {X, Y) = ^J2 e- lK(X - Y) K ( R ) > ( 2 - 42c ) 

K 



E (K) = / e lKZ E(Z) , (2.42d) 
Jz 

where I assumed translational invariance of propagators and self-energies, Z = X — Y . In energy- 
momentum space, the Dyson-Schwinger equation for the quarks reads, 

S- 1 (K) = So 1 (K) + E (K) . (2.43) 

The inverse free quark propagator in energy-momentum space is given by 



(2.44) 



where the massless inverse Dirac propagator for quarks and charge-conjugate quarks, respectively, 
reads, 

[G^]- 1 (^)=7o(fco±A)-7-k. (2.45) 
The quark self-energy is obtained by 

E= (£+ D = - g2 v E T « s (o) T * D t> > ( 2 - 46 ) 

where 



$+ I » V 

7 Q 



Q Q 
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are the gap equations which are independent of K because I approximate the gluon-exchange 
interaction between quarks by a point-like four-fermion coupling. The inverse quark propagator is 



S- 1 (K) = 



[Gt] 



[Go] 



(2.48) 



The regular self-energies play an important role in the dynamics of chiral symmetry breaking, but 
they are of less importance in color-superconducting quark matter. (The effect of the regular self- 
energies was studied in Ref. [107].) Therefore, the regular self-energies in Eq. (|2.46[) are omitted. 
The quark self-energy (|2.46[) contains the Feynman gauged gluon propagator, see Sec. IB. 41 in the 
Appendix, 

(2.49) 



D 



ah 



tab 9^v 
A 2 



which represents the gluon-exchange interaction between quarks by a point-like four-fermion cou- 
pling. In this approximation, the sunset-type diagram becomes a double-bubble diagram, which 
is shown in the right panel of Fig. 12.11 

The quark propagator may be obtained by inverting Eq. (|2.48|) . 



S(K) 



(G+ 3- 



where 



G ± (K) = {[G±]- 1 -**G%*±} 1 
is the propagator for quasiquarks or charge-conjugate quasiquarks, respectively, and 

3 ± (K) = -G ¥ $ ± G ± = -G T $ ± Gj , 



(2.50) 



(2.51) 



(2.52) 



are the anomalous propagators. 

I want to study the most general ansatz of the gap matrix for the CFL phase in color-flavor 
space. Therefore, I use the following nine-parameter ansatz [110]: 
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(2.53) 



This ansatz for the gap matrix is indeed the most general one for the CFL phase because it 
can be obtained with some modifications from the antitriplct and sextet gap ansatz (|1.20p . The 
difference between the gap ansatz (|1.20[) and the gap matrix (|2.53[) is that the gap ansatz (|2.53[) 
has different entries for all nonzero color-flavor elements. The wavefunction requirements arc still 
fulfilled because the gap matrix (|2.53p is symmetric in color-flavor space. 

The effective action of quarks at the stationary point (|2.41|) in energy-momentum space reads, 



r* = i5>dct(^)-i£Tr(E5), 

K V / K 



(2.54) 



where I used Eqs. (|2.42[) , the derivation of the Fourier transformed kinetic part of the QCD grand 
partition function which is shown in Sec. IB.5l in the Appendix, and the relation lndet A = TrlnA 
which is proven in Sec. IB.6l in the Appendix. 
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2.1.3 Propagators and self-energies in projector representation 

I only represented the color-flavor structure in the gap ansatz (|2.53p . Together with the Dirac 
part, the complete gap matrices read [128], 

* + = E <^c(k) , $ - = £ 4>lVZm , (2.55) 

c,e c,e 

where <f% is only the color-flavor part of the gap ansatz (|2.53[) , and 

P c e (k) = i(l + c 7 5)(l + e7o7-k) (2.56) 

are the energy-chirality projectors, see Eq. (|A.16c[) . The gap parameters </>c//' are real- valued 
numbers. It will be very advantageous in later calculations if the (inverse) massless Dirac propa- 
gators are also written in terms of the energy-chirality projectors, 

[Gjr 1 = ioY,^r 1 v!(k) = Y,[G^}- 1 v!.^)io, (2.57a) 

c.e c.e 

Go = E G ° V c ( k ) 7o - 7o ^ G±V e c (k) . (2.57b) 

c.e c,e 

Let the matrix of quark chemical potentials // and the color-flavor part of the gap matrix <fi be 
arbitrary for the moment. I separated the Dirac part and the color- flavor part, 

[G^]- 1 = fc ±M- ek , [Go]- 1 = k ±fl + ek , (2.58) 

of the (inverse) massless Dirac propagators in Eq. (|2.57j) . Note that the color- flavor matrices in 
Eq. (I2.58|) depend on e. I did not show an extra index e in the above equations for simplicity. I 
shall omit the indices c and e at the color-flavor matrices of propagators and self-energies. 
The (inverse) regular propagator reads, 

[G + ]- 1 =7oE^ + r 1 ^ e (k) , G+=^G+P c e (k) 7 o, [G+]- 1 = [G+]- 1 - 4>G^4> ■ (2.59) 

c,e c,e 

The anomalous propagator can be written as 

S+ = Y, ^V-t (k) , 5 + = -Go4>G+ . (2.60) 

c,e 

2.1.4 The potential part of the effective action of quarks 

The potential part of the effective action of quarks (|2.54[l . 

r; ot ^-J]TTr(£S) , (2.61) 

K 

can be simplified. After performing the Nambu-Gorkov trace, one obtains for the potential part 
of the effective action of quarks, 

r ;ot = 4E Tr ( $ ~ s++<i>+s ~) • ( 2 - 62 ) 

K 

With the anomalous propagators (|2.52[) and the relation Tr (AB) = Tr (BA), one can easily check 
that 

Tr ($~S+) = -Tr ($-G(7$ + G + ) = -Tr ($+G+$"G ) = Tr ($+S _ ) . (2.63) 



50 



The phase diagram of neutral quark matter 



Therefore, the potential part of the effective action of quarks can be simplified to 

Kot - -\ E Tr = - \ E E Tr V i+ vzl (k) ] = - ^ E E Tr (^ + ) ' ( 2 - 64 ) 



K c,e 



K c,e 



where the Dirac trace has been performed in the last step. In the following, I shall further simplify 
the potential part of the effective action of quarks (|2.64p : I insert the anomalous propagator (|2.60|) 
into the gap equation (|2.47[) and multiply it from the left side with the energy-chirality projectors 
in order to project onto (jf c - Thus, the sums over c and e vanish. By performing the Dirac trace 
one obtains, 



,T 



Q e> 



V e c (h)^VZ e c , (q) 7 
V c A e (k) 7 ^A e ' (q) 7 " 



D ab 

{IV 



(2.65) 



In this expression, I explicitly show that the color-flavor part of the gap matrix carries the indices 
c and e which arise from the energy and chirality projections. In the second line, the color-flavor 
part of the Feynman gauged gluon propagator has been performed. The Dirac trace is computed 
in Sec. IB. 71 in the Appendix. One obtains 



4r = -¥^EE^ + )/F(Q)^'' v , 



A 2 V 



(2.66) 



Q e' 



where Tf i = {Tf) T . As one can see, the gap equations do not depend on the indices c and e 
because they do not appear any longer on the right-hand side of this equation. After summation 
over a one obtains [19] 



J ff> 



which can be rewritten as 



A2 V 



Q e' 



-(H+)>> (Q) - ~(E + )f f , (Q) 



(2.67) 



3 A 2 V 



(2.68) 



where only one term is left on the right-hand side. This equation can be transformed into 

3 A 2 V 
9 2 



30' 



(2.69) 



K e 



where I renamed Q in K and e' in e. I now want to write the right-hand side of this equation in 
a nicer form. Therefore, I use the relation, 



T 31 6 

a 

Multiplying by six and rearranging yields, 

a 

This expression can be inserted into Eq. (12.691) . 

3A 2 V 



(2.70) 



(2.71) 



EE^w 



K e 



4 g 2 T 



3^T Q T 0T a 



(2.72) 
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By inserting this result into the potential part of the effective action of quarks (|2.64p one obtains 



r = ^Tr 
pot A g 2 T 




(2.73) 



where the sum over c has been performed which yields an extra factor of two. 
2.1.5 The kinetic part of the effective action of quarks 

The inverse quark propagator S~ l (|2.48[) is a 72 x 72-matrix and consists of the following subspace 
structure: it is a 2 x 2-matrix in Nambu-Gorkov space, a 9 x 9-matrix in color-flavor space, and 
a 4 x 4-matrix in Dirac space. In order to calculate the kinetic (first) part of the effective action 
in Eq. (|2.54p . I transform the inverse quark propagator into block-diagonal form. This is easily 
achieved by changing the order of rows and columns in color-flavor and Nambu-Gorkov space. By 
such a transformation, the absolute value of a determinant does not change. Only the sign of the 
determinant is changed if there is an odd number of exchanges of rows and columns. In the case 
mentioned above, no corrections have to be made to the result of the determinant because there 
is an even number of exchanges of rows and columns. There are six 2 x 2-blocks, 

o-i _ ( [Go+r 1 [*5r A s -i_( [GtT 1 [*S 



t ~9\ 
T sJ 



* 3 "I M + [Goi]- 1 ) ' \ m b u} + [GoT 1 i ' 

and one 6 x 6-block in the block-diagonal structure of the inverse quark propagator (|2.48[) in 
color-flavor and Nambu-Gorkov space. This 6 x 6-block is a 2 x 2-matrix in Nambu-Gorkov space, 

g-l = ( *(3x3) \ (2 75) 

\ *(3x3) L^O J(3x3) / 

where the Dirac propagators as well as the gap matrices, 

[G±]^ 3) =di ag ([G±Q-\[G±^\[G±^ , $f 3x3) = KT Kl ± , (2-76) 

are 3 x 3-matrices in color-flavor space. 

Since det [diag (A, B)\ = det A ■ det B [129], where A and B are quadratic matrices, the kinetic 
part of the effective action of quarks (|2.54|) is, 

S~ l \ 1^, , {S7 1 



n^\T.^{^r)=\Y.T.^{^) • (2-77) 
In order to simplify the determinants, I make use of the Gauss-elimination procedure [129], 

det (c D ) " det [ D ( A ~ BD ~ lc )] , (2.78) 



where A, B, C, and D arc quadratic matrices and D is invertiblc. Hereby, all 2 x 2-blocks in 
color-flavor and Nambu-Gorkov space become 1 x 1-blocks in color-flavor space, while the original 
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6 x 6-block in color-flavor and Nambu-Gorkov space becomes a 3 x 3-block in color-flavor space 
so that the kinetic part of the effective action reads, 



r* = - V 

km r, 



K 



In det 



In det 



In det 
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(2.79) 



where 



[G + ]r 3 x 



(3x3) 



,[Go 



1 *+ 

J(3x3) ^(3x3) 



(2.80) 



J(3x3) — L^O J(3x3) *(3x3)L 

The color and flavor indices at the propagators in Eq. (|2.79| denote the respective elements of the 
propagators which are given by Eqs. (|2.45|) and (|2.5ip . The kinetic part of the effective action of 
quarks can be simplified by separating the color-flavor part and the Dirac part in the propagators 
which can be achieved by using the energy-chirality projectors, see Sec. 12.1.31 With the relations 
In det A = TrlnA and Trln^a^'P.; = ^ i lnaiTr'P, (see Sees. IB.5 1 and IB.8l in the Appendix), only 
the color-flavor part is remaining in the kinetic part of the effective action of quarks, 



kin 




Gof]' 1 [G^]- 1 

[Gp^rMG+r 1 ' 

rp2 



hi 



In 



rp2 
J<2 



(2.81) 



In this expression, also the sum over c has been performed which gives an extra factor of two. One 
can simplify the effective action of quarks, 
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hi 
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rp2 



j, 2 



In 



(k + S^Y 



{m,<t>i)Y 



j, 2 



(2.82) 



Some computations and definitions have been made in order to get this result. These are explained 
in the following: 
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5<rf 



Mrf) 



M2 



5W 
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M3 



<^2 = ~ (Ms - (K 



1 



A4) 



(2.83) 



are the averaged values and half of the differences of various pairs of quark chemical potentials 
which come from the six 2 x 2-blocks in color-flavor and Nambu-Gorkov space, and 



4>i = 


i gb 
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4> rb 

rsu > 
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•Pds ) 
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fdd i 
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(2.84) 



— Mbb 
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are useful definitions for the gap parameters in order to write Eq. (|2.82|) in a more compact form. 
The gap parameters in the first two lines correspond to the attractive antitriplet channel, while 
the gap parameters in the last line correspond to the repulsive sextet channel. The quasiparticle 
energies which come from the six 2 x 2-blocks in color-flavor and Nambu-Gorkov space are defined 
by the following equation: 



(fi, (j>) = \J{k- e/i) 2 + \(j)\ 2 , 



(2.85) 



cf. Eqs. (jl.l 2|) and (|1.24[) . The dispersion relations (energy eigenvalues) of the quasiquarks which 
come out of the six 2 x 2-blocks in color-flavor and Nambu-Gorkov space read, 



(pi,(/)i) ± 5^ 



(2.86) 



The most complicated expression arises from the determinant of the 3 x 3-matrix in color-flavor 
space in Eq. (|2.81[) which can be calculated analytically. For doing this, one can use computer 
software [130, 131]. The determinant has the following form: 



det (3 x 3) e = kl + bk% + ck% + d 



(2.87) 



The coefficients b, c, and d are rather complicated functions of the quark momentum k, the three 
chemical potentials fi 3 d , and fi h s , and the six gap parameters ipi and ct,, where i — 1,2,3, and 
e = ± which stands for quasiquarks and quasiantiquarks, respectively. The determinant (|2.87p 
can always be factorized as follows: 



det(3x3) e ^n[ fc o-(^r 



(2- 



The functions are the quasiparticle energies, dispersion relations, and energy eigenvalues of the 
three quasiquarks which come from the 6 x 6-block of the inverse Nambu-Gorkov quark propagator. 
In order to get their explicit expressions, one has to solve the cubic equation 



det (3 x 3) e = C 5 + 6£ 2 + c£ + d = , 



(2.89) 



where £ = fcg. By making use of Cardano's formulae, see Sec. IB. 91 in the Appendix, the solutions 
can be presented in an analytical form. Because of their very complicated nature I refrain from 
presenting them explicitly. 

The sum over all fermionic Matsubara frequencies has to be performed in the kinetic term 
of the effective action of quarks (|2.82|) . This is also done with the grand partition function in 
Sec. IB.lOl in the Appendix which is of the same form as the kinetic term of the effective action of 
quarks. Therefore, the result in Sec. lB.lOl can be applied to the kinetic term of the effective action 
of quarks, 
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k i=l e 

2 EEE 
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T 



1 + exp 
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1 + exp 
1 + exp 



4 (Phfa) - &JM 
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T 



(2.90) 



where irrelevant constants are neglected and the vacuum contribution is subtracted. 
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2.1.6 The pressure of color-superconducting quark matter 

Finally, by combining all contributions to the effective action of quarks (|2.54p . I derive the result 

21] 

v 



for the pressure p = tIT* of color-superconducting quark matter, 
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(2.91) 



where I have performed the trace in the potential part of the effective action of quarks (|2.73[) . I 
also converted the sums over all k into integrals by using the relation, 



-Y. 
v 4^ 



+ 00 



d 3 k 

(2^ 
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2^2 



dfc/fc 2 



(2.92) 



The first term in the pressure (|2.91|) is the contribution of leptons, i.e. electrons and muons. They 
are added in order to make color-superconducting quark matter electrically neutral. In principle, 
the contribution of neutrinos should be added as well. In this section, however, their contribution 
is neglected which is a good approximation for neutron stars after deleptonization. The dispersion 
relations of the leptons are given by Ep = (k 2 +m 2 j ) 1 / 2 , where m e = 0.51099906 MeV is the electron 
mass and — 105.658389 MeV is the muon mass [16]. In various applications, a bag constant 
could be added to the pressure if necessary. In order to render the integrals in the expression for 
the pressure finite, I introduced a three-momentum cutoff k. In QCD with dynamical gluons, of 
course, such a cutoff would not be necessary. Here, however, I use a model in which the gluon- 
exchange interaction between quarks is approximated by a point-like four-fermion coupling. Such 
a model is nonrenormalizable. 

The stationary conditions of the pressure (|2.91|) with respect to the nine gap parameters read, 



dp 



= 



d<jj 







(2.93) 



cf. Eq. (|1.7p . In order to find the values of the gap parameters, one can solve for them in Eq. (|2.68|) . 
But a much simpler and equivalent method is to solve the gap equations (|2.93|) for the nine gap 
parameters, 
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, (2.94a) 
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which has to be done numerically. In order to obtain these results, I made use of the relation, 

tanh (|J = 1 - 2n F (x) , (2.95) 

where np is the Fermi-Dirac distribution function which was defined by Eq. (|1.63p . 

Matter in neutron stars has to satisfy the conditions of charge neutrality. The condition for 
electric charge neutrality (| 1 . 54|) reads, 
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where is the number density of electrons and muons, respectively, see Eq. (|1.8ip . The conditions 
for color neutrality (I1.60P arc, 
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= .(2.97b) 



The derivatives of the averaged values of various pairs of quark chemical potentials which come out 
of the 2 x 2-blocks in color-flavor and Nambu-Gorkov space with respect to the chemical potential 
of electric charge \xq and the color chemical potentials /i3 and /ig, respectively, are, 



dm 


1 


dm 


1 


dp3 


1 




~3 ' 




~ 6 ' 


d^Q 


~ 6 ' 


djli __ 


1 


dpi 


1 


dp 3 


= 0, 


dfi3 


~4 ' 


<9/i3 


~ 4 ' 


dfJ>3 


dp,! __ 


1 


dpi 


1 


dp 3 


1 


dfi 8 


4^3 ' 


dfj, s 


4^3 ' 


du 8 


2V3 



(2.98) 



The derivatives of half of the differences of various pairs of quark chemical potentials which come 
out of the 2 x 2-blocks in color-flavor and Nambu-Gorkov space with respect to the chemical 
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potential of electric charge fiQ and the color chemical potentials fj, 3 and fig, respectively, read, 
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The number density of quarks with color i and flavor / can be obtained by 
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(2.99) 
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(2.100) 



By using this relation, the number density of red up, green down, and blue strange quarks respec- 
tively reads, 
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(2.101) 



The number density of green strange, red strange, and red down quarks (upper sign) as well as 
the number density of blue down, blue up, and green up quarks (lower sign) is, 
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(2.102) 



In this equation for the quark number densities, those fij and 8fj,j (|2.83|) have to be taken which 
contain the respective quark chemical potential that corresponds to the respective quark number 
density which one wants to calculate. Since the quasiparticle distribution relations are obtained 
analytically by solving the cubic equation (|2.89j) . the derivatives of ef with respect to the gap 
parameters and the chemical potentials can also be computed analytically. Because of their very 
complicated nature, I refrain from presenting them explicitly. 

All these thermodynamic quantities are valid for nonzero temperature. I also want to show 
the results in the special case of zero temperature. Therefore, Eqs. (|1.64|) and 1(1.68)) . and the 
following relations are useful: 
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and 
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[-e£ (ft, + M , (2.104a) 
6[-fZ(fli,<t>i) + \8iJn\]Bga(5iJn) . (2.104b) 



In this context, it is important to mention that the quasiparticle energies e^[pi,(pi) and ef are 
positive numbers. 

In the limit of zero temperature, the pressure of color-superconducting quark matter ((2.91)) 
reads, 
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where 



lif = /2i ± ^/Sf^iSf • (2-106) 
At zero temperature, the gap equations (|2.94|) are, 
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From Eq. (|2.96[) . the condition for electric charge neutrality at zero temperature can be easily 
derived, 
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where the number density of electrons and muons respectively is given by Eq. p.83p . The color 
charge neutrality conditions (12. 97|) at zero temperature read, 
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The quark number densities of red up, green down, and blue strange quarks at zero temperature 
can be easiliy derived from Eq. (|2.101|) . 
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(2.110) 



The quark number densities of green strange, red strange, and red down quarks (upper sign) as 
well as the number densities of blue down, blue up, and green up quarks (lower sign) at zero 
temperature can be obtained by Eq. (|2.102[) . 

i 1 f f^ 3 ,2 P-j - ek f K 2 pj - ek \ sgn(<5/x,-) r, ,s3 / _n3 
"-f = 7T^> / dkk z -fj- — — T-+ / dkk^-fj- — — ± fa / !T JJ Uf) - U ) 

(2.1H) 

I use the following model parameters: the strength of the diquark coupling and the value of 
the cutoff are fixed as, 

o 2 

^- = 45.1467 GeV~ 2 , k = 0.6533 GeV . (2.112) 
A 2 
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The model parameters are chosen to reproduce several key observables of vacuum QCD such as 
the pion decay constant. 

For given temperature T, strange quark mass m s , and quark chemical potential fi, I solved 
a coupled system of twelve non-linear equations, i.e. nine gap equations and three neutrality 
conditions, in order to obtain the values of the nine gap parameters and of the chemical potentials 
of electric and color charge. This was done numerically [132]. 



2.1.7 Results at zero temperature 

In this subsection, I focus on three-flavor quark matter at zero temperature. It is clear that, for 
small and moderate values of the strange quark mass, the ground state of neutral quark matter 
should correspond to either the regular (gapped) CFL phase or the gCFL phase. At very large 
strange quark mass and/or relatively weak coupling, the ground state can also be either a regular 
(gapped) or gapless 2SC color superconductor. 
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Figure 2.2: The absolute values of the gap parameters as a function of m 2 s /fjL for electrical and color 
neutral color-superconducting quark matter at T — and fi — 500 MeV. The actual values 
of the gap parameters shown in panel a) are negative. 



In order to see how the phase structure of neutral three-flavor quark matter changes with the 
strange quark mass m s , I solve a coupled set of twelve equations, i.e., nine gap equations (|2.107|) 
and three neutrality conditions (|2. 108|) and (|2. 109|) . for various values of m s , keeping the quark 
chemical potential fixed. In this calculation I take /j, = 500 MeV. The results for the absolute 
values of the gap parameters and the chemical potentials [1q, fi^, and fig are shown in Figs. I2.2I 
and !2.3l respectively. Note that, strictly speaking, the gap parameters do not coincide with the 
actual values of the gaps in the quasiparticle spectra. In the case of the CFL phase, for example, 
there is a degenerate octet of quasiparticles with a gap </> oc tct = \fa\ an d a singlet state with 
a gap </> s ingict = 3(^1 - |0i | . In the CFL phase, fa = fa = fa < 0, tpi — ip 2 = ^3 > 0, and 
<7i = (72 = 03 > 0. Also, in the CFL phase, the following relation between the gap parameters 
is satisfied: cr, = ifi — \ fa\ = 20( 6 ,6)) * = 1, 2, 3, where 0(6,6) is the sextet gap in the notation of 
Rcf [19]. 
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Figure 2.3: The electrical and color chemical potentials as a function of m 2 s /[i of electrical and color 
neutral color-superconducting quark matter at T = and ji — 500 MeV. 




Figure 2.4: The absolute values of the gap parameters as a function of \i of electrical and color neutral 
color-superconducting quark matter at T — and m s — 300 MeV. The actual values of the 
gap parameters shown in panel a) are negative. 
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The results in Figs. 12.2 1 and 12 .31 extend the results of Refs. [90,91] by considering a more general 
ansatz for the gap matrix that takes into account, in particular, the pairing in the symmetric sextet 
channel. The effect of including pairing in the symmetric channel is a splitting between the pairs 
of gaps (| 4>\ |, (fx), (| 02 1, (f2), and (|<fo|, ^3) that is also reflected in the change of the quasiparticle 
spectra. In agreement with the general arguments of Refs. [18, 19, 109], the symmetric sextet gaps 
are rather small, see Fig. 12.21 c) . This explains the fact why the splittings between the above 
mentioned pairs of gap parameters are not very large [compare the results in Figs. 12.21 a) and b)]. 




The effects of the strange quark mass are incorporated by a shift of the strange quark chemical 
potential (|2.35[) . Such an approach is certainly reliable at small values of the strange quark mass. 
However, I assumed that it is also qualitatively correct at large values of the strange quark mass. 
By comparing the results for the gap parameters shown in panel a) of Fig. 12. 21 with the results for 
the gap parameters obtained by taking the strange quark mass properly into account shown in the 
left panel of Fig. 1 in Ref. [113], one can see that both results are qualitatively in good agreement. 

I confirm that the phase transition from the CFL phase to the gCFL phase happens at a critical 
value of the parameter m 2 s / n that is in good agreement with the simple estimate of Refs. [90,91], 

2 

in 

— ^ ~ 2A ~ 190 MeV . (2.113) 



The qualitative results for the chemical potentials uq, fi 3: and u 8 in Fig. l2.3l are in agreement with 
the corresponding results obtained in Refs. [90, 91] as well. In my notation, the color chemical 
potential u% fulfills the identity, 



ml 



M = -^~, (2.114) 

in the CFL phase, cf. Eq. (|1.100|) . While the CFL phase requires no electrons to remain neutral, 
the pairing in the gCFL phase is distorted and a nonzero density of electrons appears. This is 
seen directly from the dependence of the chemical potential of electric charge uq in Fig. 12. 31 which 
becomes nonzero only in the gCFL phase. This observation led the authors of Refs. [90, 91] to 
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the conclusion that the phase transition between the CFL and gCFL phase is an insulator-metal 
phase transition, and that the value of the electron density is a convenient order parameter in 
the description of such a transition. In fact, one could also choose one of the differences between 
number densities of mutually paired quarks as an alternative choice for the order parameter [84,85] . 
In either case, there does not seem to exist any continuous symmetry that is associated with such 
an order parameter. To complete the discussion of the chemical potentials, I add that the other 
color chemical potential /X3 is zero only in the CFL phase at T = 0. 

The effects of a nonzero strange quark mass on the phase structure of neutral strange quark 
matter could be viewed from a different standpoint that, in application to stars, may look more 
natural. This is the case where the dependence on the quark chemical potential is studied at a 
fixed value of m s . The corresponding numerical results are shown in Fig. 12.41 (Note once again 
that </>i, where i = 1, 2, 3, have negative values, and I always plotted their absolute values.) In this 
particular calculation I chose m s — 300 MeV. 
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Figure 2.6: The quasiparticle dispersion relations for electrical and color neutral color-superconducting 
quark matter at T = 0, fj, = 500 MeV, and m s = 400 MeV. 



At large values of the quark chemical potential [fi > m^/(2A) ~ 475 MeV which is similar to 
the small strange quark mass limit considered before] , the ground state of quark matter is the CFL 
phase. At smaller values of the chemical potential, the ground state of dense matter is the gCFL 
phase. In this case, there are nine gap parameters all of which are different from each other. One 
could also check that the density of quarks that pair are not equal in the gCFL phase. This can 
be seen from Fig. I2.5I where all nine quark number densities are plotted for the same value of the 
strange quark mass, m s = 300 MeV. Only n d = and n r s ss n h u: all other quark number densities 
are different from each other. This agrees with the general criterion of the appearance of gapless 
phases at T — that was proposed in Refs. [84, 85] in the case of two-flavor quark matter. In the 
ordinary CFL phase, in contrast, one finds that n r u — n d = n s , and n r d — n r s = n£ = n 9 s = n h u = n b d . 

In order to see that the gCFL phase indeed describes a gapless superconductor, it is necessary 
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to show that the dispersion relations of quasiparticles contain gapless excitations. In Fig. 12.61 
the dispersion relations of all nine quasiparticles are plotted. The dispersion relations for the 
corresponding antiparticles are not shown. From Fig. 12.61 one can see that there is indeed a gapless 
mode in the green-strange-blue-down sector. This is the same that was found in Refs. [90,91]. Note 
also that, in agreement with Refs. [90,91], the red-strange-blue-up quasiparticle has a dispersion 
relation that is nearly quadratic, fc £„ — \k — k*\ 2 with k* w 400 MeV for the given choice of 
parameters, see Fig. 12.61 c). The nearly quadratic dispersion relation resembles the situation at 
the transition between the 2SC phase, where n r u = n 9 ,, and the gapless 2SC phase, where n r u and 
n 9 d are different. This explains the approximate equality n r s rs n b u mentioned above. 

2.1.8 Results at nonzero temperature 

In this subsection, I present the results for the phase structure of dense neutral three-flavor quark 
matter in the plane of temperature and m 2 / /j,, as well as in the plane of temperature and quark 
chemical potential. 

I discuss the temperature dependence of the gap parameters in the two qualitatively different 
cases of small and large values of the strange quark mass. As one could see in Sec. I2.1.7[ the 
zero-temperature properties of neutral quark matter were very different in these two limits. 




25 50 75 25 50 75 25 50 75 

T [MeV] T [MeV] T [MeV] 

Figure 2.7: The temperature dependence of the gaps in the case of a small strange quark mass, ml/ fx = 
80 MeV. Note, that the actual values of the gap parameters shown in panel a) are negative. 
The results are plotted for /j, — 500 MeV. 

The results for the temperature dependence of the gap parameters are shown in Figs. 12.71 12.81 
and 12.91 for two different values of the strange quark mass that represent the two qualitatively 
different regimes. In the case of a small strange quark mass (i.e., the case of rn 2 /n = 80 MeV 
which is shown in Figs. 12.71 and 12.8]) . the zero-temperature limit corresponds to the CFL phase. 
This is seen from the fact that the three different gaps shown in every panel of Fig. 12.71 merge as 
T — > 0. At nonzero temperature, on the other hand, the gap parameters are not the same. This 
suggests that, similar to the zero-temperature case of Figs. 12.21 and 12. 4[ a phase transition to the 
gCFL phase happens at some nonzero temperature. However, I shall show below that there is no 
phase transition between the CFL and gCFL phases at any nonzero temperature. Instead, there 
is an insulator-metal crossover transition between the CFL phase and a so-called metallic CFL 
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(mCFL) phase. At this point, all quasiparticles are still gapped. At some higher temperature, the 
mCFL phase is replaced by the gCFL phase. 




67.5 70 72.5 
T [MeV] 



75 




67.5 70 72.5 
T [MeV] 



75 




65 67.5 70 72.5 75 
T [MeV] 



Figure 2.8: 



The near-critical temperature dependence of the gaps in the case of a small strange quark 
mass, rris/fj, = 80 MeV. Note, that the actual values of the gap parameters shown in panel a) 
are negative. The results are plotted for fi = 500 MeV. 
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Figure 2.9: The temperature dependence of the gaps in the case of a large strange quark mass, m 2 s / fx = 
320 MeV. Note, that the actual values of the gap parameters shown in panel a) are negative. 
The results are plotted for /i = 500 MeV. 
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If the temperature is increased even further, there are three consecutive phase transitions. 
These correspond to the three phase transitions predicted in Ref. [Ill] in the limit of a small 
strange quark mass. In order to resolve these, I show a close-up of the near-critical region of 
Fig. 12.71 in Fig. 12.81 The three transitions that one observes are the following: 

1. transition from the gCFL phase to the so-called uSC phase, 

2. transition from the uSC phase to the 2SC phase, 

3. transition from the 2SC phase to the normal quark phase. 

Here, the notation uSC (dSC) stands for superconductivity in which all three colors of the up 
(down) quark flavor participate in diquark pairing [111]. My results differ from those of Ref. [Ill] 
in that the dSC phase is replaced by the uSC phase. The reason is that, in my case, the first 
gaps which vanish with increasing temperature are 4>i and ipi, see Fig. 12.81 while in their case A2 
(corresponding to my (f>2 and 1^2) disappears first. The authors of Ref. [112] also obtain the dSC 
phase instead of the uSC phase at temperatures close to the critical temperature and at small 
values of the strange quark mass. I also obtain the dSC phase with my numerical calculations 
when I use the model parameters of Ref. [112]. In fact, the main difference between these two 
studies is the value of the cutoff parameter in the model. From this, I conclude that the size of 
the uSC or dSC region, respectively, in the phase diagram is particularly sensitive to the choice of 
the cutoff parameter in the model [113]. Although I use the same terms for the phases that were 
introduced in Ref. [Ill], I distinguish between the gapped phases (e.g., CFL and mCFL phase) 
and the gapless phases (e.g., gCFL). Also, in order to reflect the physical properties of the mCFL 
phase, I prefer to use the term metallic CFL, instead of modified CFL as in Ref. [111]. (Note that, 
in that work, the mCFL phase also encompasses the gCFL phase.) 
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Figure 2.10: The temperature dependence of the electrical and color chemical potentials for m'/n = 
80 MeV (left panel) and for m^/fj, — 320 MeV (right panel). The quark chemical potential 
is taken to be ^ = 500 MeV. 



In the case of a large strange quark mass (i.e., the case oim^/fi = 320 MeV shown in Fig. 12. 9p . 
the zero-temperature limit corresponds to the gCFL phase. By looking at the corresponding tem- 
perature dependence of the gap parameters, one can see that this case is a natural generalization 



2.1 The phase diagram of massless quarks 



65 



of the previous limit of a small strange quark mass. There are also three consecutive phase tran- 
sitions. It is noticeable, however, that the separation between the different transitions becomes 
much wider at large m s . 




Figure 2.11: The phase diagram of massless neutral three-flavor quark matter in the plane of temperature 
and m 2 a / fj, (left panel) and in the plane of temperature and quark chemical potential (right 
panel). The results in the left panel are for a fixed value of the quark chemical potential, 
H = 500 MeV. The results in the right panel are for a fixed value of the strange quark mass, 
m 3 — 250 MeV. The dashed lines are associated with the appearance of additional gapless 
modes in the spectra. The dotted lines indicate the insulator-metal crossover. 



I now take a closer look at the transition between the CFL, the mCFL, and the gCFL phase. At 
zero temperature, there was no symmetry connected with the order parameter, i.e., the number 
density of electrons, that is associated with the CFL — > gCFL phase transition. At nonzero 
temperature, the electron density is not strictly zero in the CFL phase as soon as m s =/= 0. Indeed, 
the arguments of Ref. [20] regarding the enforced neutrality of the CFL phase do not apply at 
T ^ 0. This is why I conclude that the insulator- metal transition between the CFL and the 
mCFL phase is just a smooth crossover at T ^ 0. Of course, in principle, one can never exclude 
the existence of a first-order phase transition. My numerical analysis, however, reveals a crossover. 
The transition can only be identified by a rapid increase of the chemical potential of electric charge 
or the electron density in a relatively narrow window of temperatures, see panel a) in Fig. 12.101 
The location of the maximum of the corresponding susceptibility (i.e., driQ/dT) is then associated 
with the transition point. 

The transition between the mCFL and the gCFL phase corresponds to the appearance of 
gapless quasiparticle modes in the spectrum. There is no way of telling from the temperature 
dependences in Figs. r2~71 12.81 and 12.91 whether the corresponding CFL and/or 2SC phases are 
gapless or not. This additional piece of information can only be extracted from the behavior of 
the quasiparticle spectra. I also investigated them, but I do not show them explicitly. 

My results for the phase structure of dense neutral three-flavor quark matter are summarized 
in Fig. 12.111 I show the phase diagram in the plane of temperature and m 2 s /fj, at a fixed value of 
the quark chemical potential, /i = 500 MeV, and in the plane of temperature and quark chemical 
potential at a fixed value of the strange quark mass, m s = 250 MeV. The three solid lines denote 
the three phase transitions discussed above. The two dashed lines mark the appearance of gapless 
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modes in the mCFL and 2SC phases. I term these as the mCFL — > gCFL and 2SC — > g'2SC 
crossover transitions. In addition, as I mentioned above, there is also an insulator-metal type 
transition between the CFL and mCFL phase. This is marked by the dotted lines in Fig. 12.111 

The effects of the strange quark mass are incorporated only by a shift of the strange quark 
chemical potential (|2.35p . By comparing the phase diagram shown in panel a) of Fig. [2~TT1 with that 
one shown in the right panel of Fig. 1 in Ref. [113], where the strange quark mass is properly taken 
into account, one can see that both phase diagrams are quantitatively in good agreement for low 
strange quark masses. For large strange quark masses, the phase diagrams are still qualitatively 
in good agreement. 

2.2 The phase diagram with a self-consistent treatment of 
quark masses 

The first attempt to obtain the phase diagram of dense, locally neutral three-flavor quark matter 
as a function of the strange quark mass, the quark chemical potential, and the temperature was 
made in Ref. [110]. I presented this work in Sec. 12.11 It was shown that, at zero temperature and 
small values of the strange quark mass, the ground state of quark matter corresponds to the CFL 
phase [18, 19, 109]. At some critical value of the strange quark mass, this is replaced by the gCFL 
phase [90,91]. In addition, several other phases were found at nonzero temperature. For instance, 
it was shown that there should exist a metallic CFL (mCFL) phase, a so-called uSC phase, as well 
as the standard 2SC phase [15,32,107,108] and the g2SC phase [84,85]. 

In Ref. [110], the effect of the strange quark mass was incorporated only approximately by a 
shift of the strange quark chemical potential, /i s — ► /i s — m^/(2/i). Such an approach is reliable 
at small values of the strange quark mass but for larger values of the strange quark mass, the 
results are only qualitatively in good agreement by comparing them to the case, where the strange 
quark mass is properly taken into account [113]. The phase diagram of Ref. [110] was further 
developed in Refs. [112, 113] where the shift-approximation in dealing with the strange quark was 
not employed any more. So far, however, quark masses were treated as free parameters, rather 
than dynamically generated quantities. 

In this section, I study the phase diagram of dense, locally neutral three-flavor quark matter 
within the framework of a NJL model [14, 103-106], treating dynamically generated quark masses 
self-consistently. I introduce the model and, within this model, I derive a complete set of gap 
equations and charge neutrality conditions. The results are summarized in the phase diagram in 
the plane of temperature T and quark chemical potential /i. Some results within this approach 
were also obtained in Refs. [115, 117, 118, 133]. 

2.2.1 Model 

In this section, I use a three-flavor NJL model [103-106]. The Lagrangian density is given by [106, 



114,116], 



C = ip {ift— to) ip 



8 




- K {det, [ip (1 + 75) V] + det/ [ip (1 - 75) 1>] } 



(2.115) 



where the quark spinor field, 



(2.116) 
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carries color (a = r, g, b) and flavor (a = u, d, s) indices. The matrix of quark current masses is 
given by 

rh = diag(m u ,m d ,m s ,m u ,m d ,m s ,rn u ,m d ,m s ) , (2.117) 

and A a with a — 1, . . . , 8 are the Gell-Mann matrices in flavor space (|A.27|) . Ao = -\/2/3. In the 
Lagrangian density, I use the charge- conjugate quark spinors (|A.14j) . The Dirac conjugate quark 
spinor is given by "ip = V^7o- 

Throughout this work, I shall assume isospin symmetry on the Lagrangian level, vn u d = m u — 
m d , whereas m s will be different, thus explicitly breaking [SU(3)f] flavor symmetry. The terms 
in the Lagrangian density (|2.115p can straightforwardly be generalized to any number of flavors 
Nf. I only deal with the lighter up, down, and strange quarks because charm, bottom, and top 
quarks are too heavy to occur in the cores of neutron stars. 

The term in the first line in the Lagrangian density (|2.115p . ^Dirac, is the Dirac Lagrangian 
density of non-interacting massive quarks. The terms in the second line correspond to a U(Nf)t x 
U(Nf) r symmetric four-point interaction with the scalar coupling constant Gs- Therefore, this 
part of the Lagrangian density is termed symmetric Lagrangian density, £ S ym- The term in the 
third line of the Lagrangian density, £diquark, describes a scalar diquark interaction in the color- 
antitriplet and flavor-antitriplet channel. For symmetry reasons there also should be a pseudoscalar 
diquark interaction with the same diquark coupling constant Gd- This term would be important 
to describe Goldstone boson condensation in the CFL phase [134]. In this study, however, I neglect 
this possibility, and therefore drop the pseudoscalar diquark term. The terms in the last line of 
the Lagrangian density, £dct, are called the 't Hooft interaction, where K is the coupling constant. 
The determinants in these terms have to be taken in flavor space which means that £det is a 
maximally flavor-mixing 2V/-point interaction, involving an incoming and an outgoing quark of 
each flavor. Consequently, for three flavors, the 't Hooft term is a six-point interaction term of 
the form, 

det, $Cty) = J2 e ^ k ^O^i) (faCtyj) {i) s Oi> k ) , (2.118) 

where i, j, and k are flavor indices. 

The 't Hooft interaction is SU(Nf)e X SU(Nf) r symmetric, but it breaks the J7a(1) axial 
symmetry [135], which was left unbroken by £ sym . It thus translates the Ua(1) anomaly, which 
in QCD arises at quantum level from the gluon sector, to a tree-level interaction in a pure quark 
model. The 't Hooft term is phenomenologically important to get the correct mass splitting of 
the rj and rj mesons. In the chiral limit (m u — m d — m s — 0), the rf mass is lifted to a nonzero 
value by Cdet, while the other pseudoscalar mesons, including the ry, remain massless. There are 
many other terms which are consistent with the symmetries and which could be added to the 
Lagrangian. But for simplicity, however, I shall restrict myself to the Lagrangian density (|2.115p . 
This Lagrangian density should be viewed as an effective model of strongly interacting matter 
that captures at least some key features of QCD dynamics. The Lagrangian density contains 
three interaction terms which are chosen to respect the symmetries of QCD. 

In the Lagrangian density (|2.115|1 . there are six parameters: the quark current masses m u ,d 
and m s , the coupling constants Gs, Gd, and K, and the cutoff A. I use the following set of model 
parameters [136]: 

m Ui d — 5.5 MeV , 

m s = 140.7 MeV , 
G S A 2 = 1.835 , (2.119) 
KA 5 = 12.36 , 
A = 602.3 MeV . 

After fixing the current masses of the up and down quarks at equal values, m Ul d = 5.5 MeV, the 
other four parameters in Eq. (|2.119p are chosen to reproduce the following four observables of 
vacuum QCD [136]: = 135.0 MeV, m K = 497.7 MeV, m n > = 957.8 MeV, and U = 92.4 MeV. 
This parameter set gives m v = 514.8 MeV [136]. In Ref. [136], the diquark coupling constant Gd 
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was not fixed by the fit of the meson spectrum in vacuum. In general, it is expected to be of the 
same order as the quark-antiquark coupling constant Gs- In this section, I study the following 
two cases in detail: 

1. the regime of intermediate diquark coupling strength with Gd — § Gs, 

2. the regime of strong diquark coupling strength with Gd = Gs- 

In fact, the value of the diquark coupling constant for the regime of intermediate coupling strength 
is used in many studies because this value for the diquark coupling constant is predicted by the 
Fierz transformation in vacuum. Therefore, it is a standard value, and, if not mentioned otherwise, 
throughout this thesis, the standard value of the diquark coupling constant, 



Gd = 



(2.120) 



is used. 

In order to obtain the Lagrangian density (|2.115|) in mean-field approximation, one can adopt 
the rules [105], 



^iipiipjipjipkipk 



^ {i>iipi)('ipjipj)'ipkipk ~ 2 {$i$i)$ji 3 j)('<Pk'ipk) 



(2.121a) 
(2.121b) 



(cyclic) 



where the last terms on the right-hand side are indispensable to avoid a double counting of the 
interaction energy. This is why there is also a factor of two in the six-fermion interaction. One 
should note that in this approach, the condensates are treated as variational parameters which 
are determined so as to give a maximum point of the pressure. 

The diquark condensates A c are analogously defined as in Eq. p.5[) . The quark-antiquark 
condensates are defined as, 

a a cx $M . (2.122) 
Herewith, I arrive at the following Langangian density in mean-field approximation: 



c = $ M)tp - 2G S J2 °l + \¥r$ c + \i>c* + il> - 7J- £ I Ac < 

a— 1 c— 1 

where the constituent quark mass matrix is defined as, 

M = diag {M u , M d , M s , M u ,M d , M s , M u ,M d , M s ) 

The constituent mass of each quark flavor can be obtained by 

M a = m a — AGs<J a + 2Ka^a 1 , 

where the set of flavor indices (a, (3, 7) is a permutation of (u,d,s). The Lagrangian density in 
mean-field approximation contains the gap matrices in color-flavor and Dirac space, 

$ ± = ±75$ , (2.126) 

which fulfill the relation <I>~ = 7o( ( I )+ )^7o- The color-flavor part of the gap matrices reads, 



4Ka u a d a s , (2.123) 

(2.124) 
(2.125) 
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where A c are real-valued gap parameters. Here, as before, a and b refer to the color components, 
and a and (3 refer to the flavor components. The gap parameters Ai, A 2 , and A 3 correspond to 
the down-strange, the up-strange, and the up-down diquark condensates, respectively. All three of 
them originate from the color- antitriplet, flavor- antitriplet diquark pairing channel. For simplicity, 
the color and flavor symmetric condensates are neglected in this study. They were shown to be 
small [see Figs. 12.21 c) . 12. 41 c). 12.71 c) . andHHJc)] and not crucial for the qualitative understanding 
of the phase diagram [110]. 

I should note that I have restricted myself to field contractions corresponding to the Hartree 
approximation in Eq. (|2.123[) . In a more complete treatment, among others, the 't Hooft inter- 
action term gives also rise to mixed contributions containing both diquark and quark-antiquark 
condensates, i.e., cx J2a=i a a\A a \ 2 [137]. In this study, as in Refs. [106,133], I neglected such 
terms for simplicity. While their presence may change the results quantitatively, one does not 
expect them to modify the qualitative structure of the phase diagram. 

Up to irrelevant constants, the grand partition function is given by 

Z = J V$Vipexp{l [V>,V]} , (2.128) 

where 

I$,ll>]= [ (C + (jj\f + MqA/q + HaNa) , (2.129) 



A' 

is the action. The conserved quantities, 

AA=^7o^, A/q = ^ 7 oW, N a = ^M, (2.130) 

are the quark number density operator, the operator of electric charge density of the quarks, and 
the operators of color charge densities of the quarks, respectively. The action can be split into a 
kinetic and a potential part, 

l[-0,V]=Ikm[<M]+/pot. (2-131) 

By inserting the Lagrangian density in mean- field approximation Q2.123j) into the action (|2.129j) . 
one obtains for the kinetic part of the action, 

/kin^V]/ (^Po + rV + ^#^c + ^"c^) , (2.132) 
and for the potential part of the action, 



Ipot — 



j x (yZGs^al- ^Y^\& c \ 2 +AKa u a d a)j . (2.133) 



The massive inverse Dirac propagator for quarks and charge-conjugate quarks, respectively, is 
given by 

[Gj]- 1 =iP±p.i Q -M, (2.134) 

where 

fi = diag (/^, fi d , fi r s , n B u , fi 9 d , n g s , n h u , fi b d , nl) , (2.135) 

is the matrix of quark chemical potentials. The chemical potentials for each quark color and flavor 
are defined by Eq. (|1.52|) because quark matter inside neutron stars is in (3 equilibrium. 
By introducing the Nambu-Gorkov basis, 

*=(/ c )> ( 2 - 136 ) 

the kinetic part of the action can be rewritten as 

/ki„ [*,*] = ~ f ^S-H , (2.137) 
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where I have used the inverse quark propagator in Nambu-Gorkov space, 



In momentum space, the massive inverse Dirac propagator for quarks and charge-conjugate 
quarks, respectively, reads, 

[Gjr 1 = 7o(fco±A) -7-k-M. (2.139) 

The kinetic part of the grand partition function (|2.128j) is Fourier transformed in Sec. IB. 51 in the 
Appendix. With that result, one obtains for the pressure of color-superconducting quark matter, 
p = Sln2, 



w ]T lndct (^) -ZGsYsrt-^-Y, l A -| 2 + + Pp , (2-140) 

K ^ ^ ol—1 ® c—1 



P = 



where I also added the contribution of leptons, pp, which will be specified later. 

The most complicated expression in the pressure is the determinant because the inverse quark 
propagator (|2.138[) is a 72 x 72-matrix with the following substructure: it is a 2 x 2-matrix in 
Nambu-Gorkov space, a 9 x 9-matrix in color-flavor space and a 4 x 4-matrix in Dirac space. By 
using the relation det (cA) — c d det A [129], where c is a factor and A is a quadratic matrix with 
dimension d, one obtains for the kinetic part of the pressure, 



K v 7 K 



where it is useful to rewrite the determinant of the inverse quark propagator with the relation 
det (AB) = det A ■ det B, where the matrices A and B are quadratic [129], 

detS"" 1 = det (7070S'" 1 ) = det 70 det (-/qS" 1 ) = det (70 S -1 ) . (2.142) 

By using the spin projectors (j A.4|) . 

P, = i(l + a<T.k), (2.143) 

where s = ± stands for projections onto states with spin up or spin down, respectively, one can 
rewrite the matrices 7o[Gq ] _1 and 7o$ ± as [114], 

loicjr 1 = E ( "* *' to ±7+ m ) p * • w 

^ = ±^(.4 o) P " (2 - 145) 

The spin projectors are 2 x 2-matrices in spin space. Thus, I have separated the spin space from 
the color-flavor space and particle-antiparticle space. The subdivision of the Dirac space into the 
spin space and the particle-antiparticle space is necessary in order to take the different quark 
masses properly into account. 

By making use of the definition in Eq. (|2.138|> . as well as Eqs. (12.1441) and (|2.145[) . I obtain 
the following representation: 

^S~ X =Y J S- X V S , (2.146) 

S 

where 

§- 1 = k -M a , (2.147) 
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and 



M s = 



f-jl + M sk $ 
sk -ft - M -$ 
fi + M sk 

\ 4 sk fi-Mj 



(2.148) 



(with s = ±) is real and symmetric. By using the spin projectors, the spin space is separated 
from the Nambu-Gorkov space, color-flavor space, and particle-antiparticle space. Therefore, 
the matrices S~ 1 and A4 S are 36 x 36-matrices and have the following substructure: they are 
2 x 2-matrices in Nambu-Gorkov space, 9 x 9-matrices in color-flavor space, and 2 x 2-matrices in 
particle-antiparticle space. Each clement of the matrices which are shown in Eqs. (|2.144j) . (|2.145|) . 
and (I2.148|) are 9 x 9-matrices in color-flavor space. 

Since there is no explicit energy dependence in M Sl its eigenvalues e; determine the quasipar- 
ticle dispersion relations, kg = €i(k). By using the properties of the projectors V s as well as the 
matrix relations, In det A = TrhiA, and Trln^a^'P; = ^jlnaiTr'Pi, see Sees. IB.61 and IB~8l in 
the Appendix, I derive, 



lndet( 7o S'- 1 ) = Trln( 7o 5'- 1 ) ^Tcbx^S^V. 

8 

= Tr m S- 1 = ^ In det S' 1 = In (det SZ 1 ■ det SZ 1 



(2.149) 



where the traces in the first line run over Nambu-Gorkov, color-flavor, particle-antiparticle, and 
spin indices, while the trace in the second line only runs over Nambu-Gorkov, color-flavor, and 
particle-antiparticle indices. It turns out that the two determinants appearing on the right-hand 
side of this equation are equal, i.e., det SZ 1 = det SZ 1 ■ From the physics viewpoint, this identity 
reflects the twofold spin degeneracy of the spectrum of quark quasiparticles. The formal proof of 
this degeneracy is straightforward after noticing that the following matrix relation is satisfied: 



where 



szl = nsz 1 ^ 1 



ft = diag (1,-1, -1,1) 



(2.150) 



(2.151) 



is a matrix with the properties 1Z = 1Z T = TZ) = TlT 1 and det 1Z = 1. 

Another observation which turns out to be helpful in the calculation is that the determinant 
detS'j 1 (/co) is an even function of fco, i.e., det S'j' 1 (— kg) = det 5"~ 1 (/co) ■ This is a formal conse- 
quence of the following matrix relation: 



S- 1 (-ko) = -BS- 1 (k )B- 



where 



B = 



/ i\ 

-i 

i 

\-i 0/ 



(2.152) 



(2.153) 



is a matrix with the properties B = B** = B^ 1 and det 2? = 1. The invariance of the determinant 
detS'7 1 (A;o) with respect to the change of the energy sign, k — > — k , is directly related to the 
use of the Nambu-Gorkov basis for quark fields. In this basis, for each quasiparticle excitation 
with a positive energy ko = e(k), there exists a corresponding excitation with a negative energy 
ko = —e(k). Therefore, the result for the determinant should read, 



18 



dcts-^n^o-^ 



(2.154) 
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The expression k$ — e f i n Eq. (|2.154p has to be squared because of the spin degeneracy. Because of 
the artificial Nambu-Gorkov degeneracy, fcg — ef = (k — ej) (k + e*), where u are the 18 positive 
energy eigenvalues of Eq. (|2.148|) . 

In order to simplify the numerical calculation of the eigenvalues of the matrix defined in 
Eq. (|2.148p . I first write it in a block-diagonal form. With a proper ordering of its 36 rows and 
36 columns, it decomposes into six diagonal blocks of dimension 4x4 and one diagonal block of 
dimension 12 x 12. The explicit form of these blocks reads, 

f-p A + M d k -A 3 \ 

k -fjT d -M d A 3 
A 3 nl + Mv k 

-A 3 k ni -M,, J 

-A 3 \ 


k 

-Hi + M u j 
k -A 2 \ 

-p B - Ms A 2 
A 2 p u + M u 
k 



M 



M 



(2) 



(3) 



V 

(tfi-Mi k 

k Hi + M d 
A 3 






A 3 



ifi 



H b u 



M { f> = 





A 2 

-Hu - M u 



M 



(6) 




M u ) 



\ 



A 2 


k 

+ M U / 

-Aj 


k 

-M d J 

Ai 


k 



P d + M d J 



(2.155a) 



(2.155b) 



(2.155c) 



(2.155d) 



(2.155e) 



(2.155f) 



(2.155g) 



(7) 

Here, I have only shown the color-flavor structure of the matrix A4\ because the repesentation 
of the total structure of the 12 x 12-matrix .M? would be too large. Each color-flavor element in 
M.\ is a 4 x 4-matrix in Nambu-Gorkov and particle-antiparticle space, 



( 



Ml = 



Hu ~ Mu 
k 
Q 




M d 



k 



-H b s ~ M s 
k 





k 

-Hi + Mu 



k 






Hi - M u 
k 



Hv 



M u j 



-Pi 



M d 






k 

-H b s + M s 





Pi 







-M d k 
k P d + M d J 

\ 



- M, k 



Ms J 



Ai = 



A 2 = 



-Ai\ 
Ai 
Ai 



V-Ai 





A 2 

-A 2 





A 3 

-A 3 



0/ 



-A 2 
A 2 



-A 3 \ 
A 3 





(2.156) 
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The inverse quark propagator is a 72 x 72-matrix, and therefore, it has 72 energy eigenvalues. Be- 
cause of spin degeneracy half of them are equal. Out of 36 eigenvalues from all seven blocks (|2.155|) , 
there are 18 positive and 18 negative eigenvalues because of the Nambu-Gorkov degeneracy. Out 
of total 18 positive eigenvalues, nine of them correspond to the nine quasiquarks and the other 
nine correspond to the nine quasiantiquarks. I obtained the 36 eigenvalues of Eq. (|2.155[) numeri- 
cally [132]. 

Here, it might be interesting to note that the eigenvalues of the first six 4x4 matrices in 
Eq. l|2.155j) can be calculated analytically in the limit when two quark masses appearing in each 
of them are equal. For example, when M c i = M u . the four eigenvalues of matrix M.+ are given 
by 

A« = ±J^±KZ/M^T¥^Al - , (2 .i57) 

(2) 

while the eigenvalues of M\_ differ only by the sign in front of the second term, 

Af > = ±^(^±^±VWTv) 2 + Aj + . (3U58) 

When the value of 6M = Md — M u is nonzero but small, the corrections to the above eigenvalues 
are ±M U 8M / '(2-v/ 'M 2 + k 2 ) with the plus sign in the case of antiparticle modes, and the minus 
sign in the case of particle modes. The eigenvalues of M.^? and in the limit M s — M u , as 

well as the eigenvalues of A4^ and in the limit M s = Md, are similar. 

It turns out that in regular color-superconducting quark phases, there are two positive-energy 
eigenvalues and two negative-energy eigenvalues in each of the first six matrices in Eq. (|2.155|) . 
The two negative-energy eigenvalues are a consequence of the artificial Nambu-Gorkov degeneracy. 
These negative-energy eigenvalues do not have the same absolute value as the positive-energy 
eigenvalues in the same matrix. That it why these negative-energy eigenvalues have their positive 
counterparts in another of these first six matrices. It is easy to extract the 18 positive-energy 
eigenvalues out of the total 36 energy eigenvalues. This is done by sorting the energy eigenvalues 
in descending order [132] in each of the first six matrices. After that, I take the first two energy 
eigenvalues from each of the first six matrices. The lower one of these positive energy eigenvalues 
corresponds to a quasiquark and the larger one corresponds to a quasiantiquark. It is necessary 
to remember to which of the first six matrices they belong to in order to determine from which 
of the matrices a gapless mode comes from. In gapless color-superconducting quark phases, an 
energy eigenvalue of at least one of the first six matrices becomes positive (negative) which has 
been negative (positive) in regular color-superconducting quark phases. This means in the case 
of one gapless mode that one of the first six matrices has three positive but only one negative 
energy eigenvalues, and another one has only one positive but three negative energy eigenvalues. 
The matrix from which arise three positive-energy eigenvalues shows a gapless mode. Therefore, 
if one sorts the energy eigenvalues in each of the first six matrices in descending order and then 
extracts from each of these matrices the first two energy eigenvalues, one obtains at least one 
negative-energy eigenvalue in a gapless color-superconducting quark phase by this method. Since 
I do not need such negative energy eigenvalues, I take the absolute value of them. Because of the 
Nambu-Gorkov degeneracy, they can be found by using the following relation: 

ei (M%* = -e, (m™) , (2.159) 

where n — 1,2,3. This relation means that for example the matrices M.+ and have the 

same four absolute values of energy eigenvalues but in each of these two matrices the four energy 
eigenvalues have different signs. So, if one finds a negative-energy eigenvalue in matrix which 

has been positive in a regular phase before, then the gapless mode arises from M + . 
The relation (|2.159|) is easily proven by introducing the matrix, 



« = diag (1,-1, 1,-1) 



(2.160) 



74 



The phase diagram of neutral quark matter 



with the properties U = U T =W = U 1 and detU = 1, 



- fco - M ( + l) = -U (fc - M% n - X) ) U- 1 , (2.161) 

where n = 1, 2, 3. 

The 12 x 12-matrix in Eq. (|2.155p . M.+ , always has six positive and six negative eigenvalues. 

(7) 

Therefore, it is easy to extract the six positive-energy eigenvalues out of M + . I sort the energy 
eigenvalues in descending order [132] and then take the first six one. The first three out of the six 
positive-energy eigenvalues correspond to the three quasiantiquarks while the second three out of 
the six positive-energy eigenvalues correspond to three quasiparticles. The remaining six negative- 
energy eigenvalues arise because of the artificial Nambu-Gorkov degeneracy. Their absolute values 

(7) 

are equal to the positive-energy eigenvalues of M + . 

By using Eq. ((2.154[) . the kinetic part of the pressure (|2.141|) can be rewritten as 



Pku 



EE ln 



i=l K 



k 2 



rp2 



(2.162) 



After performing the sum over all fcrmionic Matsubara frequencies, see Sec. lB.lOl in the Appendix, 
and by noting that InZkin = yPkm and Sfi = 0, one obtains, 



Pku 
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V 



18 



2Tln 



i=l k 



1 + exp i 



(-?)]} ■ 



(2.163) 



which can be converted into an integral representation by using Eq. (|2.92p . Therefore, the pressure 
of color-superconducting quark matter reads, 
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AK(j u o d a s 



(2.164) 



where the term in the last line is the contribution of leptons, pp, i.e. electrons and unions. They 
are added in order to make color-superconducting quark matter electrically neutral. In principle, 
the contribution of neutrinos should be added as well. In this section, however, their contribution 
is neglected which is a good approximation for neutron stars after deleptonization. The effect 
of neutrino trapping is included in Ref. [116] and is the subject in Sec. 12.31 of my thesis. The 
dispersion relations of the leptons are given by Ep = (k 2 +m|) 1 / 2 , where m e = 0.51099906 MeV is 
the electron mass and to m = 105.658389 MeV is the muon mass [16]. In various applications, a bag 
constant could be added to the pressure, and also the vacuum contribution could be substracted 
from the pressure if necessary. In order to render the integrals in the expression for the pressure 
finite, I used the three-momentum cutoff A. 

In the limit of zero temperature, the pressure of color-superconducting quark matter reads, 



1 



P = 



18 

2r^E 

i=l 



dk k 2 €i 



1 



3 



4G Z 



4Ka u a d a s 



1 



F "dk^. 



a—1 c—1 

(2.165) 

In order to obtain the values for the quark-antiquark condensates a a and the color-superconducting 
gap parameters A c , I solve the following six stationary conditions: 

dp dp 



da a 



= 



= 



(2.166) 



After the quark-antiquark condensates are known, the quark constituent masses can be determined 
by Eq. (|2.125[) . In order to enforce the conditions of local electric and color charge neutrality in 
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color-superconducting quark matter, it is necessary to require the three Eqs. (|1.54p and p.60|) to 
be satisfied. These fix the values of the three corresponding chemical potentials, /iq, ^ and fig,. 
After these are fixed, only the quark chemical potential /i and the temperature T are left as free 
parameters. 

The number density of each quark is given by 



n) - j| , (2.167) 



and the number density of each quark flavor reads, 

b 

*/ = !>/■ (2-168) 

i—r 

Since the energy eigenvalues are computed numerically, the derivatives, i.e. the six gap equations, 
the three neutrality conditions, and the quark number densities have to be computed numerically. 
The number densities of the leptons can be calculated using Eq. (|1.81[) in the case of nonzero 
temperature and Eq. (| 1 .83[) in the case of zero temperature. 



2.2.2 Results 

In order to obtain the phase diagram, I have to find the ground state of quark matter for each 
given set of the parameters in the model. In the case of locally neutral quark matter, there are 
two parameters that should be specified: temperature T and quark chemical potential [i. After 
these are fixed, one has to compare the values of the pressure in all competing neutral phases of 
quark matter. The ground state corresponds to the phase with the highest pressure. 

Before calculating the pressure, given by Eq. (|2.164|) for nonzero temperature and by Eq. (|2.165|) 
for zero temperature, one has to find the values of the chiral and the color-superconducting order 
parameters, a a and A c , as well as the values of the three charge chemical potentials, fiQ, /13 and 
/ig. These are obtained by solving the coupled set of six gap equations (|2.166[) together with the 
three neutrality conditions (|1.54[) and (|1.60p . By using standard numerical recipes [132], it is not 
extremely difficult to find a solution to the given set of nine nonlinear equations. Complications 
arise, however, due to the fact that often the solution is not unique. 

The existence of different solutions to the same set of equations, (|1.54p . (| 1 .60(1 . and ()2.166|) . 
reflects the physical fact that there could exist several competing neutral phases with different 
physical properties. Among these phases, all but one are unstable or metastable. In order to 
take this into account in my study, I look for the solutions of the eight types which are listed in 
Table O 



Phase 


Ax 


A 2 


A 3 


NQ 








2SC 






/ 


2SCus 




/ 




2SCds 


/ 






uSC 




/ 


/ 


dSC 


/ 




/ 


sSC 


/ 


/ 




CFL 


/ 


/ 


/ 



Table 2.1: The classification of eight color-superconducting quark phases. The unmarked gap parameters 
(— ) are zero while the checkmarked gap parameters (/) are nonzero in the respective color- 
superconducting quark phases. Note, that I do not distinguish here between the CFL and the 
mCFL phase [110]. 
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Figure 2.12: The phase diagram of neutral quark matter in the regime of intermediate diquark cou- 
pling strength, Gd = jGs- First-order phase boundaries are indicated by bold solid 
lines, whereas the thin solid lines mark second-order phase boundaries between two phases 
which differ by one or more nonzero diquark condensates. The dashed lines indicate the 
(dis-)appearance of gapless modes in different phases, and they do not correspond to phase 
transitions. 




Figure 2.13: The phase diagram of neutral quark matter in the regime of strong diquark coupling, 
Gd = Gs- The meaning of the various line types is the same as in Fig. 12.121 
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Figure 2.14: The dependence of the quark masses, of the gap parameters, and of the electric and color 
charge chemical potentials on the quark chemical potential at a fixed temperature, T — 
MeV (three upper panels), T = 20 MeV (three middle panels), and T = 40 MeV (three 
lower panels). The diquark coupling strength is Gd = f Gs- 



Then, I calculate the values of the pressure in all nonequivalent phases, and determine the 
ground state as the phase with the highest pressure. After this is done, I additionally study the 
spectrum of low-energy quasiparticles in search for the existence of gapless modes. This allows me 
to refine the specific nature of the ground state. 

In the definition of the eight phases in terms of A c in Table 12. 1[ I have ignored the quark- 
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antiquark condensates a a . In fact, in the chiral limit (m a = 0), the quantities a a are good order 
parameters and I could define additional sub-phases characterized by nonvanishing values of one 
or more a a . With the model parameters at hand, however, chiral symmetry is broken explicitly 
by the nonzero current quark masses, and the values of a a never vanish. Hence, in a strict sense 
it is impossible to define any new phases in terms of a a . 

Of course, this does not exclude the possibility of discontinuous changes in a a at some line in 
the plane of temperature and quark chemical potential, thereby constituting a first-order phase 
transition line. It is generally expected that the chiral phase transition remains first order at low 
temperatures, even for nonzero quark masses. Above some critical temperature, however, this line 
could end in a critical endpoint, and there is only a smooth crossover at higher temperatures. 
Among others, this picture emerges from NJL-model studies, both, without [138] and with [139] 
diquark pairing (see also Ref. [106]). Therefore, I expect a similar behavior in my analysis. 

My numerical results for neutral quark matter are summarized in Figs. [2TT21 and [2. 131 These 
are the phase diagrams in the plane of temperature T and quark chemical potential fj, in the case 
of an intermediate diquark coupling strength, Gd = jGs, and in the case of a strong coupling, 
Gd — Gs, respectively. The corresponding dynamical quark masses, gap parameters, and three 
charge chemical potentials are displayed in Figs. 12.141 and 12.151 respectively. All quantities are 
plotted as functions of /i for three different fixed values of the temperature: T — 0,20,40 MeV 
in the case of G D = \G S (see Fig. [2~T4ll and T = 0,40,60 MeV in the case of G D = G s (see 
Fig. MB- 

Let me begin with the results in the case of the diquark coupling being Gd = jG$- In the 
region of small quark chemical potentials and low temperatures, the phase diagram is dominated 
by the normal phase in which the approximate chiral symmetry is broken, and in which quarks 
have relatively large constituent masses. This is denoted by %SB in Fig. 12.121 With increasing 
the temperature, this phase changes smoothly into a normal quark phase (NQ) in which quark 
masses are relatively small. Because of explicit breaking of the chiral symmetry in the model at 
hand, there is no need for a phase transition between the two regimes. 

However, as pointed out above, the symmetry argument does not exclude the possibility of 
a first-order chiral phase transition. As expected, at lower temperatures I find a line of first- 
order chiral phase transitions. It is located within a relatively narrow window of the quark 
chemical potentials (336 MeV < /x < 368 MeV) which are of the order of the vacuum values 
of the light-quark constituent masses. (For the parameters used in my calculations one obtains 
M u = M d = 367.7 MeV and M s = 549.5 MeV in vacuum [136].) At this critical line, the quark 
chiral condensates, as well as the quark constituent masses, change discontinuously. With increas- 
ing temperature, the size of the discontinuity decreases, and the line terminates at the endpoint 
located at (T cr ,^ cr ) w (56,336) MeV, see Fig. EH 

The location of the critical endpoint is consistent with other mean-field studies of NJL models 
with similar sets of parameters [106, 138, 139]. This agreement does not need to be exact because, 
in contrast to the studies in Refs. [106, 138, 139], here I imposed the condition of electric charge 
neutrality in quark matter. (Note that the color neutrality is satisfied automatically in the normal 
quark phase.) One may argue, however, that the additional constraint of neutrality is unlikely to 
play a big role in the vicinity of the endpoint. 

It is appropriate to mention here that the location of the critical endpoint might be affected 
very much by fluctuations of the composite chiral fields. These are not included in the mean- field 
studies of the NJL model. In fact, this is probably the main reason for their inability to pin 
down the location of the critical endpoint consistent with lattice QCD calculations [5]. It is fair 
to mention that the current lattice QCD calculations are not very reliable at nonzero \x either. 
Therefore, the predictions of this study, as well as of those in Refs. [106, 138, 139], regarding the 
critical endpoint cannot be considered as very reliable. 

When the quark chemical potential exceeds some critical value and the temperature is not too 
large, a Cooper instability with respect to diquark condensation should develop in the system. 
Without enforcing neutrality, i.e., if the chemical potentials of up and down quarks are equal, this 
happens immediately after the chiral phase transition, when the density becomes nonzero [139]. 
In the present model, this is not the case at low temperatures. 
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Figure 2.15: The dependence of the quark masses, of the gap parameters, and of the electric and color 
charge chemical potentials on the quark chemical potential at a fixed temperature, T — 
MeV (three upper panels), T = 40 MeV (three middle panels), and T = 60 MeV (three 
lower panels). The diquark coupling strength is Gd = Gs- 



In order to understand this, let me inspect the various quantities at T = which are displayed 
in the upper three panels of Fig. 12.141 At the chiral phase boundary, the up and down quark 
masses become relatively small, whereas the strange quark mass experiences only a moderate 
drop of about 84 MeV induced by the 't Hooft interaction. This is not sufficient to populate any 
strange quark states at the given chemical potential, and the system mainly consists of up and 
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down quarks together with a small fraction of electrons, see Fig. 12.161 The electric charge chemical 
potential which is needed to maintain neutrality in this regime is between about —73 MeV and 
—94 MeV. It turns out that the resulting splitting of the up and down quark Fermi momenta is 
too large for the given diquark coupling strength to enable diquark pairing and the system stays 
in the normal quark phase. 
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Figure 2.16: 



The dependence of the number densities of quarks and electrons on the quark chemical 
potential at T — MeV for the diquark coupling strength Gd — jGs- Note that the 
densities of all three quark flavors coincide above /j, = 457 MeV. The density of muons 
vanishes for all values of /j. 



At n w 432 MeV, the chemical potential felt by the strange quarks, /i — A*q/3, reaches the 
strange quark mass and the density of strange quarks becomes nonzero. At first, this density is 
too small to play a sizeable role in neutralizing the quark matter, or in enabling strange-nonstrange 
cross- flavor diquark pairing, see Fig. 12. 161 The first instability happens at /x g cFL ~ 443 MeV, where 
a first-order phase transition from the NQ phase to the gCFL phase takes place. This is directly 
related to a drop of the strange quark mass by about 121 MeV. As a consequence, strange quarks 
become more abundant and pairing gets easier. Yet, in the gCFL phase, the strange quark mass is 
still relatively large, and the standard BCS pairing between strange and light (i.e., up and down) 
quarks is not possible. In contrast to the regular CFL phase, the gCFL phase requires a nonzero 
density of electrons to stay electrically neutral. At T = 0, therefore, one could use the value of 
the electron density as a formal order parameter that distinguishes these two phases [90,91]. 

With increasing the chemical potential further (still at T = 0) , the strange quark mass decreases 
and the cross-flavor Cooper pairing gets stronger. Thus, the gCFL phase eventually turns into 
the regular CFL phase at /^cfl ~ 457 MeV. The electron density goes to zero at this point, as it 
should. This is indicated by the vanishing value of fiQ in the CFL phase, see the upper right panel 
in Fig. 12.141 I remind that the CFL phase is neutral because of having equal number densities of 
all three quark flavors, n u = rid — n s , see Figs. 12.161 and [2~T71 This equality is enforced by the 
pairing mechanism, and this is true even when the quark masses are not exactly equal [20]. 
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Figure 2.17: The dependence of the number densities of quarks, electrons, and muons on the quark 
chemical potential at T = MeV for the diquark coupling strength Gd = Gs- Note, that 
the densities of all three quark flavors coincide above fi — 414 MeV. 



The same NJL model at zero temperature was studied previously in Ref. [133]. My results 
agree qualitatively with those of Ref. [133] only when the quark chemical potential is larger than 
the critical value for the transition to the CFL phase at 457 MeV. The appearance of the gCFL 
phase for 443 MeV < fi < 457 MeV was not recognized in Ref. [133]. Instead it was suggested 
that there exists a narrow (about 12 MeV wide) window of values of the quark chemical potential 
around // w 450 MeV in which the 2SC phase is the ground state. By carefully checking the same 
region, I find that the 2SC phase does not appear there. 

This is illustrated in Fig. I2.18l where the pressure of three different solutions is displayed. Had 
I ignored the gCFL solution (thin solid line), the 2SC solution (dashed line) would indeed be the 
most favored one in the interval between fi w 445 MeV and \i w 457 MeV. After including the 
gCFL phase in the analysis, this is no longer the case. 

Now let me turn to the case of nonzero temperature. One might suggest that this should be 
analogous to the zero temperature case, except that Cooper pairing is somewhat suppressed by 
thermal effects. In contrast to this naive expectation, the thermal distributions of quasiparticles 
together with the local neutrality conditions open qualitatively new possibilities that were absent 
at T = 0. As in the case of the two-flavor model of Refs. [84, 85], a moderate thermal smearing 
of mismatched Fermi surfaces could increase the probability of creating zero-momentum Cooper 
pairs without running into a conflict with Pauli blocking. This leads to the appearance of several 
stable color-superconducting phases that could not exist at zero temperatures. 

With increasing the temperature, the first qualitatively new feature in the phase diagram 
appears when 5 MeV < T < 10 MeV. In this temperature interval, the NQ phase is replaced by 
the uSC phase when the quark chemical potential exceeds the critical value of about 444 MeV. The 
corresponding transition is a first-order phase transition, see Fig. 12.121 Increasing the chemical 
potential further by several MeV, the uSC phase is then replaced by the gCFL phase, and the gCFL 
phase later turns gradually into the (m)CFL phase. (In this study, I do not distinguish between 
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the CFL phase and the mCFL phase [110].) Note that, in the model at hand, the transition 
between the uSC and the gCFL phase is of second order in the following two temperature intervals: 
5 MeV < T < 9 MeV and T < 24 MeV. On the other hand, it is a first-order phase transition when 
9 MeV < T < 24 MeV. Leaving aside its unusual appearance, this is likely to be an accidental 
property in the model for a given set of parameters. For a larger value of the diquark coupling, in 
particular, such a feature does not appear, see Fig. 12.131 
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Figure 2.18: The pressures of different phases of neutral color-superconducting quark matter divided by 
H 4 as a function of the quark chemical potential fj, at T = 0: regular CFL phase (bold solid 
line), gap less CFL phase (thin solid line), 2SC phase (dashed line), normal quark phase 
(dotted line). The diquark coupling strength is Go — jGs- 



The transition from the gCFL to the CFL phase is a smooth crossover at T > 0. The reason is 
that the electron density is not a good order parameter that could be used to distinguish the gCFL 
from the CFL phase when the temperature is nonzero. This is also confirmed by my numerical 
results for the electric charge chemical potential [iq in Fig. 12.141 While at zero temperature the 
value of /Iq vanishes identically in the CFL phase, this is not the case at nonzero temperatures. 

Another new feature in the phase diagram appears when the temperature is above about 
11 MeV. In this case, with increasing the quark chemical potential, the Cooper instability happens 
immediately after the %SB phase. The corresponding critical value of the quark chemical potential 
is rather low, about 365 MeV. The first color-superconducting phase is the g2SC phase [84,85]. 
This phase is replaced with the 2SC phase in a crossover transition only when /j, > 445 MeV. The 
2SC phase is then followed by the gapless uSC (guSC) phase, by the uSC phase, by the gCFL 
phase and, eventually, by the CFL phase (see Fig. I2.12p . 

In the NJL model at hand, determined by the parameters in Eq. (|2.119p . I do not find the 
dSC phase as the ground state anywhere in the phase diagram. This is similar to the conclusion 
of Refs. [110, 113], but differs from that of Refs. [Ill, 112]. This should not be surprising because, 
as was noted earlier [113], the appearence of the dSC phase is rather sensitive to a specific choice 
of parameters in the NJL model. 
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The phase diagram in Fig. 12. 121 has a very specific ordering of quark phases. One might ask if 
this ordering is robust against the modification of the parameters of the model at hand. Below, I 
argue that some features are indeed quite robust, while others are not. 

It should be clear that the appearance of color-superconducting phases under the stress of 
neutrality constraints is very sensitive to the strength of diquark coupling. In the case of two- 
flavor quark matter, this was demonstrated very clearly in Refs. [84, 85] at zero as well as at 
nonzero temperatures. Similar conclusions were also reached in the study of three-flavor quark 
matter at zero temperature [117]. 

In the model at hand, it is instructive to study the phase diagram in the regime of strong diquark 
coupling, Gd = Gs- The corresponding results are summarized in the diagram in Fig. 12.131 As 
one can see, the main qualitative difference between the diagrams in Figs. 12. 121 and l2~T"3l occurs at 
intermediate values of the quark chemical potential. While at Gd = f Gs, there is a large region 
of the g2SC phase sandwiched between the low-temperature and high-temperature NQ phases, 
this is not the case at stronger coupling, Go = Gs- The regions of the gapless phases shrink with 
increasing diquark coupling constant. Above some value of the diquark coupling constant in the 
regime of (very) strong coupling, all gapless phases disappear [118]. 

The last observation can easily be explained by the fact that with increasing diquark coupling 
strength, the condensation energy also increases and therefore Cooper pairing is favorable, even if 
there is a larger mismatch of the Fermi surfaces due to charge neutrality constraints. Moreover, 
in the presence of large gaps, the Fermi surfaces are smeared over a region which is of the order of 
the values of the gaps. Therefore additional thermal smearing is of no further help, and it is not 
surprising that the thermal effects in a model with sufficiently strong coupling are qualitatively the 
same as in models without neutrality constraints imposed: thermal fluctuations can only destroy 
the pairing. In the model with a not very strong coupling, on the other hand, the interplay of the 
charge neutrality and thermal fluctuations is more subtle. The normal phase of cold quark matter 
develops a Cooper instability and becomes a color superconductor only after a moderate thermal 
smearing of the quark Fermi surfaces is introduced [84,85]. 

Other than this, the qualitative features of the phase diagrams in Figs. 12.121 and 12.131 are 
similar. Of course, in the case of the stronger coupling, the critical lines lie systematically at 
higher values of the temperature and at lower values of the quark chemical potential. In this 
context one should note that the first-order phase boundary between the two normal regimes %SB 
and NQ is insensitive to the diquark coupling. Therefore, upon increasing Gd it stays at its place 
until it is eventually displaced by the expanding 2SC phase. As a result, there is no longer a 
critical endpoint in Fig. 12.131 but only a critical point where the first-order normal(xSB)-2SC 
phase boundary changes into second order. 

2.3 The phase diagram with the effect of neutrino trapping 

Quark matter is expected to be color-superconducting. At extremely large densities, namely 
when the quark chemical potential [i is much larger than the constituent, medium- modified quark 
masses, the ground state of matter is given by the CFL phase [18] (for studies of QCD at large 
densities, see also Refs. [19,109]). At the highest densities existing in stars, however, the chemical 
potential is unlikely to be much larger than 500 MeV, while the constituent mass of the strange 
quarks is not smaller than the current mass, which is about 100 MeV. In stellar matter, therefore, 
the heavier strange quarks may not be able to participate in diquark Cooper pairing as easily 
as the light up and down quarks. Then, the pairing of light quarks can lead to the two-flavor 
color-superconducting ground state [15,32,107,108]. It should be pointed out, however, that the 
2SC phase is subject to large penalties after imposing the charge neutrality and /3-equilibrium 
conditions [83]. Indeed, when the fraction of strange quarks and leptons is small, the electric 
neutrality requires roughly twice as many down quarks as up quarks. In this case, Cooper pairing of 
up quarks with down quarks of opposite momenta becomes rather difficult. Then, depending on the 
details of the interaction, the gCFL phase [90,91], the g2SC phase [84,85], or even the normal quark 
matter phase (NQ) may be more favored. Also, there exist other reasonable possibilities [35,36,39, 
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94,96,99-101,120,140,141] which, in view of the known instabilities in the gapless phases [92,94,97], 
are considered to be very promising. 

In Ref. [114], the phase diagram of neutral quark matter was obtained which was described 
by a NJL-type model with the parameter set from Ref. [136]. This work is presented in Sec. 12.21 
In contrast to previous studies in Refs. [110,112,113], dynamically generated quark masses were 
treated self-consistently in Ref. [114]. (For earlier studies on color superconductivity, treating 
quark masses as dynamical quantities, see also Refs. [139,142].) 

In this section, I follow the same approach to study the effect of a nonzero neutrino (or, more 
precisely, lepton- number) chemical potential on the structure of the phase diagram [116]. This is 
expected to have a potential relevance for the physics of protoneutron stars where neutrinos are 
trapped during the first few seconds of the stellar evolution. In application to protoneutron stars, 
it is of interest to cover a range of parameters that could provide a total lepton fraction in quark 
matter of up to about 0.4. This is the value of the lepton-to-baryon charge ratio in iron cores of 
progenitor stars. Because of the conservation of both, lepton and baryon charges, this value is also 
close to the lepton fraction in protoneutron stars at early times, when the leptons do not have a 
chance to diffuse through dense matter and escape from the star. 

The effect of neutrino trapping on color-superconducting quark matter has been previously 
discussed in Ref. [133]. It was found that a nonzero neutrino chemical potential favors the 2SC 
phase and disfavors the CFL phase. This is not unexpected because the neutrino chemical potential 
is related to the conserved lepton number in the system and therefore it also favors the presence 
of (negatively) charged leptons. This helps 2SC-type pairing because electrical neutrality in quark 
matter can be achieved without inducing a very large mismatch between the Fermi surfaces of 
up and down quarks. The CFL phase, on the other hand, is electrically and color neutral in the 
absence of charged leptons when T = [20]. A nonzero neutrino chemical potential can only spoil 
CFL-type pairing. 

In this section, I extend the analysis of Ref. [133] by performing a more systematic survey of 
the phase diagram in the space of temperature, quark, and lepton-number chemical potentials. 
This also includes the possibility of gapless phases which have not been taken into account in 
Ref. [133]. 

2.3.1 Model 

I consider a system of up, down, and strange quarks, in weak equilibrium with charged leptons 
and the corresponding neutrinos. I assume that the lepton sector of the model is given by an ideal 
gas of massive electrons (m e = 0.51099906 MeV) and muons (m M = 105.658389 MeV) [16], as 
well as massless electron and muon neutrinos. I do not take into account the r lepton, which is 
too heavy to play any role. Also, I neglect the possibility of neutrino mixing, and therefore I do 
not take into account the r neutrino either. In the quark sector, I use the same three-flavor NJL 
model as in SecjUH cf. Eq. (|2~TT5l) . 

In the model at hand, there are six mutually commuting conserved charge densities. They split 
naturally into the following three classes: 

electric charge: This is related to the U(l) symmetry of electromagnetism. The corresponding 
charge density is given by 

nq = (^QV) - n e - n M , (2.169) 

where Q is the electric charge matrix of the quarks (|1.50j) , and n e and denote the number 
densities of electrons and muons, respectively, cf. Eq. (I1.53|) . 

two lepton charges: As long as the neutrinos are trapped and their oscillations are neglected, 
the lepton family numbers are conserved. The corresponding densities read, 



n Le = n e + n Ve , n L ^ =n ll + n Vtl , 



(2.170) 
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baryon number and two color charges: The SU(3)-co\or symmetry and U(l) baryon number 
symmetry imply the conservation of three independent charge densities in the quark sector, 



n 3 



n 8 = (ip^Tstp) , 



(2.171) 



where T 3 and Tg are the matrices associated with the two mutually commuting color charges 
of the [SU(3) C ] gauge group, cf. Eq. (|1.57p . Note that n is the quark number density that 
is related to the baryon number density as follows: n B = n/3. An alternative choice of the 
three conserved charges is given by the number densities of red, green and blue quarks, i.e., 



1 

n r = n B + n 3 + ~/=™8 



n B -n 3 + , 



n b = n B - ~/= n 8 • (2.172) 



The six conserved charge densities defined above are related to the six chemical potentials of the 
model. These are the quark chemical potential /1 = /is/3, the two color chemical potentials \x 3 
and us, the electric charge chemical potential [1q, and the two lepton- number chemical potentials 
^ Le and fi Lfi . 

In chemical equilibrium, the chemical potentials of all individual quark and lepton species can 
be expressed in terms of these six chemical potentials according to their content of conserved 
charges. For the quarks, which carry quark number, color and electric charge, this is related 
to the matrix of the quark chemical potentials (|1.52|) . The neutrinos, on the other hand, carry 
lepton number only, cf. Eq. (|1.43p . Finally, electrons and muons carry both, lepton number and 
electric charge, see Eq. (|1.47[) . As in Sec. 12.21 the quark part of the model is treated in the mean- 
field (Hartree) approximation, allowing for the presence of both, quark-antiquark condensates and 
scalar diquark condensates. By using the same gap ansatz (|2.126p ~ (|2.127p and by following the 
same steps in the derivation as in Sec. 12.21 and by including the ideal-gas contribution for the 
leptons, I arrive at the following expressions for the pressure at nonzero and zero temperature, 
respectively: 
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where are eighteen independent positive-energy eigenvalues, see Sec. 12.21 for details. In order 
to obtain the values for the quark-antiquark condensates a a and the color-superconducting gap 
parameters A c , I solve the following six stationary conditions: 



dp 

d<J a 







dp 

dKc 



= . 



(2.174) 



In order to enforce the conditions of local charge neutrality in quark matter, one also requires 
Eqs. (|1.54p and (| 1 .60(1 to be satisfied. By solving these, I determine the values of the three 
corresponding chemical potentials fiQ, [i 3 , and for a given set of the other chemical potentials, 
/i, /ii c , 1-tLui an d for a given temperature T. In general, therefore, the phase diagram of dense 
quark matter with neutrino trapping should span a four-dimensional parameter space. 
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Note that instead of using the chemical potentials, /x, (J-l c , and /xj, as free parameters in the 
study of the phase diagram, one may also try to utilize the quark number density and the two 
lepton fractions, 

n = — , Y Lb = 3 , Y L = 3 — , (2.175) 

o/i n n 

where the two lepton densities arc defined by 

- d P _ dp 
n Le = 75 , n L = . (2.176 

In some cases, the choice of n, Yi, e , and Yl as free parameters is indeed very useful. For instance, 
this is helpful in order to determine the initial state of matter inside protoneutron stars at very 
early times, when the lepton fractions are approximately the same as in the progenitor stars (i.e., 
Yjj e ~ 0.4 and Yl = 0). The problem is, however, that such an approach becomes ambiguous in 
the vicinity of first-order phase transitions, where the baryon number density as well as the lepton 
fractions are in general discontinuous. For this reason, it is more appropriate to study the phase 
structure of (dense) QCD at given fixed values of the chemical potentials /x, /xl c , and /it . Unlike 
densities, the chemical potentials change continuously when the system crosses a boundary of a 
first-order phase transition. (It should be noted that the chemical potentials /xq, /X3, and /xs may 
change discontinuously at a boundary of a first-order phase transition because of the long-range 
Coulomb interaction enforcing the constraints uq = 0, 123 = 0, and rig = 0.) 



2.3.2 Simplified considerations 

As mentioned before, neutrino trapping favors the 2SC phase and strongly disfavors the CFL 
phase [133]. This is a consequence of the modified /3-equilibrium condition in the system. In this 
section, I would like to emphasize that this is a model-independent effect. In order to understand 
the physics behind it, it is instructive to start my consideration from a very simple toy model. 
Later, many of its qualitative features will be also observed in my self-consistent numerical analysis 
of the NJL model. 

Let me first assume that strange quarks are very heavy and consider a gas of non-interacting 
massless up and down quarks in the normal quark phase at T = 0. As required by /3 equilibrium, 
electrons and electron neutrinos are also present in the system. (Note that in this section I neglect 
muons and muon neutrinos for simplicity.) 

In the absence of Cooper pairing, the densities of quarks and leptons are given by 

cf. Eq. (|1.76| . Expressing the chemical potentials through /x, \iq and /xl c by using Eq. ()1.49|) . 
and imposing electric charge neutrality by setting Eq. (|2.169[) equal to zero, one arrives at the 
following relation: 

2 (1 + \yf - (1 - \vf - (x - y) 3 = , (2.178) 

where I have introduced the chemical potential ratios x = fi^^/fj, and y = (J-q/h- The above cubic 
equation can be solved for y (electric chemical potential) at any given x (lepton-number chemical 
potential) with Cardano's formulae which are shown in Sec. IB.9l in the Appendix. The result can 
be used to calculate the ratio of quark chemical potentials, /Xd//x« = (3 — y)/(3 + 2y). 

The ratio /Xd//x„ as a function of /Xi e //x is shown in Fig. 12.191 At vanishing /xj, e , one finds 
y w —0.219 and, thus, /itrf//x u w 1.256 (note that this value is very close to 2 1 / 3 sa 1.26). This 
result corresponds to the following ratios of the number densities in the system: n u /rid ~ 0.504 
and n e /rid ss 0.003, reflecting that the density of electrons is tiny and the charge of the up quarks 
has to be balanced by approximately twice as many down quarks, cf. Sec. 11.4.41 
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At fiL e = fi, on the other hand, the real solution to Eq. ()2. 178|) is y — 0, i.e., the up and 
down Fermi momenta become equal. This can be seen most easily if one inverts the problem and 
solves Eq. (|2.178|) for x at given y. When y = one finds x = 1, meaning that fid = M« an d, m 
turn, suggesting that pairing between up and down quarks is unobstructed at fiL e — M- This is in 
contrast to the case of vanishing fJ,L e , when the two Fermi surfaces are split by about 25%, and 
pairing is difficult. 

It is appropriate to mention that many features of the above considerations would not change 
much even when Cooper pairing is taken into account. The reason is that the corresponding 
corrections to the quark densities are parametrically suppressed by a factor of order (A//z) 2 , 
where A is the gap parameter. 

In order to estimate the magnitude of the effect in the case of quark matter in protoneutron 
stars, I indicate several typical values of the lepton fractions Yl b in Fig. 12.191 As mentioned earlier, 
Yl c is expected to be of order 0.4 right after the collapse of the iron core of the progenitor star. 
According to Fig. 12.191 this corresponds to fid/ i^u ~ 1-1, i.e., while the splitting between the up 
and down Fermi surfaces does not disappear completely, it gets reduced considerably compared 
to its value in the absence of trapped neutrinos. This reduction substantially facilitates the cross- 
flavor pairing of up and down quarks. The effect is gradually washed out during about a dozen of 
seconds of the deleptonization period when the value of Yj Je decreases to zero. 

The toy model is easily modified to the opposite extreme of three massless quark flavors, where 
the number density of strange quarks reads, 




(2.179) 



Basically, this corresponds to replacing Eq. (|2.178[) by 



2(l + ^) 3 -2(l-ly) 3 ~ (x-yf = 



(2.180) 
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In the absence of neutrino trapping, x — 0, the only real solution to this equation is y = 0, 
indicating that the chemical potentials (which also coincide with the Fermi momenta) of up, down 
and strange quarks are equal, cf. Sec. 11.4.41 This reflects the fact that the system with equal 
densities of up, down and strange quarks is neutral by itself, without electrons. With increasing 
x oc /i£ e , the solution requires a nonzero y cx [Xq, suggesting that up-down and up-strange pairing 
becomes more difficult. To see this more clearly, I can go one step further in the analysis of the 
toy model. 

Let me assume that the quarks are paired in a regular, i.e., fully gapped, CFL phase at T = 0. 
Then, as shown in Ref. [20], the quark part of the matter is automatically electrically neutral. 
Hence, if I want to keep the whole system electrically and color neutral, there must be no electrons. 
Obviously, this is easily realized without trapped neutrinos by setting \iq equal to zero. At non- 
vanishing /xx, e the situation is more complicated. The quark part is still neutral by itself and 
therefore no electrons are admitted. Hence, the electron chemical potential fi e = — /iq must 
vanish, and consequently (1q should be nonzero and equal (ii, e . It is natural to ask what should 
be the values of the color chemical potentials /13 and /xs in the CFL phase when fii, e =/= 0. 

In order to analyze the stress on the CFL phase due to nonzero (iL e , I follow the same approach 
as in Sec. 11.5.11 In this analytical consideration, I also account for the effect of the strange quark 
mass simply by shifting the strange quark chemical potential by —Mg/(2fi). In my notation, 
CFL- type pairing requires the Fermi momenta defined by Eq. (|1.96|) . These are used to calculate 
the pressure in the toy model, 



Ptoy 



a—r a—u 



where in contrast to the toy model for the CFL phase defined by Eq. (|1.95[) . the contributions of 
electrons and electron neutrinos are included and the bag pressure is neglected. With Eqs. (|1.98[) 
and (|2.181|) . one easily derives the neutrality conditions (|1.54p and p.60|) . Thus, it becomes 
obvious that charge neutrality requires /ig = The neutrality condition, 713 = 0, requires 

that fig, = —fiQ which means that ^3 = — fiL e - Finally, one can check that the third neutrality 
condition, ns = 0, requires, 

7f-f- <" 82 > 

The results for the charge chemical potentials (iq, /13, and /is imply the following magnitude of 
stress on pairing in the CFL phase: 

5»(rd,gu) = ^-y^ = ML, , (2.183a) 

o^{ gs ,bd) - 2 ~ "2/7 ' (■2.186c) 

Note that there is no mismatch between the values of the chemical potentials of the other three 
quarks, yT u = l4=j4 = M ~ ^ S 2 /(M- 

From Eq. (|2. 183ft one can see that the largest mismatch occurs in the (rs, bu) pair (for positive 
/i£ e ). The CFL phase can withstand the stress only if the value of 5^( rs .bu) is l ess than A2. A 
larger mismatch should drive a transition to a gapless phase exactly as in Refs. [84,85,90,91]. 
Thus, the critical value of the lepton-number chemical potential is, 

/4-)«A 2 -|j. (2.184) 

When (iL e > /i^ cl \ the CFL phase turns into the gCFL' phase, which is a variant of the gCFL 
phase [90,91]. By definition, the gapless mode with a linear dispersion relation in the gCFL' phase 
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is rs—bu instead of gs-bd as in the standard gCFL phase. (Let me remind that the mode aa-b[3 is 
defined by its dispersion relation which interpolates between the dispersion relations of hole-type 
excitations of aa-quarks at small momenta, k -C //^, and particle- type excitations of 6/3-quarks at 



In order to see what this means for the physics of protoneutron stars, I should again try to 
relate the value of fiL c to the lepton fraction. There are no electrons in the (regular) CFL phase at 
T = 0. Therefore the entire lepton number is carried by neutrinos. For the baryon density I may 
neglect the pairing effects to first approximation and employ the ideal-gas relations. This yields, 



Inserting typical numbers, y, = 500 MeV and y Lc < A « 50 MeV-100 MeV, one finds Y Le < 10" 4 - 
10 -3 . Thus, there is practically no chance to find a sizeable amount of leptons in the CFL phase. 
The constraint gets relaxed slightly at nonzero temperatures and/or in the gCFL phase, but the 
lepton fraction remains rather small even then (my numerical results indicate that, in general, 
Y Le < 0.05 in the CFL phase). 



The simple toy-model considerations in Sec. 12.3.21 give a qualitative understanding of the effect 
of neutrino trapping on the mismatch of the quark Fermi momenta and, thus, on the pairing 
properties of two- and three-flavor quark matter. Now I turn to a more detailed numerical analysis 
of the phase diagram in the framework of the NJL model defined in Sec. 12.3.11 

In the numerical calculations, I use the same set of model parameters as in Sec. 12.21 of my 
thesis, sec Eq. (|2.119[) . The parameters are chosen to reproduce several key observables of vacuum 
QCD [136]. In this section, I choose a diquark coupling constant Go — jGs- 

In order to obtain the phase diagram, one has to determine the ground state of matter for 
each given set of the parameters. As discussed in Sec. 12.3. 1[ in the case of locally neutral matter 
with trapped neutrinos, there are four parameters that should be specified: the temperature T, the 
quark chemical potential /x as well as the two lepton family chemical potentials /Lt£ e and . After 
these are fixed, the values of the pressure in all competing neutral phases of quark matter should 
be compared. This is determined by using the same algorithm as in Sec. 12.21 The complete set of 
equations (|1.54|) . (11.60|) . and (|2.174|) is solved for each of the eight phases allowed by symmetries. 
Then, the corresponding values of the pressure are determined from Eq. (|2. 173(1 . The phase with 
the largest pressure is the ground state. 

In this section, I always assume that the muon lepton-number chemical potential vanishes, 
i.e., Hl = 0. This is expected to be a good approximation for matter inside protoneutron stars. 
My analysis can thus be interpreted as an extension of the T-fj, phase diagram which is discussed 
in Sec. 12.21 into the fj,j Je direction. Consequently, the complete phase structure requires a three- 
dimensional presentation. 

The three-dimensional phase diagram 

The general features of the phase diagram in the three-dimensional space, spanned by the quark 
chemical potential /i, the lepton-number chemical potential /Lt/, e , and the temperature T are de- 
picted in Fig. 12.201 Because of the rather complicated structure of the diagram, only the four main 
phases (xSB, NQ, 2SC, and CFL) are shown explicitly. Although it is not labeled, a thin slice of 
a fifth phase, the uSC phase, squeezed in between the 2SC and CFL phases, can also be seen. In 
addition, at small values of the lepton-number chemical potential and small temperatures, another 
unmarked region between /i ~ 370 MeV and fi ~ 440 MeV exists. In Fig. 12.201 however, only a 
small part of the surface at its boundary with the 2SC phase can be seen. This region corresponds 
to the second piece of the normal quark matter phase which is disconnected from the main region 
(the appearance of a small disconnected region can also be deduced from the T-[i phase diagram 



large momenta, k ^> /Jp.) 




(2.185) 



2.3.3 Results 



90 



The phase diagram of neutral quark matter 



in Fig. I2.12j) . While lacking detailed information, the phase diagram in Fig. 12.201 gives a clear 
overall picture. Among other things, one sees, for example, that the CFL phase becomes strongly 
disfavored with increasing /j,L e and gets gradually replaced by the 2SC phase. 




Figure 2.20: The phase diagram of neutral three-flavor quark matter in the three-dimensional space 
spanned by the quark chemical potential fi, the lepton-number chemical potential /i£, e , and 
the temperature T. 



In order to discuss the structure of the phase diagram in more details I proceed by showing 
several two-dimensional slices of it. These are obtained by keeping one of the chemical potentials, 
fi or fixed and varying the other two parameters. 



The T fj, phase diagram 

The phase diagrams at two fixed values of the lepton-number chemical potential, fiL e = 200 MeV 
and fiL e = 400 MeV are presented in Figs. 12.211 and 12.221 The phase diagram at fiL e — MeV, 
see Fig. 12.121 has already been discussed in Sec. 12.21 The general effects of neutrino trapping can 
be understood by analyzing the similarities and differences between these three phase diagrams. 
In this section, I use the same convention for line styles as in Sec. 12.21 thick and thin solid 
lines denote first-order and second-order phase transitions, respectively; dashed lines indicate the 
(dis-) appearance of gapless modes in different phases. 

Here it is appropriate to note that, in the same model, a schematic version of the T—fj, phase 
diagram at Hl c = 200 MeV was first presented in Ref. [133], see the right panel of Fig. 4 there. 
If one ignores the complications due to the presence of the uSC phase and various gapless phases, 
the results of Ref. [133] are in qualitative agreement with my findings. 
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Figure 2.21: The phase diagram of neutral quark matter at a fixed lepton-number chemical potential 
fl Le = 200 MeV. 
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Figure 2.22: The phase diagram of neutral quark matter at a fixed lepton-number chemical potential 
fi Lc = 400 MeV. 



In order to understand the basic characteristics of different phases in the phase diagrams in 
Figs . 12 . 2 H and ET^l I also present the results for the dynamical quark masses, the gap parameters, 
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and the charge chemical potentials. These are plotted as functions of the quark chemical potential 
in Figs. 12.231 and |2~2~4"1 for two different values of the temperature in the case of hl c = 200 MeV 
and /iLe = 400 MeV, respectively. 




Figure 2.23: Dependence of the quark masses, of the gap parameters, and of the electric and color charge 
chemical potentials on the quark chemical potential at a fixed temperature, T = MeV 
(three upper panels) and T — 40 MeV (three lower panels). The lepton-number chemical 
potential is kept fixed at = 200 MeV. 



In each of the diagrams, there are roughly four distinct regimes. At low temperature and 
low quark chemical potential, there is a region in which the approximate chiral symmetry is 
spontaneously broken by large ^-^-condensates. The corresponding phase is denoted by x^B. In 
this regime, quarks have relatively large constituent masses which are close to the vacuum values, 
see Figs. 12.231 and 12.241 Here, the density of all quark flavors is very low and even vanishes at 
T = 0. There is no diquark pairing in this phase. The xSB phase is rather insensitive to the 
presence of a nonzero lepton-number chemical potential. With increasing /i^ e the phase boundary 
is only slightly shifted to lower values of fi. This is just another manifestation of the strengthening 
of the 2SC phase due to neutrino trapping. 

With increasing temperature, the ^^-condensates melt and the %SB phase turns into the NQ 
phase where the quark masses are relatively small. Because of the explicit breaking of the chiral 
symmetry by the current quark masses, there is no need for a phase transition between the two 
regimes. In fact, at low chemical potentials, I find only a smooth crossover, whereas there is a 
first-order phase transition in a limited region, 320 MeV < \i < 360 MeV. In contrast to the 
xSB regime, the high-temperature NQ phase extends to arbitrary large values of fj,. All these 
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qualitative features are little affected by the lepton-number chemical potential. 

The third regime is located at relatively low temperatures but at quark chemical potentials 
higher than in the x^B phase. In this region, the masses of the up and down quarks have already 
dropped to values well below their respective chemical potentials while the strange quark mass 
is still large, see left columns of panels in Figs. 12.231 and 12.241 As a consequence, up and down 
quarks are quite abundant but strange quarks are essentially absent. 
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Figure 2.24: Dependence of the quark masses, of the gap parameters, and of the electric and color charge 
chemical potentials on the quark chemical potential at a fixed temperature, T = MeV 
(three upper panels) and T — 40 MeV (three lower panels). The lepton-number chemical 
potential is kept fixed at = 400 MeV. 



It turns out that the detailed phase structure in this region is very sensitive to the lepton- 
number chemical potential. At fj,L e — 0, as discussed in Sec. l2.2| the pairing between up and down 
quarks is strongly hampered by the constraints of neutrality and (5 equilibrium. As a consequence, 
there is no pairing at low temperatures, T < 10 MeV, and g2SC-type pairing appears at moderate 
temperatures, 10 MeV < T < T c with the value of T c in the range of several dozen MeV, see 
Fig. 12.121 The situation changes dramatically with increasing the value of the lepton-number 
chemical potential. Eventually, the low-temperature region of the normal phase of quark matter 
is replaced by the (g)2SC phase (e.g., at fi = 400 MeV, this happens at /x Lc ~ 110 MeV). With 
fiL e increasing further, no qualitative changes happen in this part of the phase diagram, except 
that the area of the (g)2SC phase expands slightly. 

Finally, the region in the phase diagram at low temperatures and large quark chemical poten- 
tials corresponds to phases in which the cross-flavor strange-nonstrange Cooper pairing becomes 
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possible. In general, as the strength of pairing increases with the quark chemical potential, the 
system passes through regions of the gapless uSC (guSC), uSC, and gCFL phases and finally 
reaches the CFL phase. (Of course, the intermediate phases may not always be realized.) The 
effect of neutrino trapping, which grows with increasing lepton-number chemical potential, is to 
push out the location of the strange-nonstrange pairing region to larger values of /i. Of course, 
this is in agreement with the general arguments in Sec. 12.3.21 
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Table 2.2: The classification of gapless phases in color-superconducting quark matter. The unmarked 
gap parameters (— ) are zero while the checkmarked gap parameters (/) are nonzero in the 
respective gapless color-superconducting quark phases. The gapless modes aa-b/3 are defined 
by their dispersion relations which interpolate between the dispersion relations of hole-type 
excitations of aa-quarks at small momenta, k <C ^ , and particle- type excitations of fe/3-quarks 
at large momenta, k S> fJ^- 
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Figure 2.25: Dependence of the electron family lepton fraction Yl c for = 200 MeV (thick lines) 
and jUi e = 400 MeV (thin lines) on the quark chemical potential at a fixed temperature, 
T = MeV (solid lines), T = 20 MeV (dashed lines), and T = 40 MeV (dotted lines). 



It is interesting to note that the growth of the strangeness content with increasing quark chemi- 
cal potential could indirectly be deduced from the behavior of the electric charge chemical potential 
Hq at T — 0, see the solid lines in the right panels of Figs. 12.231 and [2~2"4l The value of fiQ reaches 
its minimum somewhere in a range of values of the quark chemical potential around /i ~ 440 MeV. 



2.3 The phase diagram with the effect of neutrino trapping 



95 



This corresponds to the point where the strange quark chemical potential, ji a ~ fi — ^q/3, reaches 
the value of the strange quark mass (see left panels). Hence, there are essentially no strange 
quarks at lower values of /x, and a rapidly increasing amount of strange quarks at higher values of 
fi. Since the latter contribute to the electric neutralization, this is a natural explanation for the 
drop of |/iq| above this point. 
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Figure 2.26: 




1 50 200 250 
H, [MeV] 



400 



The phase diagram of neutral quark matter in the plane of temperature and lepton-number 
chemical potential at a fixed value of the quark chemical potential, /i = 400 MeV. This 
phase diagram corresponds to the outer stellar core. 



As I mentioned earlier, the presence of the lepton-number chemical potential //L e leads to a 
change of the quark Fermi momenta. This change, in turn, affects Cooper pairing of quarks, 
facilitating the appearance of some phases and suppressing others. As it turns out, there is also 
another qualitative effect due to a nonzero value of fiL e - In particular, I find several new variants 
of gapless phases which do not exist at vanishing /i£ e . In Figs. 12.21) 12.221 12.26[ and 12.271 these 
are denoted by the same names, g2SC or gCFL, but with one or two primes added. 

I define the g2SC as the gapless two-flavor color-superconducting phase in which the gapless 
excitations correspond to rd-gu and gd-ru modes instead of the usual ru-gd and gu-rd ones, i.e., 
u and d flavors are exchanged as compared to the usual g2SC phase. The g2SC phase becomes 
possible only when the value of {pT u — /i^)/2 = (/iq + Hz)/2 is positive and larger than A3. The 
other phases are defined in a similar manner. In particular, the gCFL' phase, which was already 
introduced in Sec. 12.3.21 is indicated by the gapless rs-bu mode, while the gCFL" phase has both, 
gs-bd (as in the gCFL phase) and rs-bu gapless modes. The definitions of all gapless phases are 
summarized in Table l2~2"l 

The lepton fraction Yj, e 

The numerical results for the lepton fraction Yi, e are shown in Fig. 12.251 The thick and thin lines 
correspond to two different fixed values of the lepton-number chemical potential, /ir, e = 200 MeV 
and hl c — 400 MeV, respectively. For a fixed value of I find that the lepton fraction changes 
only slightly with temperature. This is concluded from the comparison of the results at T = MeV 
(solid lines), T = 20 MeV (dashed lines), and T = 40 MeV (dotted lines) in Fig. EH] 
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Figure 2.27: The phase diagram of neutral quark matter in the plane of temperature and lepton-number 
chemical potential at a fixed value of the quark chemical potential, n — 500 MeV. This 
phase diagram corresponds to the inner stellar core. The triangle denotes the transition 
point from the CFL phase to the gCFL' phase at T — 0. 



As is easy to check, at T — 40 MeV, i.e., when Cooper pairing is not so strong, the fj, dependence 
of Yt e does not differ very much from the prediction in the simple two-flavor model in Sec. I2.3.2I 
By saying this, of course, one should not undermine the fact that the lepton fraction in Fig. I2.25I 
has a visible structure in the dependence on fi at T = MeV and T = 20 MeV. This indicates 
that quark Cooper pairing plays a nontrivial role in determining the value of Yi, e . 

My numerical study shows that it is hard to achieve values of the lepton fraction more than 
about 0.05 in the CFL phase. Gapless versions of the CFL phases, on the other hand, could 
accommodate the lepton fraction up to about 0.2 or so, provided the quark and lepton-number 
chemical potentials are sufficiently large. 

From Fig. 12.251 one can also see that the value of the lepton fraction Yr, c w 0.4, i.e., the value 
expected at the center of the protoneutron star right after its creation, requires a lepton-number 
chemical potential /i£ c in the range somewhere between 200 MeV and 400 MeV, or slightly higher. 
The larger the quark chemical potential /j,, the larger a lepton-number chemical potential 
is needed. Then, in a realistic construction of a protoneutron star, this is likely to result in a 
noticeable gradient of the lepton-number chemical potential at the initial time. This gradient may 
play an important role in the subsequent deleptonization due to neutrino diffusion through dense 
matter. 

The T—fiL e phase diagram 

Now let me explore the phase diagrams in the plane of temperature and lepton-number chemical 
potential, keeping the quark chemical potential fixed. Two such slices of the phase diagram are 
presented in Figs. I2.2B1 and |2~2"TI The first one corresponds to a not very large value of the quark 
chemical potential, fi — 400 MeV. This could be loosely termed as the outer stellar core phase 
diagram. The second one corresponds to fj, = 500 MeV, and one could associate it with the inner 
stellar core case. Note, however, that the terms outer stellar core and inner stellar core should 
not be interpreted literally here. The central densities of (proto-)neutron stars are subject to large 
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theoretical uncertainties and, thus, are not known very well. In the model at hand, the outer 
stellar core case corresponds to a range of densities around 4no, while the inner stellar core case 
corresponds to a range of densities around lOno- These values are of the same order of magnitude 
that one typically obtains in models (see, e.g., Ref. [48]). 
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Figure 2.28: Dependence of the quark masses, of the gap parameters, and of the electric and color charge 
chemical potentials on the temperature at a fixed value of the lepton-number chemical 
potential, fj,L e = MeV (three upper panels), fj,L e = 200 MeV (three middle panels), and 
/iL e = 400 MeV (three lower panels). The quark chemical potential is fj, = 400 MeV (outer 
stellar core case). 
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Figure 2.29: Dependence of the quark masses, of the gap parameters, and of the electric and color charge 
chemical potentials on the temperature at a fixed value of the lepton-number chemical 
potential, fiL c — MeV (three upper panels), ~ 200 MeV (three middle panels), and 
HL e — 400 MeV (three lower panels). The quark chemical potential is [i = 500 MeV (inner 
stellar core case). 



In addition to the phase diagrams, I also show the results for the dynamical quark masses, the 
gap parameters, and the charge chemical potentials. These are plotted as functions of temperature 
in Fig. I2.28I (n — 400 MeV, outer stellar core case) and in Fig. I2.29I (fi = 500 MeV, inner stellar 
core case), for three different values of the lepton-number chemical potential, hl c = MeV (upper 



2.3 The phase diagram with the effect of neutrino trapping 



99 



panels), hl c = 200 MeV (middle panels), and fi^ e = 400 MeV (lower panels). 

At first sight, the two phase diagrams in Figs. 12.261 and l2~27l look so different that no obvious 
connection between them could be made. It is natural to ask, therefore, how such a dramatic 
change could happen with increasing the value of the quark chemical potential from fi — 400 MeV 
to n = 500 MeV. In order to understand this, it is useful to place the corresponding slices of the 
phase diagram in the three-dimensional phase diagram in Fig. 12.201 

The fi = 500 MeV phase diagram corresponds to the right-hand side surface of the bounding 
box in Fig. 12.201 This contains almost all complicated phases with strange-nonstrange cross-flavor 
pairing. The fi — 400 MeV phase diagram, on the other hand, is obtained by cutting the three- 
dimensional phase diagram with a plane parallel to the bounding surface, but going through the 
middle of the phase diagram. This part of the phase diagram is dominated by the 2SC and NQ 
phases. Keeping in mind the general structure of the three-dimensional phase diagram, it is also 
not difficult to understand how the two phase diagrams in Figs. 12.261 and l2~27l transform into each 
other. 

Several comments are in order regarding the zero-temperature phase transition from the CFL to 
gCFL' phase, shown by a small black triangle in the phase diagram at fi = 500 MeV, see Fig. 12.271 
The appearance of this transition is in agreement with the analytical result in Sec. 12.3.21 Moreover, 
the critical value of the lepton-number chemical potential also turns out to be very close to the 
estimate in Eq. I|2.184[) . Indeed, by taking into account that M s « 214 MeV and A 2 « 76 MeV, I 
obtain /i^ r) = A 2 - M^/(2fi) » 30 MeV which agrees well with the numerical value. 

In order to estimate how the critical value of changes with decreasing the quark chemical 
potential below fx = 500 MeV, one can use the zero-temperature numerical results for M s and A 2 
from Sec. 12.21 Then, one arrives at the following power-law fit for the /^-dependence of the critical 
value: 

(Jg* « 0.00575 MeV -1 (p - 457 .4 MeV) (622.1 MeV - fi) , (2.186) 

for 457.4 MeV < fi < 500 MeV. Note that the CFL phase does not appear at T = when 
fi < 457.4 MeV, cf. Fig. 12.121 With increasing the values of the quark chemical potential above 
fi = 500 MeV, one expects that the critical value of fiL e should continue to increase for a while, and 
then decrease when the effects of the cutoff start to suppress the size of the gap A 2 . However, the 
validity of the fit in Eq. (|2. 186|) is questionable there because no numerical data for fi > 500 MeV 
was used in its derivation. 

A schematic version of the phase diagram in the T—fii, e plane was presented earlier in Ref. [133], 
see the left panel in Fig. 4 there. In Ref. [133], the value of the quark chemical potential was 
fi = 460 MeV, and therefore a direct comparison with my results is not easy. One can see, however, 
that the diagram of Ref. [133] fits naturally into the three-dimensional diagram in Fig. 12.201 Also, 
the diagram of Ref. [133] is topologically close to my version of the phase diagram at ft = 500 MeV 
which is shown in Fig. 12.271 The quantitative difference is not surprising: the region of the (g)CFL 
phase is considerably larger at fi = 500 MeV than at fi = 460 MeV. 
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Chapter 3 

Conclusions 



The phase diagram of neutral quark matter was poorly understood as I began with the research 
on this topic in 2003, see Fig. 11.11 The task of my thesis is to illuminate the phase structure of 
neutral quark matter. From the phase diagram of neutral quark matter, one can predict in which 
state is neutral quark matter in the cores of neutron stars which are most probably the candidates 
for neutral color-superconducting quark matter in nature. 

3.1 Summary and discussion 

In Chapter[U I gave an introduction to quark matter, color superconductivity, and stellar evolution. 
I presented the most important color-superconducting quark phases and explained how stars evolve 
and become neutron stars. At extremely large densities which are probably present in the cores 
of neutron stars, there is a dominant attractive interaction between the quarks that causes the 
formation of quark Cooper pairs [10-12]. This causes color superconductivity. Since quarks are 
spin-i fermions which appear in various colors and flavors, there are many possibilities of forming 
quark Cooper pairs. This is why there exist many different color-superconducting quark phases. 

In order to predict which of these phases are energetically preferred in nature, it is necessary to 
compute the phase diagram of neutral quark matter and to determine which phases are dominant 
at a given temperature and quark chemical potential. The phase diagram of neutral quark matter 
consisting of up, down, and strange quarks is presented in Chapter [21 In this context, it is im- 
portant to consider neutral quark matter because otherwise there would occur repulsive Coulomb 
forces which lead to the explosion of the neutron stars. Quark matter has to be color neutral be- 
cause color-charged objects have never been seen in nature. Besides, color-charged neutron stars 
would be unstable. Matter inside neutron stars is in [3 equilibrium which means that neutron-star 
matter is in equilibrium with respect to weak interactions. 

In Sec. l2.1[ I studied massless neutral three-flavor quark matter at large baryon densities within 
an NJL model. The effects of the strange quark mass were incorporated by a shift of the chemical 
potential of strange quarks, see Eq. (I2.35|) . This shift reflects the reduction of the Fermi momenta 
of strange quarks due to their mass. Such an approach is certainly reliable at small values of the 
strange quark mass. In Ref. [113], where the strange quark mass was properly taken into account, I 
confirmed that it is also qualitatively correct at large values of the strange quark mass. I obtained 
a very rich phase structure by varying the strange quark mass, the quark chemical potential, and 
the temperature. 

At T = 0, there are two main possibilities for the strange quark matter ground state: the CFL 
and gCFL phases. These findings confirm the results of Refs. [90, 91] concerning the existence 
of the gCFL phase, the estimate of the critical value of the strange quark mass m s , and the 
dependence of the chemical potentials on m s . I also confirmed that it is the color- neutrality 
condition, controlled by the color chemical potential which drives the transition from the CFL 
to the gCFL phase [90,91]. This is in contrast to gapless two-flavor color superconductivity which 
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results from electrical neutrality [84,85]. 

Because I use a nine-parameter ansatz for the gap matrix, see Eq. (|2.53p . the results of Sec. 12. II 
are more general than those of Refs. [90,91]. For example, I was able to explicitly study the effects 
of the symmetric pairing channel, described by the sextet gap parameters, that were neglected 
in Refs. [90,91]. As one might have expected, these latter modify the quasiparticle dispersion 
relations only slightly. This check was important, however, to see that the zero-temperature phase 
transition from the CFL phase to the gCFL phase, which is not associated with any symmetry, is 
robust against such a deformation of the quark system. 

In Sec. 12. 1[ I also studied the temperature dependence of the gap parameters and the quasi- 
particle spectra. In particular, this study revealed that there exist several different phases of 
neutral three-flavor quark matter that have been predicted in the framework of the Ginzburg- 
Landau-type effective theory in Ref. [111]. My results extend the near-critical behavior discussed 
in Ref. [Ill] to all temperatures. Also, I show how this behavior evolves with changing the value 
of the strange quark mass. The only real qualitative difference between my results and the results 
of Refs. [Ill, 112] is that, instead of the dSC phase, I find the uSC phase in the phase diagram. 
In Ref. [113], I confirmed the results of Sec. [2J] regarding the existence of several different phases 
of neutral three-flavour quark matter at nonzero temperature. I also confirmed the order in which 
they appear. In particular, I observed the appearance of the uSC phase as an intermediate state 
in the melting of the (g)CFL phase into the 2SC phase. Formally, this is different from the pre- 
diction of Ref. [111]. I find, however, that the difference is connected with the choice of the model 
parameters. In the NJL model with a cutoff parameter of 800 MeV used in Ref. [112], there is a 
non-vanishing (although rather small) region of the dSC phase. On the other hand, in the NJL 
model with a relatively small cutoff parameter of 653.3 MeV, no sizeable window of the dSC phase 
can be found. Therefore, I conclude that the appearence of the uSC or dSC phase, repectively, is 
very sensitive to the choice of the model parameters, or more precisely, it depends strongly on the 
cutoff parameter. 

The main result of Sec. 12. H is the complete phase diagram of massless neutral three- flavor quark 
matter in the T-m 2 s / \i and T-\i plane, shown in Fig. 12.111 In this figure, all symmetry- related 
phase transitions are denoted by solid lines. (The symmetries of all phases appearing in this 
figure were discussed in Ref. [111].) The dashed lines in Fig. 12. Ill separate the mCFL and regular 
2SC phases from the gCFL and g2SC phases. These cannot be real phase transitions, but are at 
most smooth crossovers. At T = 0, there is an insulator-metal phase transition between the CFL 
and the gCFL phase [90,91]. At nonzero temperature, there exists a similar insulator-metal-type 
transition between the CFL and the mCFL phase, given by the dotted lines in Fig. 12.111 

In Sec. 12. 2| I studied the T-\i phase diagram of neutral three-flavor quark matter within the 
NJL model of Ref. [136] in which chiral symmetry is broken explicitly by small but nonzero current 
quark masses. As in the previous study in Sec. 12.11 I used the mean-field approximation in the 
analysis. In contrast to Sec. 12.11 in Sec. 12.21 the constituent quark masses were treated self- 
consistently as dynamically generated quantities. The main results were summarized in Figs. [2~T2"1 
and[H31 

By comparing the phase diagram of massless quarks shown in the right panel of Fig. 12.111 
with the phase diagram shown in Fig. 12.121 I noticed several important differences. First of all, 
I observed that a self-consistent treatment of quark masses strongly influences the competition 
between different quark phases. As was noticed earlier in Ref. [139], there exists a subtle interplay 
between the two main effects. On the one hand, the actual values of the quark masses directly 
influence the competition between different normal and color-superconducting phases. On the 
other hand, competing phases themselves determine the magnitude of the masses. Very often, this 
leads to first-order phase transitions in which certain regions in the mass-parameter space become 
inaccessible. 

Some differences to the results in Sec. 12. H and Sec. 12.21 are related to a different choice of model 
parameters. Most importantly, the value of the diquark coupling Go = jGs is considerably 
weaker in the NJL model of Sec. 12.21 This can be easily seen by comparing the magnitude of the 
zero-temperature gap at a given value of the quark chemical potential, say at /i = 500 MeV, in the 
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two models. It is A 5UU) rs 140 MeV in Sec. O and A 5UU) « 76 MeV in Sec. O (Note that the 
strength of the diquark pairing in Ref. [112] is even weaker, corresponding to A 5 00 ^ « 20 MeV.) 

It should be noted that even the strong coupling Gd — Gs which leads to Ag 500 '' w 120 MeV, is 
still slightly weaker than that in Sec. 12.11 In this case, however, the corresponding results differ 
mostly because the quark masses are treated very differently. 

Because of the weaker diquark coupling strength, the Cooper instabilities in Fig. 12.121 occur 
systematically at higher values of the quark chemical potential than in the right panel of Fig. 12.111 
In particular, this is most clearly seen from the critical lines of the transition to the (g)CFL phase. 
Another consequence of the weaker interaction is the possibility of a thermal enhancement of 
the (g)2SC Cooper pairing at intermediate values of the quark chemical potential. This kind of 
enhancement was studied in detail in Refs. [84,85]. Making use of the same arguments, one can 
tell immediately how the phase diagram in Fig. 12.121 should change with increasing or decreasing 
the diquark coupling strength. 

In particular, with increasing (decreasing) the diquark coupling strength, the region of the 
(g)2SC phase at intermediate values of the quark chemical potential should expand (shrink) along 
the temperature direction. The regions covered by the other (i.e., uSC and CFL) phases should 
have qualitatively the same shape, but shift to lower (higher) values of the quark chemical potential 
and to higher (lower) values of the temperature. In the case of strong coupling, in particular, these 
general arguments are confirmed by my numerical calculations. The corresponding phase diagram 
is shown in Fig. 12.131 

Several comments are in order regarding the choice of the NJL model used in Sec. 12.21 The 
model is defined by the set of parameters in Eq. (|2.119[) which were fitted to reproduce several 
important QCD properties in vacuum [136]. (Note that the same model also was used in Ref. [133].) 
It is expected, therefore, that this is a reasonable effective model of QCD that captures the main 
features of both, chiral and color-superconducting pairing dynamics. Also, a relatively small value 
of the cutoff parameter in the model, sec Eq. (I2.119|) . should not necessarily be viewed as a bad 
feature of the model. In fact, this might simply mimic a natural property of the full theory in 
which the coupling strength of relevant interactions is quenched at large momenta. 

In this relation, note that the approach of Ref. [112] regarding the cutoff parameter in the NJL 
model is very different. It is said there that a large value of this parameter is beneficial in order to 
extract results which are insensitive to a specific choice of the cutoff. However, I do not find any 
physical argument that would support this requirement. Instead, I insist on having an effective 
model that describes reasonably well the QCD properties at zero quark chemical potential. I do 
not pretend, of course, that a naive extrapolation of the model to large densities can be rigorously 
justified. In absence of a better alternative, however, this seems to be the only sensible choice. 

The results of Sec. 12.21 might be relevant for understanding the physics of (hybrid) neutron 
stars with quark cores, in which the deleptonization is completed. In order to obtain a phase 
diagram that could be applied to protoneutron stars, one has to generalize the analysis to take 
into account neutrino trapping. 

In Sec. 12.31 1 studied the effect of neutrino trapping on the phase diagram of neutral three-flavor 
quark matter within the NJL model of Ref. [136]. The results were obtained in the mean-field 
approximation, treating constituent quark masses as dynamically generated quantities. The overall 
structure of the phase diagram in the space of three parameters, namely temperature T, quark 
chemical potential fi and lepton-number chemical potential fJ,L e , was summarized in Fig. 12.201 
This was further detailed in several two-dimensional slices of the phase diagram, including phase 
diagrams in the plane of temperature and quark chemical potential (see Figs. 12.211 and 12. 22j) and 
in the plane of temperature and lepton-number chemical potential (see Figs. I2.2TS1 and l2~2"Tl) . 

By making use of simple model-independent arguments, as well as detailed numerical calcu- 
lations in the framework of the NJL model, I found that neutrino trapping helps Cooper pairing 
in the 2SC phase and suppresses the CFL phase. In essence, this is the consequence of satisfying 
the electric neutrality constraint in the quark system. In two-flavor quark matter, the (positive) 
lepton-number chemical potential fiL c helps to provide extra electrons without inducing a large 
mismatch between the Fermi momenta of up and down quarks. With reducing the mismatch, of 
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course, Cooper pairing gets stronger. This is in sharp contrast to the situation in the CFL phase 
of quark matter, which is neutral in the absence of electrons. Additional electrons due to large 
fii, e can only put extra stress on the system. 

In application to protoneutron stars, my findings in Sec. 12.31 strongly suggests that the CFL 
phase is very unlikely to appear during the early stage of the stellar evolution before the delep- 
tonization is completed. If color superconductivity occurs there, the 2SC phase is the best candi- 
date for the ground state. In view of this finding, it might be quite natural to suggest that matter 
inside protoneutron stars contains little or no strangeness (just as the cores of the progenitor 
stars) during the early times of their evolution. In this connection, it is appropriate to recall that 
neutrino trapping also suppresses the appearance of strangeness in the form of hyperonic matter 
and kaon condensation [48]. My finding, therefore, is a special case of a generic property. 

The authors of Ref. [83] claimed that the 2SC phase is absent in compact stars and underpinned 
their claim by using a very simple model-independent calculation. The phase diagrams which were 
created by using more precise NJL-type models and which are shown in my thesis are the best 
evidence that there indeed could exist a 2SC phase in compact stars and that the 2SC phase in 
general is not absent. This is always the case in protoneutron stars and in cold neutron stars at 
large diquark coupling constants. But in the limit of weak diquark coupling and zero temperature 
in which the simple model-independent calculation in Ref. [83] is valid, the statement seems to be 
true that the 2SC phase is absent, cf. Fig. 12.121 

After the deleptonization occurs, it is possible that the ground state of matter at high density 
in the central region of neutron stars is the CFL phase. This phase contains a large number of 
strange quarks. Therefore, an abundant production of strangeness should happen right after the 
deleptonization in protoneutron stars. If realized in nature, in principle this scenario may have 
observational signatures. 

3.2 Open questions and outlook 

The quark sector of the phase diagram of strongly interacting matter was poorly understood in 
2003, see Fig. 11.11 With my investigations of this subject, I made the first attempt to illuminate 
the phase structure in the quark sector of the phase diagram of strongly interacting matter within 
the framework of an NJL model. In Fig. 13.11 I finally show the status of knowledge about the 
phase diagram of strongly interacting matter in 2005. 

With my investigations I took the first step in order to study the phase diagram of neutral 
quark matter, see Ref. [110] and Sec. 12.11 I further improved my result for the phase diagram of 
neutral quark matter by properly incorporating the strange quark mass [113]. The next step in my 
investigation was to treat the quark masses self-consistently as dynamically generated quantities, 
see Ref. [114] and Sec. \'2.'2\ and finally to incorporate a nonzero lepton-number chemical potential, 
see Ref. [116] and Sec.O 

Here, I restricted my study to spin-zero color-superconducting phases only which seem to be 
the most preferred color-superconducting states because they have the highest pressure. Therefore, 
I do not expect that other than these color-superconducting states which I have investigated in 
this thesis will occur in the phase diagram of neutral quark matter. But if the pressures of one 
or more of these other color-superconducting quark phases are larger than those of the spin-zero 
color-superconducting quark phases so that these other color-superconducting quark phases indeed 
will occur in the phase diagram of neutral quark matter, then I expect that the regions of such 
color-superconducting phases in the phase diagram are small so that the phase diagram in Fig. 12.121 
will not be changed very much. In order to avoid misunderstandings, I want to say that I do not 
expect that such color-superconducting states occur in the phase diagram of neutral quark matter. 
I do not say that such color-superconducting phases are absent. 

An example for a possible occurrence of a color-superconducting phase which I have not con- 
sidered in my investigations is the spin-one color-superconducting transverse CSL phase (the most 
preferred spin-one color-superconducting quark phase [39] ) . This phase has good chances to occur 
very close to the axis of the quark chemical potential in the region of the normal quark phase in 
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the phase diagram of neutral quark matter in Fig. 12.121 If one will find that this is really true then 
the phase diagram of neutral quark matter shown in Fig. 12.121 will not be changed much because 
spin-one color-superconducting quark phases like the transverse CSL phase will be destroyed even 
at small temperatures. Therefore, the line width of the axis of the quark chemical potential is 
possibly thicker than the region of the spin-one color-superconducting transverse CSL phase in the 
phase diagram of neutral quark matter. Nevertheless, the CSL phase could occur in cold neutron 
stars [38]. 
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Figure 3.1: The knowledge about the phase diagram of strongly interacting matter in 2005. The upper 
left part shows the critical endpoint of the quark-hadron phase transition which is obtained by 
lattice QCD calculations [5] . The lower right part shows the phase diagram of neutral three- 
flavor quark matter [114], where quark masses are treated self-consistently, see Fig. 12.121 
The dashed line shows the discrepancy between the critical endpoint of the quark-hadron 
phase transition in the NJL model with that obtained in lattice QCD calculations. This 
discrepancy could be resolved by using full QCD. In addition, there is shown the state of 
the early universe and of cold neutron stars after deleptonization. 



In Ref. [143], another possible spin-zero color-superconducting quark phase, the A-phase, was 
checked if it occurs in the phase diagram of neutral quark matter. The result was that this phase 
always has a lower pressure than the phases which I have investigated in my thesis. Therefore, 
the A-phase does not appear in the phase diagram of neutral quark matter. 

However, it will be very important to check if other than these color-superconducting phases 
which I have investigated in this thesis indeed occur in the phase diagram of neutral quark matter 
(see Ref. [144], where some of these other color-superconducting quark phases are shown). Such 
an investigation is outside the scope of my thesis. Much work concerning this check has to be 
done in future. The pressures of all competing phases have to be compared. The phase with the 
highest value of the pressure is the preferred one and therefore appears in the phase diagram of 
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neutral quark matter. 

Despite the progress in understanding the phase diagram of neutral quark matter, there still 
exists a fundamental problem. Gapless color-superconducting quark phases are unstable in some 
regions of the phase diagram of neutral quark matter because of chromomagnetic instabilities [92] 
so that another phase will be the preferred state. Chromomagnetic instabilities occur even in 
regular color-superconducting quark phases. The author of Ref. [93] shows that chromomagnetic 
instabilities occur only at low temperatures in neutral color-superconducting quark matter. The 
author of Ref. [94] points out that the instabilities might be caused by using BCS theory in 
mean-field approximation, where phase fluctuations have been neglected. With the increase of the 
mismatch between the Fermi surfaces of paired fermions, phase fluctuations play more and more an 
important role, and soften the superconductor. Strong phase fluctuations will eventually quantum- 
disorder the superconducting state, and turn the system into a phase-decoherent pseudogap state. 
By using an effective theory of the CFL state, the author of Ref. [95] demonstrates that the 
chromomagnetic instability is resolved by the formation of an inhomogeneous meson condensate. 
The authors of Ref. [96] describe a new phase in neutral two-flavor quark matter within the 
framework of the Ginzburg-Landau approach, in which gluonic degrees of freedom play a crucial 
role. They call it a gluonic phase. In this phase, gluon condensates cure a chromomagnetic 
instability in the 2SC solution and lead to spontaneous breakdown of the color gauge symmetry, 
the [C/(l) cm ] and the rotational SO(3) group. In other words, the gluonic phase describes an 
anisotropic medium in which color and electric superconductivity coexist. In Ref. [97], it was 
suggested that the chromomagnetic instability in gapless color-superconducting phases indicates 
the formation of the LOFF phase [98] which is discussed in Ref. [99] in the context of quark 
matter. Other possibilities could be the formation of a spin-one color-superconducting quark 
phase, a mixed phase, or a completely new state. The authors of Ref. [100] suggest that a mixed 
phase composed of the 2SC phase and the normal quark phase may be more favored if the surface 
tension is sufficiently small [101]. The authors of Ref. [102] suggest a single first-order phase 
transition between CFL and nuclear matter. Such a transition, in space, could take place either 
through a mixed phase region or at a single sharp interface with electron-free CFL and electron- 
rich nuclear matter in stable contact. The authors of Ref. [102] constructed a model for such an 
interface. 

Another possibility in order to avoid gapless color-superconducting quark phases is to set the 
diquark coupling constant to such high values that gapless phases will not occur in the phase 
diagram [118]. But then it is another question if such large diquark coupling constants indeed 
occur in nature. My personal opinion is that one should not increase the diquark coupling constant 
to such high values in order to get rid of the problem but use the standard value of the diquark 
coupling constant in Eq. (|2.120p which is predicted by the Fierz transformation in vacuum. 

Whether gapless phases can indeed exist as color-superconducting quark phases or whether 
they will be unstable because of chromomagnetic instabilities is not known yet. Also, if they are 
unstable, it is not yet known which states will be formed in the respective regions of the phase 
diagram. Such an investigation is outside the scope of my thesis, and much work has to be done in 
future to solve this problem. The pressures of all competing quark phases have to be determined 
in these regions of the phase diagram. The phase with the highest value of the pressure is the 
preferred one. 

Another problem which has to be solved is the discrepancy between the critical endpoint of 
the quark-hadron phase transition which is obtained by lattice QCD calculations [5] and NJL 
models [106, 114-116, 118, 138, 139], see also Fig. 13.11 It is appropriate to mention here that the 
location of the critical endpoint might be affected very much by fluctuations of the composite 
chiral fields. These are not included in the mean-field studies of the NJL model. In fact, this 
is probably the main reason for their inability to pin down the location of the critical endpoint 
consistent for example, with lattice QCD calculations [5]. It is fair to mention that the current 
lattice QCD calculations are not very reliable at nonzero /i either. Therefore, the predictions of 
my thesis, as well as of those in Refs. [106, 114-116, 118, 138, 139], regarding the critical endpoint 
cannot be considered as very reliable. One definitely will be able to say where exactly the critical 
endpoint is located and how the phase diagram of neutral quark matter looks like if one uses 
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full QCD. Such an investigation is outside the scope of my thesis and has to be done in future. 
Since the NJL-type models used in my thesis fulfill at least some important key features of QCD, I 
expect that the phase diagram of neutral quark matter which is shown in Fig. 12.121 will be changed 
only quantitatively so that the overall phase structure does not change if one computes a more 
realistic phase diagram of neutral quark matter by using full QCD. 

In the context of compact stars, it is not known yet if their cores reach chemical potentials which 
make color superconductivity possible. It is of great importance to resolve this issue. Neutron stars 
are the most likely candidates for color superconductivity in nature and therefore, one expects that 
they contain color-superconducting quark cores. Therefore, one has to find an experimental proof 
for color-superconducting quark matter. This can be done by detecting neutrino and 7-emissions 
of protoneutron stars and in colliders. 

Since one does not know yet how large the lepton-number chemical potential in protoneutron 
stars is during their evolution from the supernova explosion until they become a cold neutron star, 
it is difficult to say exactly, through which phases they go in the three-dimensional phase diagram 
shown in Fig. 12.201 Therefore, it would be very interesting to study the precise evolution path of 
protoneutron stars in the three-dimensional phase diagram. 

I note that the analysis in this thesis is restricted to locally neutral phases only. This automat- 
ically excludes, for example, crystalline [99] and mixed [100, 101, 120] phases. Also, in the mean 
field approximation utilized here, I cannot get any phases with meson condensates [28-30]. 

In nature, phase transitions can take place either through sharp boundaries between pure 
phases which are located next to each other or through mixed-phase regions. In this thesis, I 
considered only the former possibility, where each phase has the highest pressure in its region 
of the phase diagram and where the pressures of the phases are equal only at the sharp phase 
boundaries. Such phase transitions are obtained by so-called Maxwell constructions. It would be 
interesting to study the possibility of mixed phases by making so-called Gibbs constructions [145], 
in which the pressures and the chemical potentials of the phases are equal and continuous in 
contrast to Maxwell constructions. 

Another interesting investigation in the context of compact stars would be to create a phase- 
radius diagram that shows which phase occurs at which radius of the compact star. Such a diagram 
would directly show the composition and structure of a compact star. 
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Appendix A 

Definitions of matrices 



In this chapter of the Appendix, I present some definitions such as the Pauli and 7-matrices, and 
the generatores of the 51/ (3) group. I also show the chirality, energy, energy-chirality, and the 
spin projectors as well as their properties. 

A.l The Pauli matrices 

The Pauli matrices are defined as: 



n f -i 



- 1 1 0) ' n " v. - 
1 1 \ 



(A.1) 



The Pauli matrices arc hcrmitian, 



oj = a t . (A.2) 
The following relations are valid: 

of = 1 , Oi&j = in k , [at, Uj] = 2ia k , {o-j, Uj} = 2S lj . (A. 3) 

A. 1.1 Spin projectors 

The spin projectors are defined as 

V s (k) = i (1 + sa ■ k) , (A.4) 

where s = ± corresponds to the projections onto spin up or spin down states, respectively, and a 
are the Pauli matrices. The spin projectors fulfill the properties of completeness and orthogonality, 

V+ (k) + V- (k) = 1 , (A.5a) 

P a (k)7V(k) - ^Py(k) . (A.5b) 

The spin projectors are hermitian, 

7>+(k) = 7> a (k) . (A.6) 



The following relations are valid: 



(k) - P_ (k) = a ■ k , (A.7) 
TrP s (k) = l. (A.8) 
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A. 2 Matrices in Dirac space 

The Dirac definition of the 7-matrices reads: 



_ / 1 0\ ! ( tri 



I .. 1 j • 7 -v- ffx 

2 _ 1 o- 2 \ 3 _ ( a 3 



The following relations are valid: 

7 M = (7°, 7 1 , 7 2 , 7 3 ) = (7°, 7) , 70 (7 M ) f 70 = 7 M , {l^ l"} = W 



(A.10) 



where 

<T = diag (1,-1, -1,-1) (A.ll) 
is the metric tensor. The charge-conjugation matrix in Dirac representation is defined as 

C = z 7 2 7o , (A.12) 

where 7 2 is given by Eq. (|A.9j) . The following relations are valid: 

C = -C- 1 = -C T = -C f , (A.13a) 
C^C- 1 = - ( 7 ") T , (A.13b) 
Ci^fc- 1 = -7" . (A.13c) 

With the charge-conjugation matrix, one can introduce the charge-conjugate spinors in space-time, 

i, c (A) = (X) , 
i>c (A) = V T (A) C , 

t/> (A) = CVc (X) , 

$ (A) = 1% (A) C , 



(A.14) 



and in energy-momentum space, 



4>c (K) = CV (-K) , 
4>c (K) = V T (-K) C , 
V (- A) = CV C (A) , 
i> (- A) = (K) C . 



(A.15) 



A. 2.1 Projectors in Dirac space 

In the limit of vanishing mass, the chirality, energy, and energy-chirality projectors are defined as: 

^0=3(1 + 075), (A.16a) 
A e (k) = i(l + e 70 7-k) , (A.16b) 
VI (k) = V c A e (k) = |(1 + c 75 ) (1 + e 7 o7 • k) , (A.16c) 



A. 3 The generators of the SU(3) group 
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where c = ± or c = r, £, respectively, stands for right-handed and left-handed projections, and 
e = ± for projections onto states of positive and negative energy. The chirality and energy 
projectors fulfill the properties of completeness and orthogonality, 

V r +Vi = l, (A. 17a) 

A+(k) + A-(k) = 1, (A.17b) 

V C V C > = 8 c . d V d , (A. 17c) 

A e (k) A e '(k) = 6 e < e 'A e '(k) . (A.17d) 

These projectors are hermitian, 

V\ = T c , (A. 18a) 

[A<f (k) = A e (k) , (A.18b) 

and they commute with each other, 

[7> c ,A e (k)]=0. (A.19) 
The energy-chirality projectors fulfill the properties of completeness and orthogonality, 

J2K(k) = l, (A.20a) 

c,e 

V e c (k) Vi (k) = 5 c , c ,5^'vi (k) . (A.20b) 
The energy-chirality projectors are hermitian, 

[P c f (k)=P c e (k) . (A.21) 

The following relations are valid: 

PZ c e (k) = 7o P c e (k)7o, (A.22) 

A+ (k) - A- (k) = 7o7 • k , (A.23) 

Tr-P c e (k) = 1 , (A.24) 

Vcl n = ^ v _ c (A _25) 



A. 3 The generators of the SU(3) group 

The generators of the SU(3) group are defined as: 



Ti = — 



n = - 



T7 = — 








-1 


-1 
000 

1 

0' 
-i 
i 



£2 = 77 



r 6 = x 




(A.26) 
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The Gell-Mann matrices are defined by 

A a = 2T a . (A.27) 

The following relations are valid: 

8 

TrT a = 0, Tl = T a , [T a ,T b ] = ^^ f abc T c . (A.28) 

c=l 

All fabc = 0, except, 

/i23 = 1 , (A.29a) 

/l47 = — /l56 = /246 = ./W = /345 = ~ /367 = | , (A. 29b) 

/458 - /678 - . (A.29C) 



Appendix B 

Useful formulae 



Several relations which are used in this thesis are derived and proven in this chapter of the 
Appendix. 



B.l Non-interacting massless fermions and antifermions at 
nonzero temperature 

The pressure of non-interacting massless fermions and antifermions at nonzero temperature T is 
given by 



P - 



2tt 2 



J dkk 2 | In 1 + cxp ^ 



k — n 



+ ln 



1 + exp 



k + fi 



(B.l) 



where g is the degeneracy factor, k = |k| is the momentum, and /i is the chemical potential of the 
fermions. The first term corresponds to the pressure of fermions while the second one corresponds 
to the pressure of antifermions. By an integration by parts one obtains, 



P 



9 1 dkk 3 



6tt 2 



no- 



where 



n F (x) 



k — [i 



1 



n F 



k + n 
T 



e x + 1 

is the Fermi-Dirac distribution function. 

With the substitutions x ± = (k± fi)/T one computes, 



(B.2) 
(B.3) 



6tt 2 



d X - (Tx~ + n) 3 n F (or) + / da; + (Tx+ - n) 3 n F (x+) 
-iL Jr. 



(B.4) 



The integrals can be split so that 



P = 



6tt 2 



/ dx- (Tx~ + fj) 3 n F (x~) + / dx+ (Tx+ - /i) 3 n F (x+) 
Jo Jo 

/ da;" (Tx~ + fi) n F (x~) + / dx+ (Tx+ - fi) n F (x + 



(B.5) 



The first two integrals can be joined. In the last integral, one can use the substitution x + = —x + . 
Then, also the last two integrals can be joined, 



P = 



6tt 2 



dx 



2T i x d + 6Txfi 2 
e* + 1 



J ^ dx (Tx + fif [n F (x)+n F (-x)\ | , (B.6) 



114 



Useful formulae 



where I renamed x + , x + , and x in x. With the relation np (x) + np (—x) = 1, the integrals [130, 
131], 

00 dx -^t = -U 4 , /'°° dx = ±-ir 2 , (B.7) 



o 



c x + 1 120 7 e* + 1 12 



and the substitution, z — Tx + fi, one obtains for the pressure of non-interacting massless fermions 
and antifcrmions at nonzero temperature [80], 

9 I .A , r>_2,,2r r 2 , ' „4rri4 



B.2 The inverse Dirac propagator 

The inverse tree-level propagator for quarks and charge-conjugate quarks, respectively, is given by 
[G+f 1 (X, Y) = {iJ3 x + A7o - rh) 5^ (X - Y) , (B.9a) 

[$„ ] _1 (x, r) = {%ja% - Ato - m) <5 (4) (x - y) , (B.9b) 

where — — igA^T a and = + igA a ^Tj are the covariant derivatives for quarks and 
charge-conjugate quarks, respectively, jl is the color-flavor matrix of the chemical potentials, and 
rh is the quark-mass matrix. By using the charge-conjugate spinors (|A.14I) . the relation 

$(X) [g+]- 1 {X,Y)^{Y)=^ c {Y) [G ]- 1 (Y,X)^c(X) (B.10) 

can be proven: 

^(X) [Qt\- l {X,Y)^{Y) = ^(X)C[G+}- 1 (X,Y)C^(Y) 

= -rPl (X) C [%J3 X + A7o - m] (X - Y) C' 1 ^ (Y) 

= -^l (X) [iC^C^D^ + fiCjoC- 1 ~ rh] (X - Y) ^ (Y) 

= (X) H (rf D, - A 7 o - m] 6^ (X - Y) ^ (Y) 

= (Y) { [-i (r f - A7o ~ m]} T 6^ (Y - X) yj c (X) 

= ^ c (Y) [iJ3 c x - A 7 o - rh] S (4) (Y - X) ^ c (X) 

= (Y) [g^r 1 (Y,X)tPc (X) . (B.ll) 



B.3 The tree-level quark propagator 

The inverse tree-level quark propagator in Nambu-Gorkov space is given by 

S [ $+ [Qo^ + Jl-) ' (R12) 

where [C^j 1 ] -1 is the inverse tree-level propagator for quarks or charge-conjugate quarks, respec- 
tively, are the regular quark self-energies, and <I )± are the gap matrices. From this expression, 
the tree-level quark propagator in Nambu-Gorkov space, 

S^(t+r~-), (B.13) 



(B.14) 



e+ g 

can be obtained by using the relation SS^ 1 = = 1, 

g+ e- \ /[g+]- 1 + s+ a- o 
s+ g-J { $+ [e-j-i + ir J i 



B.4 The Feynman gauged gluon propagator 
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where is the propagator for quasiquarks or charge- conjugate quasiquarks, respectively, and "E^ 1 
are the anomalous propagators. This can also be written in the form, 

Q + (Kir 1 + S + ) + = i , ( B .15a) 

g+<S>- + 3- ([Go]- 1 + £") = , (B.15b) 

E+ ([G+]- 1 +^+)+g-^+ = 0, (B.15c) 

S+$-+^-([g -]- 1 + E-) = 1. (B.15d) 

From Eqs. (|B.15bjl and (|B.15c|l one gets 

S ± --^ (& ±([g±]-i +I] ±)^ 1 . (B.16) 

By inserting these relations into Eqs. (|B.15a[) and (|B.15d[) . one obtains 

q ± = {[otr 1 + s ± - $ t ([Qf]- 1 + ^r 1 ^y 1 ■ (B.17) 

In the same way, one proceeds with 

[^o + ]- 1 + s + *- \(g+E-\_(io\ (R18) 



and gets 

S ± = -([^]- 1 +S^)" 1 $ ± ^ ± . (B.19) 

By inserting these results into Eq. (|B.18[) . one again obtains for Q + and Q~ the relations (|B.17[) . 
For E , one obtains 

S ± = -([^]- 1 + ST)- 1 $±g± = -^ <i) ±([g± r i + E ±)- 1 . (B .20) 

B.4 The Feynman gauged gluon propagator 

The QCD Lagrangian density is given by 



C = i, {iJ3 - m) 4, - -F£ v Fr + £„ ■ (B.21) 



In this section, I focus on the contribution of gluons to the QCD Lagrangian density, 



C A ee --F^ v Fr ■ (B.22) 



The gluon field-strength tensor is defined as 

Fr = d»Al - d v A» + gf abc A%A» . (B.23) 
By neglecting the gluon self-interaction, the gluon field-strength tensor simplifies to 

Fr = d»A» a - d v A» , (B.24) 
which breaks [ST/(3) C ] gauge invariance. The Abelian part of the gluon Lagrangian density is 

£a = -\f^F^ = -~ (d,Al d u A«) {&>A" a - d»A») 

= -- (d,Al d»A» a - d,A a v <TA» - duAl flM£ + d„A; d»A») 



--(d^d^A^-d^d'A^) . (B.25) 
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After integration by parts, and discarding surface terms, the Abelian part of the gluon Lagrangian 
density reads [146], 

C A = \ {Ald^Al - AldpFA*) = \A» (g^D d,d u ) A v a . (B.26) 
By choosing Lorentz gauge, d^A^ = 0, one obtains, 

C A = l -A^ v UA v a . (B.27) 
One gets the same result if one sets the gauge-fixing term of the QCD Lagrangian density to 

£ g au gc = ~\ (dM 2 , (B.28) 

which leads to the QCD Lagrangian density in Feynman gauge. The Feynman gauged gluon 
Lagrangian density can be rewritten as 

C a = \a^\D-X,K- (B.29) 

The quantity 

[D-Xl = <5 a Vn (B.30) 
is called the inverse Feynman gauged gluon propagator. Inverting yields 

Df v = 5 ab 9llv U- 1 . (B.31) 

In momentum space, the Feynman gauged gluon propagator reads 

D ab = s ab ^ . (B.32) 

B.5 The determinant of the inverse quark propagator 

In space-time, the contribution of the kinetic term to the grand partition function reads up to 
irrelevant constants [7, 128] 

Zkin = J V^ViP exp {/ kin [*, tf] } , (B.33) 
where the action of the kinetic term is 

4m [*,*]=£/ * (X) S- 1 (X, Y) * (Y) . (B.34) 

Z JX,Y 

The inverse quark propagator (without regular self-energies) 



(B.35) 



is a 2 x 2-matrix in Nambu-Gorkov space, where 

*=( 

ipc 



*= (V^c) , *= (. t ) , (B-36) 
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are the Nambu-Gorkov quark spinors. The Fourier transforms of the quark spinors read 



ib ( X\ — 


v v K 


l U . Ota) 


${X) = 




(B.37b) 


1>c (X) = 


V V K 


(B.37c) 


4>c (X) = 


^E^c(A-). 


(B.37d) 



By assuming translational invariance, the Fourier transform of the inverse quark propagator is 
given by 

S- 1 (X, Y) = ^ e- lK{x - Y) S- x (K) . (B.38) 



K 



By Fourier transforming the kinetic term of the action, the integration measure has to be rewritten, 
V$Vip = Y[ dtp {K) dip (K) 



K 



= JJ d{p{K)d^(-K)dip{K)dip{-K) 
(K-K) 

= Af Yl d ^ (K) d V>c (K) d^P (K) dtpc (K) 



(B.39) 



(K-K) 



where AT is an irrelevant constant Jacobian which arises from charge conjugation of quark spinors. 
The Fourier-transformed kinetic term of the action is 



/kin *, * 



T 

T 
W 1 



J2 f ^(K)S- 1 (Q)^(P)^ KX e- lPY e-^ x - Y ^ 

C,P,Q Jx ^ Y 

V f S K -Q) X [ e i ^- p ^ Y ^(K)S- 1 (Q)^(P) 



= E ^ { k] Q 1^{K)S-\Q)^{P) 

K,P,C J 



K ' (K,-K) 



(K-K) 



S-'jK) 
T 



V(K) . (B.40) 



Therefore, the contribution of the kinetic term to the grand partition function in energy-momentum 
space reads [146] 



Zki„ = J Xtytftyexp {/ kin [*,*]} = ^ { det^ f ^_ j = 



det 

k \ T 



1/2 



(B.41) 



Up to irrelevant constants, the logarithm of the kinetic part of the grand partition function is 



lnZ, 



kin 



In 



det 

K 



s- 1 

T 



1/2 



- In det , 

2 k \ T 



2 ln 



11** 



K 



1 ^ In det (— ) = ~Y 

2 4^ \ T J 2 4^ 



Tr In 



s- 1 

T 



(B.42) 



K x 7 if 

The relation in the last of line of this equation is proven in Sec. IB. 61 
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B.6 The logarithm of the determinant 

The following relation is to be proven: 

IndetAf = TrlnM . (B.43) 

The quadratic matrix M with dimension d is diagonalizable if a matrix S exists, which fulfills the 
property [129], 

SMS' 1 = diagAi , (B.44) 
where Xi are the eigenvalues of the matrix M. Therefore, one can write, 

d d 

IndetAf = lndet(SAf S^ 1 ) = In ]J A, = ^ In A, =TrlnA, . (B.45) 

»=i *=i 

The logarithm can be expanded as 

°° (— 1 1 1 

ln(l + x) = Y - x n = x - -x 2 + -x 3 --x 4 ±... , (B.46) 



n=l 



so that Eq. (|B.45|) becomes, 



IndetM = V L-^- Trdiag{(A» - 1)"} 



n-1 



n 

n—l 



j2 t ^-^{[s(M-i ) s-r} 



n 

n=l 



= — — Tr [s ( M - l T ■ ( R47 ) 

n=l n 

Since Tr (AS) = Tr (IL4.) and because of Eq. (|B.46[) . one obtains that 

In det M = V ^— ' Tr (M - 1)" = Tr In M . (B.48) 



n-1 



n—l 



B.7 The Dirac trace 

By using the Minkowski space part of the gluon propagator (|B.32I) , 



one obtains 



= , (B.49) 



Tr [V c A e (k) -f&T* (P) l v ] = ± {-Tr [V c A e (k) A e ' (p)] + Tr [V c A e (k) f A" e ' (p) f ]} . 

(B.50) 

Expanding the projectors yields, 

V c A e (k) A ±e ' (p) = | (1 ± e'707 • p + e^-y ■ k qp ee'7 • k • 7 • p 

+ C75 ± ce'75707 • p + ce757o7 ■ k =p cee'757 ■ k ■ 7 • p) . (B.51) 

With the relation, 

- 7 • k • 7 • p = - 7 ! 7 4y = -| ({ 7 \ 7 J } + [ 7 \ 7^]) tip 1 = k • p - i [ 7 \ 7 J ] fc'pj , (B.52) 



B.8 The trace of the logarithm 
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one computes the trace of the product of these projectors, 

Tr [V c A e (k) A ±e ' (p)] = | (1 ± ee'k • p) . (B.53) 
Herewith, the first trace in Eq. (|B.50p is solved. From the second trace in Eq. (|B.50| . one gets, 
Tr [P c A e (k) f A- e ' (p) 7 *] = iTr [P c A e (k) (7* - e y 7 o7 ■ P)f} 

= iTr [V c A e (k) (f + e' 7 o7Vp J )7l , (B.54) 

which can be rewritten as 

iTr [P c A e (k) ( 7 l + e'7o{-7 J 7' i ~ ^ijWh'] 
= iTr [V c A e (k) ( 7 V - e' 70 7 • P7Y - 2e' 7o %pV)] . (B.55) 
Summing over all i yields, 
Tr[P c A e (k) (-3A- e ' ( p )- e ' 7o7 .p)] = —3 Tr [P c A e (k) A~ e ' (p)] -Tr [P c A e (k) e' 7o7 .p] ■ (B.56) 

The first of these two traces is already known by Eq. (|B.53|) . The second one is determined by 
using the relation (|B.52|) . 

Tr [V c A e (k) e' 70 7 ■ p] 
= 3 Tr [(1 + c 75) (1 + e7o7 • k)e' 7o 7 • p] 
= |Tr [(1 + e7o7 • k + 075 + ce7 5 7 o7 • k)e'7o7 • p] 
= iTr [(e'707 ■ p - ee'7 • k ■ 7 • p + ce' 757o 7 • p - cee'j 5 ~f • k • 7 ■ p)] 



= ee'k • p . 

Herewith, one obtains the second trace in Eq. (|B.50|) . 

Tr [P c A e (k) 7 4 A- e ' (p) 7 l ] = -i (3 - ee'k • p) . 
Therefore, one calculates for the Dirac trace (IB.50j) . 

2 

A 1 



Tr [V c A e (k) 7 ^A" e (p) 7 "] (K - P) = - — 



(B.57) 
(B.58) 

(B.59) 
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The relation, 

tci m 

Tr In Y a iPi = Y ln a * Tr Vl ' ( R6 °) 

i=l i=i 

is to be proven. Each trace in this section is taken over a d x d matrix. The Vi are projectors, 
which fulfill the properties of completeness, Yli^i = 1) an< ^ orthogonality, V{Pj — SijVj, and the 
di are factors, where 1 < i < m and 1 < j < m. Because of the completeness of the projectors, 
one gets, 



Tr ln Y citVi = Tr ln axVx + Y a{P l = Tr ln 

»=1 V i=2 / 

which can be rewritten as 

m 

Trln^a^i = Trln 



01 1 Y v - ■ Y" ;r - 



i=2 



i=2 



(B.61) 



i=i 



a 1 



Tr ln < ai 



Y ( a i - a l)^i 



i=2 



i=2 



= Trlnai + Trln 



i=2 v 1 7 



(B.62) 
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With the expansion of the logarithm, 



00 (— •n n_1 



Eq. (|BT62|) reads, 



Tr In ^ a-iVi = d In a x + Tr ^ ^— ^ 



«i 



n=l i=2 v ' 

00 / , n _i m , 

= h a f + V^Vp-:i | Tr A 

This equation can be rewritten as 

m m 

Tr In ^ a l V l = In af + J^Tr V t . In — 

i=l i=2 
m 

= lnaf + ^ln 

i=1 

Further calculations yield 

m / m 

Tr In £ ai Vi = In ^ JJ * 



ai 



hi 



»=2 v 



Tr "Pi 



(B.63) 



(B.64) 



(B.65) 



= in LK Ti ^^n«"^ 



, / d -Tr (1-7*0 TT TrPi 



i=2 



Finally, one obtains, 



Tr In £ 0,7^ = In Pl J] a? Pi = In ]J V ' = ]T ln «* Tr ?i 



(B.66) 



(B.67) 



B.9 Cubic equations 

In this section, I show Cardano's formulae which are the roots of the general cubic equation, 

A 3 x 3 + A 2 x 2 + A lX + A = , (B.68) 

where A 3 f [132,147,148]. 

Recall that a polynomial of degree n has n roots. The roots can be real or complex, and they 
might not be distinct. If the coefficients of the polynomial are real, then complex roots will occur 
in pairs that are complex conjugates, i.e., if xi = a + ib is a root then X2 — a — ib will also be a 
root. When the coefficients are complex, the complex roots need not be related. 

The solution to the cubic (as well as the quartic) equation was published by Gerolamo Cardano 
(1501-1576) in his treatise Ars Magna. However, Cardano was not the original discoverer of either 
of these results. The hint for the cubic equation had been provided by Niccolo Tartaglia, while the 
quartic equation had been solved by Ludovico Ferrari. However, Tartaglia himself had probably 
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caught wind of the solution from another source. The solution was apparently first arrived at 
by a little-remembered professor of mathematics at the University of Bologna by the name of 
Scipione del Ferro (ca. 1465-1526). While del Ferro did not publish his solution, he disclosed it to 
his student Antonio Maria Fior [149] . This is apparently where Tartaglia learned of the solution 
around 1541. 

The general cubic equation always can be brought into the normal form, 

x 3 + a 2 x 2 + a x x + a = , (B.69) 

where 

a 2 = -r- , ai = — , a = — . (B.70) 

A 3 A 3 A 3 

With the substitution x = y + s, the cubic equation in normal form assumes the form, 

y 3 + (3s + a 2 ) y 2 + (3s 2 + 2a 2 s + ai) y + s 3 + a 2 s 2 + a lS + a = . (B.71) 
By setting s = —a 2 /3 the quadratic term disappears. Thereby, one obtains with the substitution, 

*=y-?, (B.72) 



the cubic equation in standard form, 



where 



3 



v A + qy + r = , (B.73) 



q = a x - ^aj. , (B.74a) 
2 1 

r = —al - -a 2 ai + a . (B.74b) 
The cubic equation in standard form can be rewritten as 

y 3 + 3Qy - 2R = , (B.75) 

where 

Q = |, R = ~- (B-76) 

Let B and C be, for the moment, arbitrary constants. The cubic equation in standard form (|B.75jl 
can be decomposed into a linear and a quadratic term, 

(y - B) (y 2 + By + C)=y 3 + y(C-B 2 )-BC = 0, (B.77) 

so that one easiliy recognizes by comparing the coefficients that 

3Q = C-B 2 , (B.78) 

2R = BC . (B.79) 

Solving Eq. dBTSjl for C yields, 

C = B 2 + 3Q. (B.80) 
Inserting this result into Eq. (|B. 79|) gives, 

B 3 + 3QB = 2R . (B.81) 

Therefore, if one can find an expression for B which satisfies this identity, one has factored a 
linear term from the cubic equation, thus reducing it to a quadratic equation. The trial solution 
accomplishing this miracle turns out to be the symmetrical expression, 

B = S+ + S~ , (B.82) 
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where 

S ± =\fR±VD. (B.83) 
The discriminant is given by 

D = R 2 + Q 3 . (B.84) 

Now, the remaining quadratic equation has to be solved. Inserting Eq. (|B.80|) into the quadratic 
part of Eq. (|B.77P leads to 

y 2 + By + B 2 + 3Q = . (B.85) 
Solving this quadratic equation gives the solutions, 

B 1 



X\ = 


a 2 
3 


+ B, 


X2 = 


a-2 
3 


B 
~ ~2 


%3 = 


3 


B 
~ ~2 



y = --±-V3WB 2 +4Q . (B.86) 
With the definition, 

A = S + - S~ , (B.87) 

one can find out that 

A 2 = B 2 +AQ . (B.88) 

Therefore, with the substitution (|B.72[) . the roots of the general cubic equation (|B.68|) . which are 
called Cardano's formulae, are given by 

(B.89a) 

^iV3A , (B.89b) 

-iVZA . (B.89c) 

In the case of real- valued coefficients in the general cubic equation (|B.68|I , one can distinguish the 
following three cases: 

D > 0: one real and two complex solutions, 

ii el, £2,2:3 G C , x 3 = x* 2 , (B.90) 

D — 0: three real solutions, 

R =4 0: one single and one double solution, 

xi = -y + 2^R, x 2 = x 3 = ~-VR, (B.91) 
R — 0: one triple solution, 

xx = x 2 = £3 = -y , (B.92) 
D < 0: casus irreducibilis, three different real solutions, 

« = 2v^Qcos ( ^ + 2( ;~ 1)7r ) - I , (B.93) 

where 

R 

cos<^=^^=. (B.94) 

The equations for D < in which i = 1, 2, 3 first appear in Chapter VI of Frangois Viete's treatise 
De emendatione, published in 1615. 
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B.10 Summation over the fermionic Matsubara frequencies 

The following kinetic term of a fermionic logarithmic grand partition function is given: 



In Z kin = y^ln 
which can be split into two parts, 



K 



(fc - <S,j) 2 - e 2 



in z kin = 



In 



fco — Sfi + e 
f 



rp 2 



hi 



fco — Sfj, — e 
f 



(B.95) 



(B. 



This is up to irrelevant constants, 



Ax 



where 



One uses the relation, 



e ± 5n 



T^g* ^o) = ^~ /dfc tanh (^) g± (fc ) 

n \ / 



where 



9 ± (ko) 



ko ± Tx 



and integrates over the fermionic poles, 

fco = —i (2n + 1) 7rT = — iu> n 



(B.97) 



(B.98) 



(B.99) 



(B.100) 



(B.101) 



on the imaginary axis. The ui n are called fermionic Matsubara frequencies. The factor e vx in 
Eq. (|B.100j) is needed as a damping factor so that the sum over all n does not diverge, rj > 0. In 
the end of the calculation, one sets 7/ = 0. Because of the symmetry around the imaginary axis, 
see left panel in Fig. IB. 11 one obtains, 



T Ef ± W = 4^| +s dfcotanh^J [g± (k ) + g± (-ko)] . (B.102) 



With the relation, 



where 



tanh ^— J = 1 — 2np (x) 
n F [x) = r 1 - 



is the Fermi-Dirac distribution function, one gets 



dfc 



fc 



ko ± Tx —fco ± Tx 



T^g* (fco) = — 

~ 2tti j_ ioa+s 

In this equation, one uses the theorem of residues of first order, 

2^ //(*) = E J™ [/ {z) (z ~ z * )] = E Res / (*) 



(B.103) 
(B.104) 

(B.105) 
(B.106) 
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Figure B.l: Integration over the fermionic poles. Left panel: integration over the fermionic poles on the 
imaginary axis. Right panel: integration over the fermionic pole on the positive real axis. 



In this connection, one has to note that contour integrations are done anticlockwise. But the 
contour integration in Eq. (|B. 105|1 has to be done clockwise, see right panel in Fig. IB. II Therefore, 
the result has to be multiplied by minus one. One also has to note that the respective pole is on 
the positive real axis. Therefore, an extra Heaviside function has to be included in the result, 



Res5 ± (fc ) 



\k =Tx 



= -6(Tx) 



1 



fco ± Tx 
1 



— Up (x) 



(fco - Tx) + 



-fco ± Tx 



(fc - Tx) 



±9{Tx 



Up (x) 



± 



k =Tx 

- tanh (^)9 (Tx) , (B.107a) 



Res^fco)!^^ = -6{-Tx) 



1 



fco ± Tx 



— up (—x) 



(fco + Tx) + 



fco ± Tx 



(fco + Tx) 



k = -Tx 



7J->0 



TO(-Tx) 



1 



np (—x) 



= ± 



itanh(!)0(-Tz) 



T^tanh (-^6(-Tx) 



Summing up the two residua for each sign in Eq. (|B.105[) leads to the result, 

1 



T'£g ± (fc ) = ± 

n 

Therefore, one obtains for Eq. (|B.97|) 



np (x) 



(B.107b) 



(B.108) 



InZ, 



kin 



E 



dx 



2~ np ( x ) 



d:r 



- - n F (x) 



(B.109) 
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Up to irrelevant constants, one computes, 
In Z kin = ^2 | -f + ln 1 + exp ^ 



e — 5fi 



+ ln 



1 + cxp 



e + Sfi 



(B.110) 
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Anhang C 

Zusammenfassung 



Das Phasendiagramm von neutraler Quarkmaterie war, bis ich mit meinen Forschungen auf die- 
sem Gebiet im Jahr 2003 begonnen hatte, nur unzureichend bekannt, siehe Abb. 11.11 Fiir mcinc 
Dissertation war dieses Thema daher besonders gut geeignet. Meine Aufgabe war es, das Phasen- 
diagramm von neutraler Quarkmaterie zu erstellen. 

Im Kapitel [T] wird man in das Themengebiet Quarkmaterie, Farbsupraleitung und Sterne ein- 
gefuhrt. Neutronensterne sind in der Natur die wohl wahrscheinlichsten Platze, in denen farb- 
supraleitende Quarkmaterie vorkommen konnte. Daher gehe ich in Kapitel [T] nicht nur auf die 
wichtigstcn Phasen der Farbsupraleitung, sondern auch auf das Thema Sternentwicklung ein, um 
dann schlicfilich die Entstehung von Neutronensternen zu erklaren, deren Kerne moglicherweise aus 
sehr dichter, neutraler, farbsupraleitender Quarkmaterie bestehen. Bei extrem hohen Dichten, wie 
sie wahrscheinlich in Ncutronensternkernen vorkommen, herrscht eine attraktive Wechselwirkung 
in Quarkmaterie vor, die die Bildung von Quark-Cooper-Paaren bewirkt [10-12]. Dies verursacht 
Farbsupraleitung. Da Quarks als Spin-^ Fermionen in verschiedener Farbe und Flavor vorkommen, 
gibt es viele Moglichkeiten, um Quark-Cooper-Paare zu bilden. Daher gibt es viele unterschiedliche 
farbsupraleitende Quarkphasen. 

In der 2-flavor color-superconducting (2SC)-Phase paaren sich rote up mit griinen down Quarks 
und rote down mit griinen up Quarks und bilden antiblaue Quark-Cooper-Paare. Die blauen up 
und down Quarks bleiben ungepaart und sind daher gaplose Quasiteilchen, die einen grofien Beitrag 
zur spezifischen Warme, zur elektrischen Leitung und zur Warmeleitung leisten. Die gaplosen blau- 
en Quarks sind aufierdem fiir eine starke Neutrinoemission verantwortlich, die durch /3-Prozesse 
verursacht wird. 

Die Beitrage der vier gegappten Quasiteilchen zu alien Transport- und vielen thermodynami- 
schen Grofien ist bei niedrigen Temperaturen durch den exponentiell kleinen Faktor exp (— A/T) 
unterdruckt [9]. Die Gluonen sind Bosonen, und daher ist deren Teilchenzahldichte bei niedrigen 
Temperaturen klein. In der 2SC-Phase ist die [5f (3) c ] Eichsymmetrie zu [SU(2) C ] gebrochen. 
(Eichsymmetrien sind durch eckige Klammern gekennzeichnet.) Daher gibt es in der 2SC-Phase 
8 — 3 = 5 gebrochene Generatoren. Diese fuhrcn zu fiinf massiven Gluonen in der 2SC-Phase. 
Daher haben diese Gluonen nur geringe Einfliisse auf die Eigenschaften der Quarkmaterie in der 
2SC-Phase. Die ungepaartcn blauen Quarks sind fiir die Abwcscnhcit von Baryonsuperfluiditat 
verantwortlich. Nur die antiblauen Quasiteilchen tragen cine von Null verschiedene Baryonenzahl. 
In der 2SC-Phase existiert kein elektromagnetischer Meifinereffekt, und das ist der Grund dafiir, 
warum ein Magnctfcld aus der farbsupraleitenden Region nicht herausgedrangt werden wiirde. Die 
2SC-Phase ist ein sogenannter Q-Leiter, weil ihre elektrische Leitfahigkeit aufgrund der ungepaar- 
ten blauen up Quarks grofi ist. 

Wenn das chemische Potential der strange Quarks die strange Quarkmasse iibersteigt, dann 
treten auch die strange Quarks in der Quarkmaterie bei der Temperatur T = auf. Dadurch 
ist es moglich, dafi auch die strange Quarks bei der Bildung von Quark-Cooper-Paaren beteiligt 
sind. Das Kondensat bricht [SU(3) C ] x SU(3) r +e zur vektoriellen Untergruppe SU(3) c + r +e und 
ist noch immer invariant unter Vektortransformationen im Farb- und Flavorraum. Das bedeutet, 
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daB eine Transformation in Farbe eine gleichzeitige Transformation in Flavor voraussetzt, um die 
Invarianz des Kondensates zu bewahren. Deswegen nannten die Entdecker [18] diese Phase audi 
color- flavor- locked (CFL)-Phase. 

Bei hohen Dichten und niedrigcn Tcmperaturen ist die CFL-Phase der wahre Grundzustand 
von Quarkmaterie, da sie den hochsten Druck hat und weil alle Quarks gegappt sind. Im Gcgcnsatz 
zur 2SC-Phase mit ihren Antitriplettgaps A existicrcn in der CFL-Phase neben ihren Antitriplctt- 
gaps A(g g) auch symmetrische Sextettgaps A( 6j6 ). In den Quasiteilchenspektra treten sowohl der 
Singlettgap Ai, als auch der Oktettgap A2 auf. Wenn man den kleinen repulsiven Sextettgap 
vernachlassigt, dann findct man heraus, daB Ai = 2A 2 = 2A ist. In der CFL-Phase existieren 
keine gaplosen Quarks. 

Die Beitrage der Quasitcilchcn zu alien Transport- und vielen thermodynamischen GroBen ist 
bei niedrigcn Tcmperaturen durch den exponentiell kleinen Faktor exp (— A/T) unterdruckt [9]. 
Der Einflufi der Gluonen ist vernachlassigbar, da in der CFL-Phase alle Gluonen gegappt sind. 
Im Gegensatz zur 2SC-Phase ist die CFL-Phase superfluid, weil die U(1)b Baryonzahlsymmctrie 
gebrochen ist. Aber die CFL-Phase besitzt eine ungebrochene [E/(l) e m] Eichsymmetrie und ist 
daher, wie die 2SC-Phase, kein elektromagnetischer Supraleiter. Daher drangt die CFL-Phase 
auch kein Magnetfeld aus ihrem Inneren heraus. Die CFL-Phase ist ein Q-Isolator, da alle Quarks 
gegappt sind und keine elektrischc Ladung verbleibt, wie das bei den ungegappten blauen up 
Quarks in der 2SC-Phase der Fall ist. Daher ist die CFL-Phase elektrisch ladungsneutral. Bei 
T = sind in der CFL-Phase keine Elektronen vorhandcn [20]. Bei niedrigcn Tcmperaturen 
wird die elektrische Leitfahigkcit in der CFL-Phase durch thermisch angeregte Elektronen und 
Positroncn gcwahrleistet. Ein Phascniibcrgang zu eincr CFL-Phase mit eincm Mcsonkondcnsat 
ist moglich, falls M s > M^A 2 / 3 [26-30], wobei M u und M s die up bzw. strange Quarkmasse 
sind. 

Die 2SC- und die CFL-Phase sind farbsupraleitende Phasen, deren Kondensate einen Spin 
J = haben. Die Bildung von Quark-Cooper-Paaren bei J — fur nur einen Flavor ist wegen des 
Pauliprinzips nicht erlaubt. Aber dies ist moglich bei J = 1, d.h. bei Spin-Eins-Farbsupraleitung. 
Spin-Eins-Farbsupraleitung ist jedoch viel schwacher als Spin-Null-Farbsupralcitung. Wahrend 
der Gap in den Phasen der Spin-Null-Farbsupraleitung ca. 100 MeV groB ist, betragt er in 
den Phasen der Spin-Eins-Farbsupraleitung nur ca. 100 keV. Solch ein kleiner Gap hat keine 
groficn Einfliisse auf Transport- und viele thermodynamischc Eigenschaftcn von Quarkmaterie [9] . 
Spin-Eins-Farbsupraleitung ist weniger energetisch bevorzugt als Spin-Null-Farbsupralcitung, da 
die lctztcrc einen hohcrcn Druck besitzt. Daher ist es auch nicht zu erwarten, dafi Spin-Eins- 
Farbsupraleitung im Phascndiagramm von neutraler Quarkmaterie auftritt. Jedoch konnte Spin- 
Eins-Farbsupraleitung bevorzugt werden, wenn der Unterschied zwischen den Fermioberflachen 
bei Spin-Null-Farbsupraleitung zu groB ist. 

Die wichtigsten Phasen der Spin-Eins-Farbsupraleitung sind die A-Phase, die color-spin-locked 
(CSL)-Phase, die polare Phase und die planare Phase. In Spin-Eins-Farbsupralcitcrn kann, im 
Gegensatz zu Spin-Null-Farbsupraleitcrn, ein elektromagnetischer MciBncreffekt auftretcn. Dies 
ist z.B. der Fall in der CSL-Phase. Die am meisten energetisch bevorzugte Phase der Spin-Eins- 
Farbsupraleitung ist die transverse CSL-Phase, da sie den hochsten Druck hat [39]. 

Sterne entstehen in interstellaren Gaswolkcn. Sie sind nuklcarc Fusionsfabrikcn, da in ihnen 
lcichtcrc Kerne zu schwereren Kernen vcrbrannt werden. Nachdem alle Fusionsreaktionen abge- 
schlosscn sind, werden Sterne entweder zu WeiBen Zwergen, Neutronensternen oder Schwarzen 
Lochern. Ein Neutronenstern entsteht, wenn ein Roter Ubcrriese seine Fusionsreaktionen beendet 
hat. Dann kollabiert dieser, weil der Innendruck fehlt und daher die Gravitation uberwiegt, bis 
schliefilich der Druck von entarteten Neutronen den Kollaps zum Erliegen bringt. Die aufieren 
Schichten fallen auf den entstandenen Protoneutronenstern, prallen zuriick und erzeugen dadurch 
eine Supernova vom Typ II. Durch den Kollaps wird eine riesige Anzahl von Neutrinos produziert, 
die im Protoneutronenstern zunachst gefangen sind. Dies wird als Neutrinotrapping bezeichnet. 
Der zunachst heiBe Protoneutronenstern kiihlt durch Neutrinoemission und spater durch Photoe- 
mission ab. Neutronensterne werden als Pulsare detekticrt. Dicsc senden wegen des sogenannten 
Leuchtturmcffckts sehr prazise Radiopulse aus und besitzen ein starkes Magnetfeld. Ihre Rotati- 
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onsperiode liegt im Millisekunden- bis Sekundcnbereich. 

Neutronensterne bestehen nicht etwa, wie der Name vermuten lafit, nur aus Neutronen. Sie 
bestehen jedoch hauptsachlich aus Neutronen. Es gibt in Neutronensternen bei verschiedenen 
Massendichten unterschiedliche Grundzustande. Die Materie in Neutronensternen ist im Gleich- 
gewicht beziiglich der schwachen Wechselwirkung, im sogenannten /3-Gleichgewicht. Zudem sind 
Neutronensterne ladungsneutral, d.h. sie tragen keine elektrische Ladung und keine Farbladung, 
denn ansonsten wurden sie instabil sein. 

Neutronensterne besitzen eine Atmosphare, die aus Elektronen, Kernen und Atomen besteht. 
Falls der Neutronenstern eine Temperatur T > 100 eV besitzt, dann ist seine Oberfiachenschicht 
fliissig. Sie besteht aus Kernen und Elektronen. Bei einer Massendichte p ~ 10 4 g/cm 3 beginnt die 
auBere Kruste von Neutronensternen, die aus Elektronen und Kernen besteht, wobei letztere in 
einem body-centered cubic (bcc) lattice angeordnet sind, wodurch ein giinstigerer Energiezustand 
eingenommen wird. Mit zunchmcndcr Massendichte werden die Kerne immer neutronenreicher, 
bis bei p ~ 4.3 • 10 11 g/cm 3 die sogenannte Neutrondripline erreicht wird, bei der die Neutronen 
beginnen, sich von den Kernen zu losen. Das Ende der auBeren Kruste ist damit erreicht. Die 
innere Kruste von Neutronensternen besteht aus Kernen, Neutronen und Elektronen. Bei hoheren 
Massendichten zerfallen die Kerne in ihre Bestandteile, Protonen und Neutronen. Diese werden 
superfluid. Auch Myonen treten in diesen Neutronensternschichten auf. In welcher Phase sich die 
Kerne von Neutronensternen behnden, hangt von deren Zentraldichte ab. Neutronensterne mit 
geringerer Zentraldichte bestehen in ihrem Inneren nur aus Protonen, Neutronen, Elektronen und 
Myonen. Neutronensterne mit grofierer Zentraldichte konnten in ihrem Inneren aus Pionen und 
Hyperonen bestehen, wahrend sich bei Neutronensternen mit riesigen Zcntraldichtcn die Proto- 
nen, Neutronen und Hyperonen in ihre Bestandteile aufldsen: Quarks, die in den Kernen von 
Neutronensternen wahrscheinlich farbsupraleitend sind. 

Bei grofien strange Quarkmassen kann neutrale Quarkmaterie bestehend aus up und down 
Quarks im /3-Gleichgewicht einen eher ungewohnlichen Grundzustand einnehmen, den man gaplose 
2SC (g2SC)-Phase nennt. Die Symmetrie im g2SC-Grundzustand ist dieselbe wie in der regularen 
2SC-Phase. Jedoch weist das Quasiteilchenspektrum in der g2SC-Phase eine gapfreie Mode auf. 
Gapfreie Moden entstehen, wenn der Unterschied zwischen den Fermiimpulsen von sich paarenden 
Quarks grofi wird. Deswegen existiert auch eine gaplose CFL (gCFL)-Phase. Jedoch weisen sowohl 
die gaplosen, als auch die regularen farbsupraleitenden Quarkphasen bei niedrigen Temperaturen 
chromomagnetische Instabilitaten auf [93]. 

Um nun vorhersagen zu konnen, welche Phasen in der Natur energetisch bevorzugt werden, 
ist es notwendig, das Phasendiagramm von neutraler Quarkmaterie zu erstellen, in das diejenigen 
Phasen eingetragen werden, die bei vorgegebenen Temperaturen und quarkchemischen Potentia- 
len vorherrschen. In diesem Zusammenhang ist es wichtig, neutrale Quarkmaterie zu betrachten, 
da ansonsten repulsive Coulombkrafte auftreten wurden, die die Neutronensterne zur Explosion 
zwingen wurden. Zudem muB Quarkmaterie auch farbladungsneutral sein, weil farbgeladene Ob- 
jekte noch niemals in der Natur gesehen worden sind und auch farbgeladene Neutronensterne 
instabil waren. Zudem befindet sich die Materie in Neutronensternen im /3-Gleichgewicht. Das 
Phasendiagramm von neutraler Quarkmaterie, bestehend aus up, down und strange Quarks im 
/3-Gleichgewicht, habe ich im Kapitel[5]erstellt. 

Im Abschnitt 12.11 habe ich masselose neutrale Quarkmaterie bestehend aus up, down und 
strange Quarks bei hohen Dichten untersucht. Die Effekte der strange Quarkmasse m s sind durch 
eine Verschiebung im chemischen Potential der strange Quarks beriicksichtigt worden, p s — > p s — 
m 2 s / (2p). Diese Verschiebung spiegelt die Reduktion der Fermiimpulse der strange Quarks wegen 
der strange Quarkmasse wider. Solch eine Naherung ist sicherlich zuverlassig bei kleinen Werten 
der strange Quarkmasse. In Ref. [113], wo die strange Quarkmasse richtig miteinbezogen wurde, 
bestatige ich, dafi sie auch qualitativ korrekte Ergebnisse bei groBen strange Quarkmassen liefert. 
Ich habe eine sehr mannigfaltige Phasenstruktur durch Variieren der strange Quarkmasse, des 
quarkchemischen Potentials und der Temperatur erhalten. 

Bei T — gibt es zwei Moglichkeiten fur den Grundzustand von Quarkmaterie: die CFL- und 
die gCFL- Phase [90,91]. Ich bestatige, daB das farbchemische Potential dafiir verantwortlich ist, 
daB die CFL-Phase bei einem kritischen Wert der strange Quarkmasse in die gCFL-Phase iibergeht. 
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Zusammenfassung 



Da ich eincn ncunparamctrischcn Ansatz fiir die Gapmatrix benutze, bin ich in der Lage, auch die 
Effckte der repulsiven Sextettgaps zu untersuchen. Wie zu erwarten war, modifizieren diese die 
Quasiteilchendispersionsrelationen nur geringfiigig. Jedoch ist diese Untersuchung wichtig gewesen, 
um zu kontrollieren, ob die Sextettgaps vielleicht ungeahnte Auswirkungen auf den Ubergang von 
der CFL- zur gCFL-Phase haben. 

Im Abschnitt 12.11 habe ich auch die Temperaturabhangigkeit der Gapparameter und der Qua- 
siteilchenspektra untersucht. Ich habe gezeigt, dafi viele unterschiedliche Phasen existieren, die in 
Ref. [Ill] vorhergesagt wurden. Jedoch erhalte ich, im Gegensatz zu den Autoren von Ref. [Ill, 
112], die die dSC-Phase erhalten, in der alle down Quarks gepaart sind, die uSC-Phase, in der alle 
up Quarks gepaart sind. In Ref. [113] habe ich mein Resultat bestatigt, dafi die uSC-Phase und 
nicht die dSC-Phase im Phasendiagramm von neutraler Quarkmaterie auftritt. Das Auftreten der 
dSC-Phase hangt stark von den Modellparametern, genauer gesagt vom Cutoff, ab. Wahlt man 
namlich einen grofieren Cutoff als in Ref. [Ill, 112], dann erhalt man in der Tat die dSC-Phase, 
und die uSC-Phase kommt im Phasendiagramm von neutraler Quarkmaterie im /3-Gleichgewicht 
nicht mehr vor. 

Das wichtigste Resultat in Abschnitt 12.11 ist das Phasendiagramm von masseloser neutraler 
Quarkmaterie bestehend aus up, down und strange Quarks im /3-Gleichgewicht, siehe Abb. 12.111 
In dieser Abbildung sind alle symmetriebezogenen Phasenubergange mit durchgezogenen Linicn 
gekennzcichnet. Die gestrichelten Linien kennzeichnen die Ubergange von regularen zu gaplosen 
Phasen. Dies sind jedoch keine echten Phasenubergange, sondern weiche Ubergange. Bei T = 
existiert ein Isolator-Metall-Phasenubergang zwischen der CFL- und der gCFL-Phase [90,91]. Bei 
hohcren Temperaturen ist ein ahnlicher Isolator-Metall-Phasenubergang zwischen der CFL- und 
der metallischen CFL (mCFL)-Phase vorhanden. Dieser ist in Abb. 12.111 durch die gepunkteten 
Linien gekennzeichnet. 

Im Abschnitt 12.21 habe ich das Phasendiagramm von neutraler Quarkmaterie mit Hilfe des 
Nambu-Jona-Lasinio (NJL)-Modells aus Ref. [136] untersucht. Wie in der vorangegangenen Un- 
tersuchung benutze ich eine Mean-Field-Naherung in der Analyse. Im Gegensatz zu Abschnitt 12. II 
behandele ich jedoch in Abschnitt 12.21 die Quarkmassen selbstkonsistent als dynamisch generierte 
GroBen. Die wichtigsten Resultate sind in Abb. 12.121 und Abb. I2.!3"l zusammcngcfafit. 

Beim Vergleichen des Phasendiagramms fiir masselose Quarks, siehe Abb. l2.lT1 rechts. mit dem 
Phasendiagramm in Abb. 12.121 bemerkt man einige wichtige Unterschiede. Es fallt auf, dafi die 
selbstkonsistente Behandlung von Quarkmaterie das Auftreten der unterschiedlichen Quarkphasen 
stark beeinflufit [139]. Zudem bestimmen aber auch die Phasen ihrerseits wiederum die Grofie der 
Quarkmassen. Dies fuhrt sehr oft zu Phaseniibergangen erster Ordnung. 

Viele Unterschiede zu den Resultaten in den Abschnitten 12.11 und 12.21 sind auf das Verwenden 
von unterschiedlichen Modellparametern zuriickzufuhren. Der Wert der Diquarkkopplungskon- 
stante Gd — jGs in Abschnitt 12.21 ist schwacher als jener in Abschnitt 12.11 was man sehr leicht 
sehen kann, indem man die Werte der Gaps bei T = und /i = 500 MeV vergleicht. Das sind 
A 500) w 140 MeV im Abschnitt O und A 5 00) w 76 MeV in Abschnitt O Im Vergleich hierzu ist 
der Wert der Diquarkkopplungskonstante in Ref. [112] noch viel schwacher, wo Ag 500 "* » 20 MeV 
ist. Im Fall von starker Kopplung im Abschnitt 12.21 mit Gd = Gs, wodurch Aq 500 -* sa 120 MeV 
ist, ist Aq 500 "* immer noch schwacher als im Abschnitt 12.11 In diesem Fall unterscheiden sich die 
Resultate hauptsachlich dadurch, dafi die Quarkmassen unterschiedlich behandelt werden. 

Wegen der niedrigen Diquarkkopplung geschehen die Cooperinstabilitaten in Abb. 12.121 bei 
grofieren Werten des quarkchemischen Potentials als in Abb. l2.lT1 rechts. Das kann man am besten 
am Ubergang zur (g)CFL-Phase erkennen. Eine andere Konsequenz der schwacheren Wechsel- 
wirkung ist die Moglichkeit eines thermal verursachten Anstiegs des (g)2SC-Cooperpaarens bei 
mittleren Werten des quarkchemischen Potentials. Diese Art des Anstiegs wurde in Ref. [84, 85] 
im Detail untersucht. Mit denselben Argumenten kann man sofort sagen, wie sich das Phasendia- 
gramm in Abb. 12.121 durch die Zu- oder Abnahme der Diquarkkopplungskonstante andern solltc. 

Insbcsondcrc solltc durch Zunahme (Abnahme) der Diquarkkopplungskonstante die Region der 
(g)2SC-Phase bei mittleren Werten des quarkchemischen Potentials expandieren (schrumpfen). Die 
Regionen, die durch die anderen (d.h. uSC- und CFL-) Phasen abgcdcckt werden, sollten qualitativ 
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dieselbe Form haben, sich jedoch zu niedrigeren (hoheren) Werten des quarkchemischen Potentials 
und zu h5heren (niedrigeren) Temperaturwerten verschieben. 

Das NJL-Modell, das ich in Abschnitt 12.21 verwende, ist definiert durch einen Parametersatz, 
dessen Werte so angepafit wurden, dafi sie wichtige QCD-Eigenschaften im Vakuum reproduzieren. 
(Das gleiche Modell wurde auch in Ref. [133] verwendet.) Es ist daher zu erwarten, dafi es ein gutes 
effektives Modell der QCD ist, das die beiden wichtigen Merkmale, chirale und farbsupraleitende 
Dynamik, abdeckt. Auch der relativ niedrige Wert des Cutoffs sollte nicht notwendigerweise als 
ein schlechtes Merkmal des Modells angesehen werden. Vielmehr konnte das eine Nachahmung der 
natiirlichen Eigenschaft der vollstandigen Theorie sein, in der die Kopplungsstarke von Wechsel- 
wirkungen bei grofien fmpulsen gering ist. 

In diesem Zusammenhang ist es wichtig darauf hinzuweisen, dafi die Naherung in Ref. [112] 
beziiglich des Cutoffs im NJL-Modell im Vergleich zu meiner sehr unterschiedlich ist. Dort wird 
namlich behauptet, dafi ein grofier Wert fur den Cutoff vorteilhaft sei, um Resultate zu erhalten, 
die unempfindlich auf eine spezielle Wahl des Cutoffs sind. Ich finde jedoch kein physikalisches 
Argument, das diese Voraussetzung unterstutzt. Vielmehr bestehe ich darauf, ein effektives Modell 
zu verwenden, das die QCD-Eigenschaften bei verschwindendem quarkchemischen Potential gut 
beschreibt. Natiirlich will ich damit nicht behaupten, dafi eine Extrapolation des Modells zu grofien 
Dichten damit zu rechtfertigen ist. Jedoch scheint mir das momentan die beste Losung zu sein. 

Die Resultate aus Abschnitt 12.21 konnten wichtig fur das physikalische Verstandnis von (hy- 
briden) Neutronensternen mit Quarkkern sein, in denen die Dclcptonisierung stattgefunden hat. 
Um ein Phasendiagramm zu erhalten, das auf Protoncutroncnsterne angewandt werden kann, mufi 
man den Effekt von Neutrinotrapping beriicksichtigen. 

Im Abschnitt 12.31 habe ich den Effekt von Neutrinotrapping auf das Phasendiagramm von neu- 
traler Quarkmaterie bestehend aus up, down und strange Quarks im /3-Gleichgewicht mit Hilfc 
des NJL-Modells aus Ref. [136] untersucht. Die Resultate habe ich in der Mean-Field-Naherung 
erhalten. Die Quarkmassen sind als dynamisch generierte Grofien behandelt worden. Die Gesamt- 
struktur der Phasen wird als dreidimensionales Diagramm im Raum von Temperatur T, quark- 
chemischem Potential fi und dem chemischen Potential der Leptonen \ii, e in Abb. l2.2Ul dargestellt. 
Dieses Phasendiagramm wird zudem anhand von samtlichen zweidimensionalen Teildiagrammen 
desselben detailliert diskutiert, siehe Abb. |272"H l2~22l l2~26l und [27271 

Sowohl anhand von einfachen modellunabhangigen Argumenten, als auch anhand von detail- 
lierten numerischen Bcrcchnungcn mit Hilfe des NJL-Modells habe ich herausgefunden, dafi Neu- 
trinotrapping das Cooper-Paaren in der 2SC-Phase begiinstigt und die CFL-Phase unterdriickt. 
Hauptsachlich kommt das durch das Erfullen der elektrischen Neutralitatsbcdingung zustande. In 
Quarkmaterie, die nur aus up und down Quarks im /3-Gleichgewicht besteht, hilft das chcmischc 
Potential der Leptonen zusatzliche Elektronen zur Verfugung zu stellen, ohne einen grofien 
Unterschied zwischen den Fermiimpulsen zu verursachen. Durch Verringern der Fermiimpulsun- 
terschiede wird das Cooperpaaren natiirlich starker. Dies ist eine vollig andere Situation als in der 
CFL-Phase, in der die Quarkmaterie in Abwesenheit von Elektronen neutral ist. Zusatzliche Elek- 
tronen wurden, wegen des grofien chemischen Potentials der Leptonen, die CFL-Phase belasten. 

Wendet man dieses Resultat auf Protoneutronensterne an, die farbsupraleitende Quarkkerne 
besitzen, dann bedeutet das, dafi die CFL-Phase wohl kaum in Protoneutronensternen auftritt 
bevor die Deleptonisierung komplett ist. Daher ist die 2SC-Phase der wahrscheinlichste Grundzu- 
stand in Protoneutronensternen mit farbsupraleitenden Quarkkernen. Zudem habe ich festgestellt, 
dafi Protoneutronensterne kaum strange Quarks enthalten. 

Die Autoren von Ref. [83] behaupten, dafi die 2SC-Phase in kompakten Sternen nicht vor- 
kommt und untermauern ihre Behauptung mit Hilfe einer einfachen modellunabhangigen Berech- 
nung. Die Phasendiagramme, die ich mit Hilfe des genaueren NJL-Modells erstellt habe und in 
meiner Dissertation zeige, sind der beste Beweis dafiir, dafi die 2SC-Phase wirklich in kompakten 
Sternen vorkommen kann. Falls kompakte Sterne farbsupraleitende Kerne besitzen, dann kame die 
2SC-Phase auf jeden Fall in Protoneutronensternen und in kalten Neutronensternen mit grofier 
Diquarkkopplungskonstante vor. Fur kleine oder mittlere Werte der Diquarkkopplungskonstan- 
te und bei verschwindender Temperatur, fur die die cinfache modcllunabhangigc Bcrcchnung in 
Ref. [83] giiltig ist, scheint es jedoch richtig zu sein, dafi die 2SC-Phase nicht auftritt. Wie dem 
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auch sei, man sollte mit der Aussage, daB eine Phase nicht auftritt, sehr vorsichtig sein. 

Nachdem die Deleptonisierung stattgefunden hat, ist es moglich, dafi der Grundzustand im 
Inncrcn von farbsupraleitenden Neutronensternen die CFL-Phase ist. Diese Phase besitzt eine 
grofie Anzahl von strange Quarks. Deswegen sollte eine grofie Produktion von strange Quarks 
sofort nach der Deleptonisierung einsetzen. Dieser Vorgang konnte beobachtbar sein. 

Im Kapitel [3] fasse ich meine Resultate zusammen und diskutiere sie. Zudem zeige ich offenge- 
bliebene Fragen auf und gebe einen Ausblick auf mogliche zukiinftige Untersuchungen. 

Im Anhang befinden sich niitzliche Formeln und Definitionen, die in meiner Dissertation ver- 
wendet werden. 
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